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A  MAP  FOR  IMPACT  INDUCED  DUCTILE  FAILURE 


Ian  M.  Fyfe 

University  of  Washington,  Seattle,  WA  (18195  U.S.A. 

Because  ductile  failure  is  the  nucleation,  growth  and  coalesence  of  voids  or  microcracks, 
the  paper  focuses  on  developing  a  failure  map  in  which  experimental  data  and  theories 
related  to  plastic  instabilites  and  failure  can  be  presented  in  terms  of  void  nucleation 
concepts. 

For  failure  maps  in  general  the  choice  of  variables  is  not  obvious,  but  depends  to  n 
large  extent  on  why  the  map  is  being  used  in  the  first  place.  Some  useful  guidelines 
have  been  presented  by  Ashby  (1977),  for  use  in  earlier  failure  map  construction,  where 
the  emphasis  was  on  creating  fracture  maps  showing  the  boundaries  between  different 
mechanisms  of  fracture.  In  the  area  of  high  strain  rate  behavior  Lindholm  (1974),  in 
a  review  of  dynamic  testing  methods,  was  more  concerned  with  constitutive  model¬ 
ing,  and  so  used  effective  stress  and  effective  strain,  with  each  curve,  designated  by 
a  constant  temperature  strain  rate  parameter.  However,  void  growth  and  spallation 
strongly  suggest  that  failure  maps  for  impact  loading  should  also  include  the  mean 
stress  variations.  The  variables  used  in  this  paper  are  those  which  describe  nucleation, 
but  also  cover  the  high  strain  rate  conditions  associated  with  shear  bands  and  spalla¬ 
tion.  A  case  is  made  for  using  a  normalized  mean  stress,  effective  plastic  strain  and 
temperature. 

With  this  choice  of  variables  it  is  possible  to  include  results  from  three  experimental 
configurations.  Considered  are  the  high  mean  stress  case  associated  with  the  plate 
impact  test,  the  intermediate  condition  of  simple  tension,  and  the  zero  mean  case 
which  results  from  the  torsion  configuration. 

As  one  of  the  main  objectives  of  impact  testing  is  to  determine  the  strain  rate  sensitiv¬ 
ity  of  the  process,  and  the  influence  of  the  different  micromechanic  mechanisms,  both 
dynamic  and  quasi-static  data  is  used.  Use  is  also  made  of  models  developed  from 
both  a  microscopic  and  a  macroscopic  viewpoint.  Although,  this  paper  deals  mainly 
with  the  macro-level  of  damage  initiation  as  characterized  by  localization  and  the  loss 
of  structural  integrity. 

With  nucleation  occurring  at  second  phase  particles  or  other  inclusions,  the  process 
then  depends  on  particle  size.  For  this  reason  it  has  been  determined  that  nucleation 
can  occur  almost  immediately  after  yield,  and  can  continue  throughout  the  deforma¬ 
tion  process.  However,  at  some  stage  of  this  deformation  it  has  been  observed  that 
a  proliferation  point  is  reached  where  the  nucleation  increases  dramatically,  and  it 
might  be  expected  that  this  proliferation  occurs  when  the  macroscopic  condition  of 
localization  and  necking  occurs.  Thig  possibility  is  discussed  in  relation  to  the  two 
materials  used  in  this  study. 

It  will  be  shown  that  the  Hill  (1952)  theory  for  localized  necking  in  a  thin  sheet  is 
quite  compatible  with  plastic  instabilities  under  high  strain  rate  loading  conditions, 
and  that  the  formation  of  instabilities  is  essentially  independent  of  the  mean  stress.  In 
contrast  ductile  failure  was  found  to  be  highly  dependent  on  the  mean  stress,  and  the 
dislocation  based  nucleation  model  of  Goods  and  Brown  (1979)  was  also  found  to  be 
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applicable  under  both  quasi-static  and  the  high  strain  rate  loading  conditions.  These 
results  suggest  that  plustic  instabilities  may  not  be  triggered  by  void  nucleation,  even 
though  ductile  failure  is  in  many  instances  dependent  on  this  mechanism,  at  least  in 
the  finnl  phase  of  the  failure  process.  These  results  as  they  apply  to  VAR  4340  steel 
are  presented  in  Figure  I. 


Figure  1  Effective  plastic  strain  as  a  function  of  normalized  tensile  stress 
dynamic  -}  and  static  (— O— *)  loading  conditions 
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Y  1,  Bai 

Institute  Of  Mechanics ,  Chinese  Academy  of  Sciences, 

Beijing,  Chinn  100080 

ABSTKACT 

This  pnpm'  reviews  the  recent  advances  in  the  study  of  evolution  of  shear 
deformation  into  local  i  zed  band- 1  ike  si  ruet.ure  in  thermo-v  i  sco-pl  astir 
materials,  Particularly,  the  progress  made  in  the  impact  laboratory  of 
the  Institute  of  Mechanics  is  summarized.  All  the  experimental  observa¬ 
tions  and  analytical  and  numerical  results  show  that  the  evolution  of 
shear  from  homogeneous  field  into  so-called  adiabatic  shear  band  is  a 
multi-stage  process.  The  band-like  structure  is  asymptotically  stable. 
So,  to  some  extend,  no  matter  how  the  disturbances  or  microscopical 
structures  are,  the  eventual  lea tores  of  the  shear  bands  look  roughly 
alike  in  nature,  although  complicated  details  do  exist  in  early  and 
i n termed) ate  stages . 

In  dealing  with  the  concerned  problem,  a  we  I  I -posed  model  and  appropri¬ 
ately  simplified  governiiiii  egn.it  ions  are  i  lie  necessity.  A  one-dimen¬ 
sional  model  of  simple  shear  was  leviewed.  Km  t hermere ,  sealing  is  oi 
llloi.il  sigiiil  iiMiii'e  in  the  .inilviiis,  .  i  i  nee  vat  in, is  I  i  me  and  length  scales 
are  involved  in  various  stages  in  the  evolution.  It  was  pointed  out 
th.il  the  effective  I’r.iruli  I  number  appears  to  he  l  he  key  dimensionless 
parameter  in  approximations,  which  describe  various  stages  of  the  shear 
evolution.  In  particular,  'lie  implications  of  gua.s  i  t  a  t  i  e  and  adiabatic 
approximal  ions  are  di  iiciiss.  d . 

Instability  criteria  reuniting  finm  linear  peri urbu t i  mi  analysis  provide 
helpful  i  n  I  oriiiat  in,  i  about  the  ending  of  uni  bum  shea,  mg  and  growth  rate 
of  disturbances  at.  early  stage.  However,  i  ns  Labi  1  i  t.y  is  not.  the  synonym 
of  localization.  Therefore,  the  concept  and  implication  ot  localization 
are  reviewed, 


Experiments  show  that,  a  number  of  very  fine  shear  bands  appear,  when 
instability  has  just,  ni'cured,  although  macroscopic  uni  form  shearing  seems 
to  remain.  Following  this,  a  localized  shear  zone  emerges  and  develops 
in  the  centre  of  test  piece.  Afterwards,  it  drastically  predominates  the 
whole  field  of  the  shear  deformation.  Finally,  a  fairly  steady  shear 
band  forms.  Numerical  simulations  were  carried  out  with  different,  kinds 
of  disturbances  iri  thermo-visco-plastic  material  and  revealed  all  those 
details  of  the  process  very  dearly.  In  order  to  understand  the  under¬ 
lying  mechanism,  which  governs  the  transition  In  predominating  shear  band 
analytical  interpretation  was  provided.  It  seems  th.d  the  hand— like 


si  riicture  riwul  Is  I  row  Mir  l>.» I anee  of  the  input  plust  if  work  rate  and 
the  heat  diffusion.  This  median  ism  can  support,  a  quasi -steady  shear  band 
in  t.lir  thormo-visoo-plusl  ir  material  as  well.  A  close  form  quasi -steady 
solution  and  art  estimation  of  I  he  shear  band  width  are  qiven.  A  variety 
ol  observe!  I  ions  appears  to  hr  in  fairly  good  agreement,  with  the  suggested 
predict i on . 

It  is  well  known  that  the  so-called  adiabatic  shear  bands  usually  serve 
as  precusor  to  fracture.  But  how  the  eventual  shear  fracture  forms  and 
why  the  eventual  lracture  sui'iace  is  almost  always  parallel  to  shear 
force  remain  unclear.  It,  is  reported  that  the  cumulation  of  microdamage, 
which  accompanies  evolution  of  shear  deformation,  plays  crucially 
important  role  in  the  formation  of  the  shear  fracture.  Optical  and  TEM 
observations  reveal  some  microscopic  details  of  the  two  processes  - 
shear  localization  and  shear  livicfure,  Perhaps,  crystallographic  slip  on 
dominate  slip  system  within  grains  should  be  responsible  for  the  formation 
of  microscopic  shear  bands.  Whereas,  large  strains  in  localized  shear 
zone  and  the  stress  concent  rai ions  resulting  from  dislocation  pile  ups 
at  grain  boundaries  and  phase  boundaries  initiate  mi crooracks . 

T  t.  is  suggested  that  the  theory  concerning  shear  band  induced  fracture 
is  badly  wanted,  especially  in  engineering  practice.  So  far,  to  the 
author's  knowledge,  there  is  nearly  no  such  a  kind  ol  theories  available. 
Consideration  of  the  cumulation  of  microdamage  and  its  variation  with 
the  shear  band  evolution  may  be  of  crucial  importance.  The  patterns  of 
shear  bands  under  combined  stresses  also  interest  engineers  in  industries. 
It  is  hoped  that  more  practical  models  and  theories  can  be  developed 
tor  this  east.  Finally,  the  role  of  material  transformations,  such  as 
phase  transformation,  vanishing  strain  hardening,  He.,  iri  I  lie  fully 
developed  shear  band,  needs  to  tie  clarified. 


EXPERIMENTAL  METHODS  OF  MATERIAL  CHARACTERIZATION 


J.L.  LATAILLADE 

Laboraloire  de  M&anique  Physique,  Unitd  de  Recherche  Associde  au  CNRS  n®  867 
351  Cours  de  la  Liberation,  33405  TALENCB  Ccdex  -  France 


Summary 

The  operating  conditions  -  more  and  more  severe  -  of  the  mechanical  structures,  as  well  as 
the  security  rules  being  Increasingly  constraining  •  especially  in  nuclear  industries  explain  the 
subtantial  development  of  the  experimental  research  is  the  field  of  high  strain  rate  testing.  During 
the  fifteen  past  years  the  experimental  methods  and  the  related  testing  techniques  for  materials 
becoming  more  and  more  tough  and  strong  have  received  a  particular  attention  of  the  researchers. 
The  general  purpose  remains,  of  course,  the  materials  characterization  by  means  of  reduced  scales 
configurations;  but  a  new  approach  -  due  to  the  introduction  of  microcomputer  and  hence  of  the 
appearance  of  hybrid  methods  -  are  noticeable. 

The  experimental  techniques  based  on  the  well  known  principle  of  the  Iiopkinson-Kolsky 
bare  are  the  most  popular.  Some  modifications  and  improvements  are  reported  in  this  conference. 
They  allow  a  better  control  of  the  loading  paths  and  consequently  a  more  appropriate  modelization 
of  the  viscoplastic  behaviour  of  materials  such  as  metallic  alloys  or  a  finer  approach  of  the  damage 
processes  of  advanced  materials.  In  order  to  extend  the  range  of  the  strain  rates  people  use  methods 
based  on  the  stress  waves  propagation.  The  techniques  allow  them  to  generate  biaxial  stresses,  and 
the  lagrangian  analysis  helped  by  high  performance  guages  lead  to  a  modem  tool  for  the 
constitutive  relationship  at  very  high  strain  rates. 

On  an  other  hand,  the  increasing  variety  of  structural  advanced  materials  (composites, 
thcrmomechanical,  ceramics,  engineering  polymers),  or  the  new  requirements  from  the  civil 
engineering  (concrete,  rocks)  have  suggested  new  methods,  made  easier,  owing  to  new 
instruments  for  the  measurements,  or  the  observation  of  phenomena. 

From  this  point  of  view  the  optical  methods  (laser,  high  speed  cinematography)  become 
more  and  more  employed,  since  they  afford  mechanical  parameters  of  a  great  importance  like 
strains,  COD,  stress  intensity  factors,  etc  ... 

Because  high  strain  rate  testing  generally  involves  adiabatic  processes  it  can  be  interesting 
or  necessary  to  measure  the  temperature  history  of  sample,  in  order  to  set  up  a  thermomechanical 
approach,  and  to  appreciate  its  influence  on  the  micros tructural  modifications  ;  the  best  example  is 
that  one  of  the  shear  banding.  Due  to  the  very  short  rise  times  the  only  experimental  methods  is  the 
infrared  thermography. 

The  advancements  of  the  physics,  microelectronics  opto  electronics  have  led  to  new 
transducers  very  sensitive  exhibiting  very  small  response  times  :  some  of  their  applications  are 
mentioned  in  this  conference. 
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THE  DEVELOPMENT  OF  CONSTITUTIVE  RELATIONSHIPS  FOR  MATERIAL 
BEHAVIOUR  AT  HIGH  RATES  OF  STRAIN 

J.  Harding 

Department  of  Engineering  Science 
University  of  Oxford 
Oxford,  U.K. 

ABSTRACT 

The  increasing  use  of  computer  codes  for  modelling  the  mechanical  response 
of  engineering  structures  under  impact  loading  has  highlighted  the  need  for 
appropriate  forms  of  constitutive  relation,  or  mechanical  equation  of  state, 
capable  of  describing  material  behaviour  over  a  wide  range  of  strain  rate  and 
temperature.  Attempts  to  derive  such  relationships  have  varied  from,  at  one 
extreme,  the  purely  empirical,  effectively  a  curve  fitting  exercise  to  avail¬ 
able  experimental  data  for  the  given  material  to,  at  the  other,  the  develop¬ 
ment  of  purely  theoretical  expressions  based  on  the  micromechanical  processes 
governing  plastic  flow  at  the  atomic  level  under  the  given  loading  condit¬ 
ions. 

Ideally  these  two  contrasting  approaches  should  come  together  in  a  single 
constitutive  relationship  which  both  describes,  with  good  accuracy,  the 
actual  macroscopic  behaviour  of  the  given  material  and  can  also  be  related 
to,  and  hence  derived  from,  the  physical  processes  which  control  plastic  de¬ 
formation.  However,  while  this  ultimate  goal  is  still  some  way  from  being 
achieved  it  remains  necessary  to  consider  both  approaches  when  seeking  to  ob¬ 
tain  the  most  suitable  form  for  the  constitutive  relation.  As  a  result  sev¬ 
eral  semi-empirical  expressions  have  been  proposed,  the  functional  form  of 
which  is  based  on  theoretical  considerations  but  involving  empirical  const¬ 
ants  which  cannot  be  related  directly  to  precise  physical  processes  and  which 
have  to  be  determined  for  the  particular  material  being  studied. 

The  present  review  is  mainly  concerned  with  this  semi-empirical  approach. 
Attempts  to  model  material  behaviour  in  terms  of  the  micromechanisms  of  plas¬ 
tic  flow  will  be  the  concern  of  a  subsequent  review  by  Klepsczko  while  a  dis¬ 
cussion  of  the  use  of  computer  codes  in  modelling  the  actual  impact  response 
of  engineering  structures  will  be  the  subject  of  the  review  by  Corran  in  the 
parallel  session. 

Experimental  results  for  the  temperature  and  strain-rate  dependence  of  the 
flow  stress  at  different  rates  of  strain  are  briefly  reviewed  for  a  variety 
of  materials.  Functional  relationships  between  the  stress,  the  temperature 
and  the  strain  rate  are  presented  and  their  theoretical  basis  is  discussed. 
The  problems  involved  in  developing  a  full  equation  of  state,  which  includes 
some  measure  of  the  extent  of  deformation,  are  considered  in  the  light  of  the 
frequently  observed  dependence  of  the  flow  stress  on  the  previous  strain  rate 
history  and  the  need  for  some  internal  structure  sensitive  parameter  is  made 
clear.  Two  attempts  at  defining  the  structural  state  of  the  material  in 
terms  of  such  a  parameter,  either  the  rate-sensitivity  of  work-hardening 
(Klepaczko  and  Chiem,  1986)  or  the  mechanical  threshold  stress  (Follansbea  et 
al.  1985)  are  described  and  compared.  In  neither  case,  however,  has  an  ex¬ 
plicit  functional  relationship  yet  bean  developed  of  s  form  suitable  for  use 


in  computer  codes  for  the  numerical  modelling  of  the  impact  response  of  a 
structural  component. 

In  practice,  the  specific  constitutive  relationships  which,  although  less 
satisfactory  on  theoretical  grounds,  have  been  used  in  computer  codes  are,  in 
general,  of  two  types  -  either  the  empirical,  for  example  Johnson  and  Cook 
(1983),  or  the  serai-empirical,  for  example  Zerilli  and  Armstrong  (1987).  The 
results  obtained  when  these  two  materials  models  are  used  in  the  computer 
program  EPIC2  to  describe  the  structural  response  in  the  cylinder  impact  test 
are  compared  and  the  differences  are  seen  to  be  small.  It  is  concluded  that 
in  this  application  the  numerical  solution  is  relatively  insensitive  to  the 
details  of  the  material  model.  While  refinements  to  the  Zerilli-Arrastrong 
model  lead  to  a  second  order  improvement  in  the  agreement  between  the  exper¬ 
imental  results  and  the  numerical  analysis  the  corresponding  form  of  the  con¬ 
stitutive  relation  is  of  considerably  greater  complexity. 

In  general,  therefore,  the  use  of  relatively  simple  semi-empirical  const¬ 
itutive  relationships  would  seem  to  be  adequate  for  the  computer  modelling  of 
many  impact  problems.  However,  in  situations  where  the  overall  impact  re¬ 
sponse  of  the  given  component  may  be  critically  determined  by  the  deformation 
processes  and  complex  internal  microstructural  states  arising  in  very  local¬ 
ised  regions  then  the  details  of  the  material  model  may  become  much  more  imp¬ 
ortant.  Thus  in  studying  the  critical  conditions  for  the  initiation  of 
adiabatic  shear  in  aluminium  and  steel  Klepaczko  (1988)  found  it  necessary  to 
use  a  "fully  temperature  coupled"  form  of  constitutive  relation  which,  in  the 
case  of  steel,  involved  a  total  of  fifteen  different  material  constants.  In 
such  circumstances  it  might  prove  more  profitable,  it  is  suggested,  to  pursue 
further  the  development  of  material  models  based  on  the  concept  of  an  inter¬ 
nal  structure  sensitive  parameter  for  which  there  is  strong  theoretical  sup¬ 
port.  A  further  problem  arises,  however,  if  the  critical  condition  arises 
from  a  change  in  the  controlling  mode  of  deformation  from,  say,  constant  vol¬ 
ume  plasticity  to  the  growth  and  coalescence  of  voids  or  the  propagation  of 
cracks,  requiring  in  addition  to  the  appropriate  constitutive  relation  also 
the  inclusion  in  the  computer  code  of  some  form  of  failure  criterion.  This, 
however,  is  outside  the  scope  of  the  present  review. 
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This  review  addresses  the  modelling  of  dynamic  plastic 
deformation  in  solid s.  it  commences  with  a  survey  of  the 
techniques  employed  in  codes  that  are  available  either 
commer daily  or  in  the  public  domain.  Lagrangian,  Sulerian  and 
Arbitrary  Lagrangian  Sulerian  formulations  are  described  to  allow 
a  discussion  of  the  different  methods  of  problem  definition  and 
how  this  influences  the  limitations  to  the  modelling.  It  is 
found  that  the  Lagrangian  formulation  is  most  useful  in  problems 
involving  impact,  since  it  maintains  a  olear  distinction  of 
boundaries,  but  has  the  disadvantage  that  the  discretization, 
i.e.  the  finite  element  mesh  or  finite  difference  grid,  may 
become  very  distorted  leading  to  numerical  difficulties.  The 
Sulerian  formulation  avoids  this  problem  but  does  not  provide  a 
clear  model  of  boundaries  between  bodies  or  materials.  The 
arbitrary  Lagrangian  Sulerian  formulation  is  a  method  that  has 
not  been  much  used  for  solid  mechanics  problems,  although  it 
combines  some  of  the  advantages  of  the  two  methods. 

Time-stepping  through  a  solution  can  be  achieved  by  either 
explicit  or  implicit  methods.  Explicit  methods  have  the 
advantage  that  they  involve  no,  or  possibly  limited,  iteration 
and  hence  proceed  through  each  time-step  at  minimal  computational 
cost.  However  the  length  of  each  time-step  must  be  kept  very 
small  tc  avoid  instability  in  the  procedure  and  maintain  an 
accurate  modelling.  Implicit  methods  may  take  considerably 
larger  time-steps,  up  to  a  thousand  times  longer,  but  involve 
iterations  to  achieve  equilibrium  and  compatibility  at  each 
time-step.  Although  explicit  methods  are  generally  used  for 
impact  modelling  in  three  dimensions  because  of  their  relative 
efficiency,  implicit  methods  may  be  more  appropriate  for  simpler 
problems . 

Difficulties  with  the  modelling  of  continua  discussed  are 
the  control  of  hourglassing,  the  level  of  discretization  and  the 
importance  of  the  constitutive  model.  Hourglassing  is  a 
deformation  mode  that  involves  aero  strain  energy  for  the  finite 
difference  or  finite  element  in  use,  but  it  can  grow  in  an 
unstable  manner  to  dominate  the  observed  behaviour.  An  example 
is  discussed  and  methods  of  controlling  the  hourglassing  are 
reviewed.  It  is  found  that  the  hourglassing  can  occur  not  only 
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on  an  element  or  tonal  level/  but  in  oertain  problems  at  a 
structural  level.  Unfortunately  the  hourglass  modes  of 
deformation  may  be  modes  that  occur  due  to  the  dynamic  problem 
being  modelled.  Hence  it  is  necessary  to  control  them  without 
eliminating  them  completely. 

An  example  of  the  effects  of  inadequate  discretization  is 
given  showing  how  the  provision  of  too  few  zones  in  a  finite 
difference  modelling  leads  to  unrealistic  impact  behaviour  in 
predicting  the  contact  time  between  the  missile  and  the  target. 
Further  the  coarse  discretization  leads  to  more  hourglassing  and 
causes  an  unstable  development  of  deformation  away  from  the 
impact.  This  implies  that  although  there  may  be  little  interest 
in  studying  the  missile  behaviour,  as  muoh  care  must  be  taken  in 
its  modelling  as  in  the  target. 

A  limitation  of  simple  constitutive  relations  which  do  not 
include  rate  sensitivity  is  that  any  localization  of  the 
deformation  is  related  to  the  discretization.  Hence  the  width 
of  any  shear  bands  that  may  be  identified  during  the  modelling 
will  be  of  one  element’s  width  in  a  simple  finite  element 
description.  However  if  the  rate  sensitivity  is  included  the 
governing  equations  show  that  the  shear  band  will  have  a 
dimension  related  to  the  constitutive  model.  Although  most 
codes  do  not  have  any  method  for  modelling  such  a  concentration 
of  the  deformation,  a  discussion  of  the  various  approaches  to 
shearing  is  given,  based  on  solutions  to  a  more  ideallised 
one-dimensional  problem.  It  is  found  that  for  a  rate-sensitive 
constitutive  model  there  Is  indeed  a  natural  sizing  of  the  shear 
band,  a  result  that  is  apparent  from  both  finite  element  studies 
and  an  analytical  solution. 

Although  this  review  concentrates  on  limitations  to  the  use 
of  numerical  modelling,  it  is  concluded  that  with  care  successful 
modelling  of  dynamic  plastic  behaviour  is  possible.  The  major 
limitation  at  the  present  is  in  the  modelling  of  failure, 
especially  ductile  failure.  As  studies  of  shear  bands  have 
shown,  however,  this  is  particularly  difficult  to  model  although 
it  is,  paradoxically,  theoretically  easier  if  a  more  complicated 
constitutive  model  is  used.  Nevertheless  it  is  possible  to  use 
some  form  of  damage  mechanics  to  give  useful  estimates  of  duatile 
fracture  without  modelling  localization.  Finally  an  example  of 
an  impact  problem  of  industrial  relevance  is  shown  in  whioh  rate 
sensitivity  was  found  to  be  important  although  strain  hardening 
was  not  modelled.  In  this  case  the  object  is  to  predict  the 
final  deformed  shape  after  impact  for  which  dynamic  relaxation 
was  used  to  damp  out  the  residual  velocities  after  all  plastic 
deformation  had  ceased. 
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EXTENDED  ABSTRACT  /  '  f  '' r 

During  last  decades  a  vast  experimental  evidence  has  been  accumulated 
that  metals  deformed  plastically  at  different  conditions  (strain  rate, : 
temperature,  pressure)  undergo  a  complicated  microstructural  evolution. 
Recently,  microstructural  aspects  have  been  more  frequently  introduced 
into  constitutive  modeling,  however,  the  emerging  picture  is  still  not 
complete.  Such  approach,  based  on  materials  science,  is  very  promising  in 
formulation  a  framework  for  rational  constitutive  modeling  with  a  possi¬ 
bility  of  reliable  extrapolations  beyond  experimental  conditions  for 
which  the  constitutive  parameters  have  been  determined. 

This  discussion  presents  a  consistent  approach  to  microstructure  evolution 
via  the  constitutive  formalism  which  have  been  constantly  improved  over 
last  few  years,  Klepaczko  (1988a).  The  formal  ism  presents  a  consistent 
approach  to  the  kinetics  of  plastic  flow  of  metals  and  alloys  with  FCC, 
BCC  and  HCP  lattices.  The  rate-sensitive  strain  hardening,  instantaneous 
rate  sensitivity  and  temperature  are  taken  into  consideration  in  terms  of 
evolution  of  microstructural  state  variables.  The  notion  is  assumed  that 
the  current  mechanical  properties  of  a  metal  depend  entirely  on  its 
current  microstructure ,  and  any  changes  in  micros  true  Lure  result  in 
changes  of  mechanical  properties.  Consequently,  the  formalism  has  been 
applied  in  which  thermal  activation  analysis  is  employed  for  both  the 
kinetics  of  glide  and  kinetics  of  microstructural  evolution,  Klepaczko 
(1974).  Although  the  formalism  is  to  some  extend  similar  in  its  construc¬ 
tion  to  that  based  on  the  mechanical  threshold  stress  as  a  state  varia¬ 
ble,  Follansbee  and  Kocks  (1988),  it  is  more  general. 

The  formalism  employed  in  this  study  as  the  background  for  discussion  of 
microstructural  evolution  is  based  on  the  stress  partitioning  at  constant 
structure.  Thus,  the  flow  stress  in  shear  t  at  constant  structure  is 
given  to  a  good  approximation  by 

T(r,T)STR  =  (r  *t)str  +  t*^,T^STR 

where  r  tand  t*  are  respectively  the  internal  and  effective  stress  compo¬ 
nents,  "rand  T  are  plastic  shear  strain  rate  and  absolute  temperature. 
The  kinetics  of  defect  (dislocation)  movements  interrelates  at  constant 
structure,  characterized  by  j  state  variables  s, ,  the  instantaneous  va¬ 
lue  of  effective  stress  t*.  Whereas,  the  internal1  stress  r  must  be  also 
rate  and  temperature  dependent  via  dynamic  recovery  processes  (relaxation 
of  long  range  internal  stresses)  occuring  due  to  annihilation  and  rearran- 
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•gement  of  strong  obstacles  to  dislocation  motion.  Since  the  microstructure 
undergoes  an  evoUition,  and  the  state  of  microstructure  is  defined  by  s, 
state  variables,  the  state  variable  evolution  is  assumed  in  the  form  of  a 
set  of  j  differential  equations  of  the  first  order,  Klepaczkb  (1987, 
1988b)  "  ii' 

ds/  -  .  ,  -  ;  '  v  ■■‘X 

-i  *  Csk  ,  r(r)  ,  T(r)]  ,  k  •  1  ...  J  (2) 

dr  J 

Assuming  that  microstructure  is  characterized  by  five  instantaneous  quan¬ 
tities  :  p,t  p  ,  d(p),  D  and  a,  where  p,  and  p  are  the  mean  immobile  and 
mobile  dislocation  densities,  d(p)  is  the  subgra irt  diameter"* which  evolues 
with  p(,  0  is  the  grain  diameter  and  a  is  the  mean  distance  between  twins 
or  microtwlns,  the  internal  and  effective  components  of  stress  are  di¬ 
rectly  related  tp  those  quantities,.'  Evolutions  of  each  of  five  introduced 
state  variables  are  discussed  in  view  of  the  papers  submitted  to  this 
session'^  ■ j, 

in  addition,  an  evolution  of  microstructure  specific  for  shock  wave 
loading  arid  microinstabilities  in  the  form  of  the  adiabatic  micro-shear 
bands  are  also  taken  into  consideration.  Finally,  a  general  discussion 
and  more  specific  conclusions  are  offered. 
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OvarvUv  of  impact  properties  of  monolithic  ceramics 
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The  generic  term  'ceramics'  embraces  a  large  variety  of  material*  that 

r?ek*  >«  u-*»*  SrSiS 

*1°' i  Thla  r,vl,w  i*  only  concerned  vith  ceramics  of  high  tensile 
cefSida'.1^ ov, porosity  and  restively  free  from  defects,  often  called  'fine 

^b*  e!*<l^?ot?r1,*t^Cin  the  mechanical  properties  of  ceramics  even 
varllhlHtv6  12,dl‘**'  pMl,,J'ts  difficulties  associated  vith  their 

!  rky‘  Practura  toughness  is  often  measured  by  means  of  indentation 
tafts‘  ?hes*  h“ve  under  considerable  criticism.  An  alternative  is 

cwtiSver  te®ts  “sln«  sln«le  edK«  notched  beams  (SEN),  double 

cantilever  beams  (DCB)  or  double  torsion  specimens.  In  all  these  cases  the 

be^lmnorr«nr  thr  lnltlatlon  notch  crack  is  introduced  has  been  shown  to 
*  °  l!°'?  ar*  Prone  t0  suffer  a  severe  strength  degradation  as 

?!  t?*\.ma5hlnl,n®  P«««  and  when  the  notch  is  produced  by  plunge 
?  *  f*th?f*dfWh!ther  the  subse<>uant  specimen  behaviourPis  * 

testa  nfd!hetT^L  £  °f  the  undama*ed  bulk  material.  Low  speed  impact 
tests  of  the  Izod  and  Charpy  types  have  been  interpreted  in  terms  of  K 

assuming  that^stetic  equilibrium  is  reached  and  considering  inertial  loading 
by  applying  the  Duhamel  integtal  method  to  the  raw  force-  time  data  provided 
by  the  impact  tester.  This  approach,  while  preferable  to  a  simple  static 

“  KmnSS^,:' ’»'I,,tJr.Sp™L",,on  ,hlc''  ,s 

Fracture  then  depends  on  the  stress  wave  velocity  and  the  stress  Dulse 

duration  or  stress  rate.  The  mechanism  consists  in  the  formation  of 

clusters  of  microcracks  and  their  subsequent  growth.  Microcracking  has  also 

loadinff°Wwht0  expla^n  the  beh«viour  of  multi-phase  ceramics  under  shock 

different  iechanir«lSiare^f°rmed  at  t!?e  boundaries  between  phases  with 
different  mechanical  impedance.  Cracks  may  be  formed  normal  to  the 

direction  of  propagation  of  the  shock  or  along  that  direction,  A  damage 

mechanics  approach,  adapted  to  the  mathematical  modelling  of  failure  In 

,1,h  n”’r,c*1  inaiysis' *•  p°>— .uTSiK  si  ;r.«u 

Accu«%8n!«d^H„USed  f°?  plating,  backed  by  a  resilient  support. 

Accurate  predictive  numerical  analysis  are  not  yet  possible  in  the  ,h..n.a 

hi???!!!  * purely  «raPirical  approach  is  normally  followed  to  obtain  the 

ballistic  efficiency  of  the  ceramic.  This  depends  on  the  HEL,  the 
compressive  strength,  the  density  and  the  wave  velocity  but  the  overall 
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efficiency  of  the  shield  will  also  be  a  function  of  tha  mechanical 
impedances  of  tha  front  layer  and  tha  backing)  of  tha  properties  of  tha 
adhtplve  and  of  tha  impact  itself.  Tha  mode  of  failure  changes  vhan  tha 
ceramic  is  unsupported. 

Erosion  and  small  particle  impact  are  briefly  mentioned.  This  review 
bonsideri  only  the  response  of  a  ceramic  to  a  single  loading  at  room 
temperature.  The  variation  of  mechanical  properties  vith  temperature  is  a 
question  of  obvious  Interest  about  which  vary  little  has  yet  been  done  as 
is  the  effect  of  repeated  impact  on  the  residual  strength.  Secondary 
causes)  such  as  machining  damage  are  known  to  have  an  Important  effect  on 
the  scatter  of  properties  as  well  as  on  the  strength,  so  much  so  that  they 
can  even  mask  the  effect  of  the  primary  loading.  The  resistance  to  impact 
is  also  impaired  by  the  presence  of  an  applied  static  load  Inducing  tensile 
stresses.  The  conclusion  ia  therefore  that  much  still  remains  to  be  done  to 
understand  fully  the  behaviour  of  high  strength  ceramics  and  build  up  a  bank 
of  design  data  that  will  permit  to  exploit  their  potential  with  confidence. 
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Composite  materials  usually  comprise  a  matrix  material  such  as  polymer, 
metal  or  ceramic,  and  reinforcing  particles  or  fibres  such  as  carbon, 
glass,  polymer  or  ceramic.  By  combining  the  properties  of  the 
constituents,  with  the  correct  interface  conditions,  compositie  materials 
offer  the  designer  many  potential  advantages  over  homogeneous  materials, 
such  as  high  strength  and  stiffness,  low  density,  elastic  anisotropy 
giving  tailored  deformations,  preferred  fracture  mechanisms  giving 
increased  toughness  and  improved  high  temperature  stability. 

Many  of  the  applications,  where  the  designer  would  like  to  make  use  of 
these  properties,  Involve  some  form  of  impact  or  shock  loading,  for 
example  in  aerospace  and  transport  structures,  iri  armour  applications,  in 
reciprocating  machinery  and  in  sports  goods.  It  is  important  to  know  the 
response  of  the  materials  at  different  loading  rates  appropriate  to  the 
appl ication. 

Increasing  the  loading  rate  on  a  component  or  test  specimen,  to  rates 
representative  of  impact  loading,  can  cause  changes  in  response  and 
damage.  These  may  be  because  of  strain  rate  effects,  which  alter  the 
behaviour  of  the  material  on  a  microstruetural  scale,  producing  effects 
such  as  increased  failure  stress  or  reduced  fracture  energy.  Alter¬ 
natively,  the  component  or  specimen  may  respond  differently  on  a 
macrostructural  scale,  so  that  higher  frequency  modps  are  excited,  or 
stress  waves  may  not  be  able  to  travel  far  in  the  ti  °  available  for  the 
impact  event  and  energy  may  be  more  concentrated  in  the  region  of  the 
impact. 

There  are  many  different  test  methods  for  investigating  impact  phenomena 
in  composite  materials,  with  strain  rates  ranging  from  10-Vs 
(quasistatic)  to  3.0 Vs,  Dropwelght  or  pendulum  impact  teste  are  useful 
for  simulating  low  velocity  impact  (a  few  m/s)  on  carbon  fibre  composites 
(Dorey  1987,  Epstein  et  al  1989),  Gas  guns  may  be  used  to  achieve  impact 
velocities  of  several  100  m/s  (Dorey  1987,  Beaumont  et  al  1989),  and 
ballistics  can  achieve  several  1.000  in/s.  High  strain  rates  can  be 
achieved  by  other  means  such  as  by  Hopkinson  bar  techniques  (Lankford 
et  al  1.989,  Harding  et  al  3  989)  and  laser  induced  shock  (Gilath  et  al 
1.989). 

In  -studying  impact  events  it  is  necessary  to  monitor'  loads  and  damage 
occurring  in  very  short  times  and  this  may  be  done  using  strain  gauges 
(Harding  et  al  1989,  Beaumont  et  al  1.989),  moire  fringes  for  strain 


(Epste'n  et  al ) ,  high  speed  cameras  to  follow  damage  (Giiath  et  al  1989, 
Beaumont  el  at  1989)  or  crack  meters  (Lankford  et  al  1989). 

The  damage  produced  by  impact  loading  can  be  similar  to  that  produced  by 
quasiatatic  loading,  for  instance  where  elastic  deformation  is  dominant, 
or  the  shock  loading  may  induce  different  failure  modes.  In  compoBiteB, 
the  damage  pattern  can  be  dominated  by  weak  fracture  paths,  such  bb  at 
fibre/matrix  interfaces,  responding  to  secondary  stresses  in  the  complex 
load  field. 

Modelling  impact  events,  either  physically  (Beaumont  el  at  1989)  or  by 
computer  techniques,  is  important  for  understanding  and  predicting 
behaviour.  It  also  gives  valuable  clues  in  the  development  of  improved 
materials.  But  modelling  impact  performance  in  compsites  is  complex  and 
difficult.  Sometimes,  to  get  adequate  predic i tions ,  it  is  necessary  to 
simulate  the  event  experimentally  as  closely  as  possible,  as  in  the  case 
of  bird  strike  on  aircraft  structures.  But  full  scale  tests  of  this 
kind  can  bo  very  expensive,  and  laboratory  impact  tests  can  be  very 
useful  if  they  are  interpreted  intelligently. 
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Thla  article  focuses  on  the  modelling  of  material  properties  for 
struoturee  subjected  to  large  dynamic  loads  causing  extensive  material 
lnelaatlo  flow,  Many  articles  have  been  published  on  the  behaviour  of 
simple  structures  subjected  to  large  dynamic  loses  [1-5],  and  the 
literature  oontains  a  rich  coverage  of  the  well-known  rigid,  perfectly 
plastic  model  [6.7]. 

The  simple  rigld-plastlc  idealisation  gives  reliable  predictions  for  some 
impaot  problems,  but  it  is  not  always  adequate  for  the  increasingly  wide 
range  of  praotical  problems  whioh  confront  a  designer,  Nevertheless,  in 
Reference  [8],  several  recent  studies  are  examined  which  confirm  that  the 
rigid,  perfectly  plastic  model  is  an  acceptable  approximation  for  an 
elastic,  perfectly  plastic  material,  provided  the  dynamic  energy  input  is 
muoh  larger  than  the  maximum  wholly  elastic  strain  energy  oapaoity  of  the 
structure,  and  that  the  pulse  duration  is  short  relative  to  the  corres¬ 
ponding  natural  period. 

The  influence  of  material  strain  rate  sensitivity  is  important  for  some 
materials  and  structural  designs  and,  in  these  cases,  should  be  incor¬ 
porated  in  the  model  of  a  material.  This  phenomenon  is  fairly  well 
understood,  and  a  considerable  body  of  experimental  data  now  exists. 
However,  the  quest  for  more  efficient  and  lighter  structural  designs, 
coupled  with  greater  environmental  demands,  often  require  estimates  for 
the  failure  of  structural  systems  under  severe  dynamic  loads.  These 
calculations  require  Information  on  the  failure  modes  c?  structures  and 
the  strain  rate  sensitive  properties  of  materials  wh^ch  undergo  large 
plastio  strains  up  to  rupture. 

Section  2  of  the  paper  di3cus3es  some  recent  experimental  and  theoretical 
studies  on  the  dynamic,  inelastio  failure  cf  beams,  and  contrasts  the 
behaviour  under  uniform  impulsive  velocities  and  local  impact  loads.  It 
is  observed  that  the  failure  of  beams  subjected  to  local  impact  loads  and 
uniform  impulsive  velocities  is  not  well  understood.  Further  experimental 
data  Is  neoessary  to  identify  more  clearly  and  model  properly  the  failure 
modes,  particularly  for  mild  steel  beams  under  looal  impaot  loads. 

Some  further  observations  on  the  dynamic,  inelastic  failure  of  structures 
are  given  in  Section  3  of  the  paper,  while  Section  4  discusses  several 
outstanding  problems  on  the  Influence  of  material  strain  rate  sensitivity. 
It  is  noted  that  the  strain  rate  sensitive  characteristics  of  materials 
*r«  wtll  undtrstood  for  fimal!  jitralnn,  whftrRRn  ntrunhurfll  nranhwnrthlnftflfl 


studies.  and  other  practical  applications,  nay  involve  large  inelastic 
strains.  A  modification  of  the  Cowper-Symonds  constitutive  law  is 
proposed  in  order  to  model  the  observed  deorease  of  strain  rate  sensi¬ 
tivity  with  increase  in  inelastic  strain. 

Finally,  a  simple  model  for  the  variation  of  rupture  strain  with  strain 
rate  is  discussed  in  Section  5  and  compared  with  some  experimental  data 
C9-133.  However,  it  is  evident  from  the  available  experimental  data  that 
the  influence  of  strain  rate  on  the  rupture  elongation  of  metals  is 
unoletr.  Further  careful  experimental  tests  are  required  to  resolve  these 
dlfficultlec  before  the  proposed  equation  or  any  other  oomplex  expressions 
tiK]  are  employed  in  design  and  for  numerical  calculations. 

It  emerges  clearly  from  this  article  that  further  experimental  studies  are 
required  in  order  to  develop  reliable  models  for  the  material  properties 
and  structural  behaviour  under  large  dynamic  loads. 
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There  is  considerable  current,  interest  in  quantifying  the  initiation  end 
propagation  of  ductile  cracks  under  high-rate  dynamic  loading.  Several 
experimental  R-curve  techniques  including  the  Chioperfield  multi-specimen 
5 -curve,  the  dynamic  DC  potential  drop  method  and  for  steel  at  certain 
temperatures,  the  transition  from  ductile  .tearing  to  cleavage,  are  now 
available  to  allow  measurements  to  be  made'w'th  reasonable  precision  and 
repeatability,  These  relatively  simple  techniques  will  be  described,, their 
advantages  and  limitations  discussed  and  compared  with  other  methods 
requiring  more  elaborate  equipment. 

The  resuiis  from  a  number  of  oynamic  experimental  programmes 
conducted  at  imperial  College  on  structural  and  pressure  vessel  steels  will 
be  presented.  These  data  have  been  obtained  mainly  from  tests  with  sma;i 
3-point  bend  specimens  and  usually  as  a  by-product  of  other  research  aims, 
They  show  a  consistent  trend  of  increased  initiation  toughness  and 
increased  slops  of  tne  J  R-curve  for  higher  loading  rates,  and  a  decrease  in 
these  properties  when  lateral  constraint  is  raised  by  means  of  side-grooves 
The  limited  results  on  specimens  of  varying  absolute  side  give  conflicting 
results.  To  further  investigate  this  area,  a  series  of  tests  is  now  m 
progress  using  specimens  of  two  steels,  BS  4360  30D  and  HY  130,  for  which 
the  static  properties  have  already  been  thoroughly  explored 


A  Modified  "Moving  Singular  Element"  Method  for 
Fast  Crack  Propagation 

Wei  Jian,  Liu  yuanyong 

Northwestern  Polytachnical  University 
Xi'an,  Shaanxi,  China 

Synopsis 

Thera  are  a  lot  of  Numerical  Methods  to  Solve  Dynamic  Crack  Propagation 
problems.  The  accuracy  and  efficieney  of  those  methods  are  greatly  based 
the  simulating  model  of  crack  propagation  and  the  basic  variational  equation 
for  dynamic  solution.  In  the  past,  the  propagation  models  are  “Nodal  Fore# 
Release  Method"  tn  and  "Moving  Singular  Element  Method"  1,1  .  The  former 
is  much  rough  in  theory,  so  can't  simulate  the  crack  propagation  accuratiy 
but  it  is  easy  to  perform  and  takes  less  CPU  time  of  computation;  The 
latter  can  express  the  crack  running  more  accuratiy,  but  it  has  some' 
serious  drawbacksidifficult  to  creat  meshs  and  a  large  amount  of  running 
time.  As  to  the  variational  equation,  Nishioka  etc.  presented  an  "energy 
consistent  variational  statement"  regarding  for  the  crack  prppagation. 

But  the  results  show  that  this  variational  eqution  is  nor  quite  success, 
it  is  ambiguous  in  theory  and  results  in  an  oscillation  in  actual 
computation. 

In  this  paper,  we  will  present  a  modified  method  of  Ref [2]  in  which 
a  new  simulation  of  propagation  and  variational  equation  are  used  in 
order  to  combine  the  advantages  of  above  two  methods. 


Basic  Theory:  In  a  cracked-body  in  which  the  crack  moves  at 

velocity  v,  we  have  the  displacement  and  stress  eigen-function  near  the 
crack  tip  as  follows: 

“—ukf-  45,  s?-als?-y(?f  1  ’ ( ' J 

)S,  if  .n|-v,<w  r-i'> 'P-  'I 


where: 


Uxxn=, 

Uyyn". 

Uxyn=. 


(To  be  ignored) 


Sf  «l-(V/c4)*  ,  s‘,  -1-(V/C,)* 

|  4S,S,/(1+S'  )  /  n=odd 
7  l  l+s?  ,  n»even 

when  v»0,  the  equation  (1)  changes  into  the  well-know  williams-functions, 
which  is  deduced  in  the  static  condition. 

From  the  oringal  definition,  the  energy  relaae  rate: 

,s 

G-lim  4-1  Cxx(x)Uy  (x-s)dx  . (2) 

I—  0  * 

If  we  use  the  propagation  eigen-functions  and  williaims-functions  as  the 
basic  displacement  and  stress  function  of  singular  elements— 'called 
propagating  singular  element  and  static  singular  element  respectively. 

We  can  derive  the  different  energy  release  rate  Gs  and  Gd  using  (2): 

•“»  ^-txYiv)  iVcc.  3.  > 


where: 


(plane  strain) 
(plans  stress) 


(3) 


u 


j T  In  the  meantime,  from  another  expression  of  G: 

If  r  _  1  ,  dw  dv  dt  . 

;  Ct-Ff3a"ai~3r) 

We  can  see  that  if  the  propagation  singular  element  is  replaced  by  the 
static  singular  element,  there  is  only  a  little  change  in  total  energy 
[  release.  Because  the  singular  element  only  takes  a  little  part  in  whole 
i  body,  we  can  give  an  approximation: 

4  Gd-*GS  . (4) 

Substitute  (3)  for  (4) ,  we  get: 

1]  . .,5, 

Now,  using  (5),  we  can  largly  reduce  the  computing  time.  Because  we 
need  no  longer  to  reproduce  singular  elements  stiffness  and  mass  matrices 
when  v  is  varying,  which  take  a  lot  of  time  to  perform,  and  the  programs 
is  much  simplified  as  well. 

As  to  the  variational  aquation,  Because  inany  genetic  time  ti  ,  wa 
can  consider  the  crack  is  not  propagating,  So  unlike  [2],  we  can  use  a 
simple  variational  equation  to  a  running  crack: 

/(tfij*xy+fuifui)dv-  /  Ti<Tuids-  /  TiJuids»0  . (6) 


Finally,  we  got: 

K-q+M*q-Q*0 

This  eq.  was  solved  by  New  mark-/B  time  integration. 


(7) 


EXAMPLE:  We  first  use  a  steel  DCB  specimen:  Kla  =7066N/mm',r 

v-lOOOm/s, At=3,«s,  All  the  experiment  data  were  got  by  J.F.Kalthoff  . 
Fig.l  shows  the  comparasion  rasults  of  different  variational  eq.  using 
propagation  singular  element.  Fig. 2  shows  the  comparasion  of  propagation 
singular  element  and  static  element  under  new  variational  eq. ;  the  latter 
takes  25%  less  CPU  time  than  former.  Second,  we  use  Hemilate-100  DCB 
specimen:  Kia  =>73 . 36N/mm"'  ,  v»295m/s,  t»6^s,  we  got  nearly  the  same 

result  compared  with  the  [2],  while  the  time  increment  at  is  30  times 
large  then  [2]. 

CONCLUSION: 

This  method  is  more  accurate  than  "Nodal  Force  Release  Method",  and 
more  efficient  then  "Moving  Singular  Element  Method".  And  compared  with 
the  Experimental  results,  it  shows  a  good  consistence.  The  stability  and 
convergency  of  this  method  is  also  attractive,  especilly,it  can  be 
performed  in  such  large  increment  of  at. 


References:  (To  be  ignored) 
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HEAT  GENERATION  DURING  FRACTURE 


By  Akira  KO0AYASH I , Pr of es s or , Dr . 

NobtJO  OHTANI,  Research  Fellow,Dr. 

Michel  GOROC,  Research  Student 

Faculty  of  Engineering, Department  of  Materials  Science, 
The  University  of  Tokyo, H  on  go  7-3-1,  Bunkyo-ku, Tokyo  113 
Dapan 


In  the  present,  report, the  validity  on  the  estimation  of  dyna 
mic  energy  release  rate  by  generated  heat  is  presented. 

The  specimen  in  interest  is  PMMA .  Experiments  are  done  under 
various  pre-stratn  conditions  in  crack  propagati on , during  which 
generated  heat  in  the  function  of  temperature  rise  is  measured  by 
thermistors  and  the  crack  velocity  by  velocity  gages.  A  crack 
is  initiated  by  an  impact  of  a  sharp  edge  on  a  notch  in  the 
specimen  to  realize  the:  crack  propagation. 


Dynamic  energy  release  rate  Gjj  for  a  finite  width  elastic 
body  with  fixed  displacement  is  expressed  as 

GId  -  1/2  f(V  )  •  E  -  h  •  £p?  ---(1) 

where  f (2/)  =  1/(1  - 1 ')  for  plane  stress 

“  1  -  V  /(I  +p)(l  -  2  L>)  for  plane  strain 
and  E^Young's  modulus , U  =Poisson ' s  ratio, h=a  distance  between 
jaws  for  gripping, and  £ p=pte-strain.  (Ref.l) 

As  for  heat  generation , assuming  the  one-dimensional  heat 
conduction , generated  heat  can  be  expressed  as 

Q  -  JTtT  fef-cy.A  Tm  ---(2) 

where  -£-*2.71828---,  f  “density , c«specific  heat, and 
T  ■  temperature  rise.  (Refs. 2  &  33 

Measurement  of  pre-strain  and  temperature  rise  in  combination 
with  Eqs.(l)  and  (2J  leads  to  Fig. 3 
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Fig.  3  Dynamic  Energy  Release  Rate/Generated  Heat  vs .  Pre-Strain 
Fig. 2  Crack  Propagation  Velocity  vs.  Pre-Strain 


Further  combination  of  Figs. 2  and  3  presents  a  solid  line 
in  Fig. 4 ,in  which  previous  experimental  results  through  a  caustic 
method  (Ref. 4)  are  shown  by  solid  circles. 


Fig. 4  Dynamic  Energy  Release  Rate/Generated  Heat:  vs.  Crack 
Veloci ty 

As  seen  from  Fig . 4 , dynamic  energy  release  rate  and 
generated  heat  are  in  good  agreement  until  the  crack  velocity 
reaches  about  550  m/s .which  is  almost  seventy  percent  of  terminal 
velocity,  Therefore , the  validity  on  estimation  of  dynamic  energy 
release  rate  by  generated  heat  holds  in  the  range  mentioned  above. 
One  reason  of  this  deviation  from  550  m/s  is  explained  by  fracto- 
graphy  of  fractured  surfaces , showing  much  microscopic  crack 
branching  leading  toward  much  energy  consumption. 


THREE-DIMENSIONAL  SEMI-INFINITE  CRACKS 
SUBJECTED  TO  DYNAMIC  LOADS 


P.  Pintado1  and  F.G.  Benitbz1 

Escuela  Superior  de  Iagenleros  Industries,  Avda.  Heina  Mercedes,  41012  Sevilla, 
Spain 

The  boundary-integral  equation  method  for  elastostatics  and  elastodynamics  is 
especially  heipfull  when  the  particular  geometry  of  the  problem  being  analyzed 
Is,  to  some  extent,  acounted  for  by  the  fundamental  Green’s  function  used.  In 
this  paper,  the  three-dimensional  solution  for  an  impact  point  load  acting  in  the 
interior  of  an  unbounded,  homogenous,  linearly-eiastic,  isotropic  thick  layer  has 
been  studied.  This  fundamental  solution  is  used  in  connection  with  the  bouadary- 
integral  equation  method  to  obtain  the  strew  and  displacement  fields  in  problems 
concerning  semi-infinite  through-cracks  in  plates  subjected  to  dynamic  loads.  The 
technique  proposed  can  be  used  for  the  analysis  of  other  geometrical  configurations 
such  as  presurised  cracks  and  circular  holes  in  plates  of  arbitrary  thickness. 
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TENSILE  INSTABILITIES  AT  HIGH  STRAIN  RATE'-.  IIO.T’n  AND  MULTIAX|AL 
EFFECTS. 


if.Fressengeas .  A.Notlr'jri 

Laboratolre  da  Physique- «>.  He can l qua  d«$  MaUrlaux.  U.A,  c.N.R.s.  1215 
L'nlversll*  de  Nets.  tie  du  Saulcy,  570i'5-  Hals  Cedex.  France. 

.  High  speed  metallic  jets  are  formed  by  tHe-ccl  lapse'  of  a  .conical  or 
wedge-shaped  l tner  under  explosive  leading.  Mypi  :al  ve Uclt ie*  are  fckti's 
at  the  tip  and  4Vpi/ 3  at  the  tall  of  the  .jet!  consequently,  the  jots  are 
submitted  to  an  enormous  stretching!  their  length  can  be  multiplied  4 
times  In  100ms.  However,  beyond  some  standoff  distance,  th«l  jets  neck 
down  In  a  series  of  locations  and  soon  break  dowr  into  closely  similar 
fragments.  -This  Is  known  to  .limit  the  Jet  penetration  ability  CBlRKHOFf 
ek  al.  1948),  Such  a  phenomenon  resembles  the  breakup  of  efcpil l arg¬ 
il  quid  jets:  therefore  it  has  been  suggested  that  surface  tension  Is  the 
destabilizing  mechanism,  and  classical  capillary  instability  theories 
have  been  applied  to  the  metallic  jets  breakup  CFRANKEL  and  WEIHS.1385). 
Vet.  the  perturbation  growth  Issued  by  such  models  Is  too  small  by 
several  orders  of  magnitude  when  It  comes  to  compare  the  predicted 
breakup  time  with  the  observed  data. 

In  this  paper,  the  metallic  jet  breakup  is  rather  "lowed  as  a  necking 
instability:  It  is  shown  that  the  necking  instability  analysis  of  a  non 
linear  vlscoplasilc  material  can  provide  perturbation  growths  large 
enough  to  account  for  the  observed  breakup  time  and  fragments  aspect. 
The  driving  mechanism  Is  the  geometric  softening  0:  the  Jet  due  to  Its 
section  reduction:  it  renders  the  uniform  stretching  jet  unatde  to 
sustain  the  axial  force  lnuol"ed.  In  the  meantime,  stabilizing  factors 
are  at  work,  which  delay  the  growth  of  non  uniform  perturbations:  in 
addition  to  the  damping  of  rite  sensitivity,  the  axial  and  lateral 
Inert  la. effects  slow  down  the  growth  of  long  wavelength  perturbs t ions. 
As  a  matter  of  fact,  the  latter  involve  larger  mass  amounts.  The 
mult l axial  aspects  of  *he  stress  field  In  the  neck  do  so  at  the  short 
wavelengths:  the  average  axial  stress  applied  to  the  neck  cro33  -section 
Is  larger  at  long  wavelengths:  consequently,  for  a  giver,  applied  to’-ce. 
the  section  area  Is  smaller,  which  means  that  the  instability  proceeds 
faster  at  long  wavelengths,  Be'ween  the  two  foregoing  limiting  cases',  an 
optimum  non  zero  finite  wavelength  Is  selected  at  each  moment,  for  which 
the  perturbation  growth  is  maximum.  This  is  Indicative  of  multiple 
necking  and  of  jet  fragmentation. 

Due  to  the  stretching,  the  basic  uniform  plash:  flow  is  time 
dependent.  Therefore,  the  problem  cannot  be  seined  by  using  ‘he 
classical  simple  eigenvalue  formulation:  the  perturbation  growth  Is  not 
that  of  a  simple  exponential,  and  there  is  no  single  dominant 
wavelength.  In  this  paper,  we  consider  a  "quasl-cleady"  uniform  plastic 
flow:  It  Is  assumed  that  the  time  scale  of  perturbation  growth  Is  much 
smaller  than  the  uniform  stretching  jet  time  scale.  This  Is  welt 
verified  In  all  the  specific  cases  we  Investigated.  Then  !t  Is  shown 
that  the  problem  can  be  addressed  through  a  series  of  Instantane  ..is 
eigenvalue  problems.  At  each  moment  a  different  dominant  wavelength  and 
a  different  pertubatlon  growth  rate  are  found.  The  procedure  provides 


predict  ions  for  the  breakup  *  1’fl*  and  >.he  frajinw v .■  aspect  r.r  to  t  length 
to  maximum  diameter) i  at  fcredkup  .  t ime  (!  lie.  when  .  ih  •■ .  non.  uniform 
perturbed  radlut  Is  zero.',  the  ln*t..intanaouf  >Tit,i:al  wavelength-  It 
Indicative  of "the  the  fragments  size. 

Tha  comparison  with  the  aval  labia  experimental  >lala  on  copper  Jets 
(CHOU  a l  al .  1977?  suggests  l hat  both  predict \ on s  are  realistic,  and 

that  the  model's  allows  a  consistent  .prediction  of  the  jet  breakup 
phenomenon.  The'  cofnpat ibi 1 1 >y  with1  the  observed  data  Is  achieved  for 
material  constants  (reference  stress,  rata  sensitivity)  l”  agreement 
with  the  accepted  experimental  trend.  Such  controlling  parameter*  of  the 
Jet  breakup  as  lateral  Inertia  and  rate  sensitivity  are  pointed  out.  and 
their  Influence  dlsjcus.sed,  -It,, is  shown  that  a  large  ra*e  sensitivity 
yields  longer  fragments ,,Sxlal  and  lateral  Inertia  effects  are  more 
complex:'  In  highly  dynamic  cases,  the  fracture  process  is  delayed  a' 
larger  plastic  deformations.  Th®  fragments  are  shorter 
dynamic  events  however,  Hlthou.gh  It  focused  on  metallic 
analysis,  plight  .have  Implication*  on  similar  processes, 
fit-re  spinning.  -  :  i  . 


than  In  less 
Jets,  such  an 
like  polymer 


References  i 


V  ;■  '  :  :•  v-.;  -  ’  ■*,  ■_  ';<■  -  .  ft  /•<;■■■ ; 

BIRKMOFP,  (3 ...  MACPOIJ0RILU . .  D.P.„  PUGH ,  E.H.  and  TAYLOR  •  '?-M  l^dS. 

Explosives  with  lined  cav,l>  ies ...  J.Gppl  .Rhys.  I?. 

CGPLEGNE.  CHOU.  r-.C.  and  .CICCAREU1 .  F  0.  i  19?''.  shaped  charge  Jet 
stability  and  pene.trat  1  on!  calcutat l on* PFL-CP.-f51  Report  . 

FRAHkel.  !  .  and,  WE  I  Ml ,  D.s.  1999.  Stability  cf,  a  jeap.i  l  U'ry  Jet  with 
linearly  increasing  axial,  veld  cl  t’y  (with  application  to  shaped  charges), 
J.Fiul  l.Mech..  IW.  £69-90?.  ,  i  ;.  :j  '  . ;  ; 


f3sj 
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A  cooperative  research  program  has  been  developed  on  the  dynamic  tensile  testing  facility  built  by 
CEA.  Results  obtained  on  smooth  bars  up  to  500  s'1  suggest  that  a  linear  strain  rate  dependency  of 
stress  can  be  used  to  represent  the  constitutive  equation  of  various  high  purity  metals  or  alloys  (such 
as  Copper,  Tantalum  or  Copper-Tungsten),  as : 

ff*Cr0  +  PF  (1)  with  Ce  =  KE"  (2). 

There  remained  a  problem  with  dynamic  tensile  testing  due  to  transient  stages  of  heterogeneous 
strains  caused  by  Inertial  efiects.  Recent  numerical  simulations  using  a  one  dimensional  or  a  two 
dimensional  Units  element  method  show  that  Inertia  I  effects  at  2000  s-'  are  significant  only  at  the  very 
beginning  of  the  test,  but  do  not  alter  the  validity  of  the  linear  stress-strain  rate  relationship  assumed 
earlier. 

In  order  to  complete  our  investigation  about  mechanical  behaviour  of  high  purity  copper  at  high  strain 
rate,  we  performed  dynamic  and  quasistatic  tensilu  tests  on  notched  bars,  with  three  different  notch 
radii.  Our  work  deals  with  the  variation  of  the  maximum  load  value  and  of  the  intrinsic  ductility  as  a 
function  of  strain  rate  or  of  the  notch  radius.  Th<>  Intrinsic  ductility  Is  characterized  by  the  length 
relative  variation  of  the  notch  part  between  the  Initial  state  and  the  fractured  one.  These  experimental 
results  are  discussed  on  the  basis  of  the  two  dimensional  finite  element  simulation  which  takes  into 
account  inertia  effects. 

Our  computation  shows  that  whatever  the  strain  rale,  for  a  given  geometry,  a  same  relative  variation  of 
the  minimum  diameter  leads  to  similar  values  of  'he  equivalent  strain  or  of  the  stress  trlaxiality.  So, 
these  parameters  which  are  in  direct  relation  to  damage  effects  are  not  influenced  by  the  test 
conditions  and  consequently  by  the  inertia  effect! , 

in  all  the  cases  (figure  1),  the  maximum  load  valus  Increases  when  the  notch  radius  decreases.  This 
phenomenon  is  well-known  for  quasistatic  test  conditions  and  named  "notch  strengthening  effect" : 
it  Increases  with  strain  rate  in  dynamic  conditions. 
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figure  1 :  Maximum  load  values  as  a  function  of  notch  radius  at  different  strain  rates 

The  numerical  simulation  leads  also  to  an  increase  of  the  maximum  load  value  with  strain  rate,  which 
appears  at  the  same  macroscopic  elongation  or  at  the  same  relative  variation  of  the  minimum  diameter 
whatever  the  dynamic  strain  rate.  So,  the  equivalent  strain  being  the  same  In  the  minimum  cross 
section,  we  are  in  position  to  assume  that  the  maximum  load  value  Increase  Is  due  to  the  linear  strain 
rate  behaviour  of  the  high  purity  copper  (equation  1).  However,  our  computation  does  not  show  an 
increase  of  the  notch  strengthening  effect  with  strain  rate. 

Whatever  the  strain  rate,  the  intrinsic  ductility  Increases  when  the  notch  radius  decreases  (figure  2). 
This  result  can  be  justified  If  we  study  the  variation  of  the  equivalent  strain  along  the  minimum  cross 
section  at  the  same  relative  variation  of  the  minimum  diameter.  The  lower  the  notch  radius  value,  the 
higher  the  equivalent  strain  gradient,  but  the  lower  the  mean  equivalent  strain. 

notched  part  rupture  elongation  (%) 
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figure  2:The  notch  part  rupture  elongation  as  a  function  of  notch  radius 
between  quasistatic  and  dynamic  strain  rates 

On  the  other  hand,  we  are  not  able  to  explain  simply  with  these  computation  results,  the  differences 
observed  between  quasistatic  and  dynamic  strain  rates  :  the  mean  equivalent  strain  decreases  from 
the  dynamic  to  the  quasistatic  strain  rates. 

Nevertheless,  we  may  remark  that  the  dynamic  constitutive  equation  leads  to  higher  strain  rate 
sensitivity  parameter  m  than  In  quaslstatlc  conditions  where  rr,  is  particularly  low  ( m  -  0.012 ). 
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Results  from  high  strain  rate  tension  and  torsion  tests  using 
Split  Hopklnson  Pressure  Bars  have  been  combined  with  data  from  the 
literature  to  determine  the  effects  of  temperature,  strain  rate  and 
stress  on  the  deformation  and  fracture  of  polycarbonate  (PC).  See 
Fig.  1. 

An  Increase  in  strain  rate  results  in  a  modest  Increase  in  yield 
stress.  Decreasing  temperatures  result  in  Increased  yield  and 

fracture  stresses.  For  temperatures  below  -60°C  bond  rotation  is 

o 

frozen,  resulting  in  small  fracture  strains.  Above  -60  C  deformation 
is  by  sliding  of  backbone  molecules  and  longer  strains  are  achieved. 

For  temperatures  in  the  range  -60  C  to  140°C  the  response  consists  of 
three  stages;  I  -  non-linear  viscoelasticity  until  yield,  II  -  strain 
softening,  III  -  strain  hardening. 

Shear  bands  develop  In  torsion  and  compression,  and  necking  occurs 

in  tension.  Above  140°C  the  van  der  Waals  bonds  melt  and  the  material 
behaves  as  a  viscous  liquid. 

The  Eyrlng  Equation  is  shown  to  fit  the  available  data  for  yield 
of  PC  for  temperatures  T  in  the  range  0.6  <  T/Tg  <  1.0,  (Tg  *  145°C  = 
glass  transition  temperature)  and  for  strain  rates  between  10  4  s  1 

and  5  x  10  s  A  simulation  of  two  Eyring  terms  is  required  for 
temperatures  below  0.6  Tg.  The  model  accounts  for  the  different 
stress  states  associated  with  tension,  shear  and  compression. 

Strain  to  fracture  for  PC  drops  linearly  with  increasing  In 
(strain  rate),  with  compression  failure  occuring  at  about  7 OX  greater 
strain  than  tensile  failure.  Fracture  is  not  by  adiabatic  softening, 
as  orientation  hardening  stabilises  the  material. 

In  shear,  fracture  is  by  the  formation  of  tensile  microcracks 
which  then  coalesce  by  plastic  collapse  in  shear.  Compression  failure 
is  thought  to  be  by  a  similar  process.  Tensile  failure  is  by  the 
propagation  of  one  or  two  tensile  crackB  from  pre-existing  surface 
flaws,  following  sub  critical  flaw  growth. 
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In  modelling  dynamical  strain  and  fracture  of  metals,  the  mo¬ 
del  of  a  Maxwell-like  type  visco-elastic  body  is  used,  in  which  the 
relation  between  the  relaxation  time  of  shear  stresses  and  the 
medium  state  parameters  is  determined  on  the  basis  of  dislocation 
representations  of  plastic  deformation.  When  obtaining  the  V  rela¬ 
tion,  the  experimental  data  are  used  on  determining  the  dynamical 
yield  limit  as  a  function  of  strain  rate.  The  model  is  supplemented 
by  the  temporal  fracture  criterion  relating  the  life  time  with  the 
load.  The  criterion  generalizing  the  well-known  approach  suggested 
by  S.N.Zhurkov, describes,  in  termB  of  dislocations,  both  "static" 
sind  "dynsunic"  branches  of  the  experimentally  observed  dependence. 

Testing  of  the  proposed  model  and  of  the  numerical  methods  is 
realized  on  the  basis  of  some  one-dimensional  problems.  The  results 
are  compared  with  the  experimental  data.  The  capabilities  of  the 
model  and  of  the  method  are  illustrated  by  the  solutions  of  two- 
-dimensional  problems  of  high  rate  strain  and  fracture.  A  good 
qualitative  and  quantitative  agreement  between  the  predicted  and 
experimental  results  has  been  established. 
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Abetraot  -  The  applicability  of  dynamic  fracture  mechanics  to  the  failure  of  heavy 
metal  alloys  subject  to  high  velocity  impact  is  a  relatively  unexplored  area.  This  is 
principally  due  to  the  lack  of  knowledge  of  the  effect  of  loading  rate  on  deformation 
and  fracture  toughness.  Recent  progress,  however  has  been  made  towards  identifying 
the  failure  modes  involved  in  impact  loading  rate  In  tungsten  alloys  [1].  As  the  loading 
rate  increases,  the  fracture  process  changes  from  a  ductile  to  a  special,  brittle  type  of 
failure.  Specifically,  the  transition  from  ductile  to  brittle  is  associated  with  a 
decrease  of  fracture  of  the  tungsten  particles  (W)  and  the  nickel  base  matrix  (M)  and 
with  an  increase  of  cracking  along  tungsten/matrix  (W/M)  and  tungsten/tungsten 
(W/W)  Interfaces  [1,2].  Unfortunately,  due  to  the  high  strain  rate  sensitivity  of  this 
material,  identification  of  the  micromechanical  processes  alone  is  insufficient  to 
predict  fracture  toughness  at  ultra  high  leading  rates. 

The  purpose  of  this  investigation  is  to  report  fracture  toughnesses  measured  over  a 
wide  range  in  loading  rate)  interpretation  of  these  results  is  based  on  SEM 
observations  of  the  fracture  specimens  and  quantitative  micromechanical  modeling.  A 
novel  technique  has  been  used  to  measure  fracture  toughness  at  dynamic  loading 
rates  lij  varying  from  10®  to  3.5  x  10®  MPa/m  s' using  small  compact  specimens 
[3].  The  quasi-static  toughness  was  evaluated  using  a  standard  compact  specimen  at  a 
loading  rate  of  1  MPa/m  s"‘.  Mloromechanieai  modeling  was  performed  using 
simple  fracture  models  for  ductile  and  brittle  fracture,  following  a  procedure 
developed  in  a  previous  study  of  the  effect  of  loading  rate  on  fracture  properties  of  a 
plain  carbon  steel  [4].  Specifically  the  modeling  is  based  on  quantification  of  the 
mieromeohanlcal  processes  involved  in  ductile  and  brittle  fracture,  and  measurement 
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of  tensile  deformation  properties  at  strain  rates  of  8  x  10  a  and  1.5  x  10  s  . 

Decrease  in  toughness  was  observed  with  increasing  loading  rate  and  was  correlated 
with  increasing  apparent  embrittlement  of  the  macroscopic  fracture  surfaces.  The 
transition  from  ductile  to  brittle  failure  is  associated  with  an  increased  frequency  in 
the  tungsten/tungsten  separation.  Good  agreement  was  obtained  between 
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Abstract 

Over  the  past  nearly  two  decades  numerous  investigations  have 
been  carried  out  on  the  formation  and  propagation  of  adiabatic 
shear  bands.  At  high  strain  rates  adiabatic  shear  localization 
is  a  very  important  phenomenon  that  has  a  profound  effect  on 
plastic  deformation  and  fracture.  This  paper  describes  obser¬ 
vations  made  of  the  nature  of  adiabatic  shear  hands  in  commer¬ 
cially  pure  Ti  shock  loaded  at  a  strain  rate  of  lO^/sec.  The 
experiments  were  carried  out  at  controlled  total  elongations 
of  5  to  11  4,  which  was  achieved  by  varying  the  pedestal 
height  in  our  shock  design.  In  all  samples  tested,  shear  bands 
were  observed.  The  quantity  and  shear  band  width  increased 
with  increasing  elongation.  Above  10  %  total  elongation 
fracture  was  dominant  and  only  local  strain  up  to  fracture  was 
obtainable.  The  strains  were  measured  from  circle  grids  photo 
etched  on  the  outer  surface  of  the  samples.  This  allowed  for 
direct  external  measurement  of  shear  band  displacement  relative 
to  the  internal  microscopic  displacements.  The  local  strain 
between  shear  band  segments,  ie,  the  plateau  steps  on  each  of 
the  samples,  irrespective  of  the  overall  elongation,  had  little 
if  any  measureable  strain.  All  local  strain  readings  were 
within  the  lower  limit  of  the  grid  technique  and  did  not  exceed 
1  %.  Strains  in  the  shear  bands  that  were  measurable  did  exceed 
30  %. 

No  phase  transformations  in  the  shear  band  were  observed,  only 
grain  refinement  and  deformation.  Shock  hardening  was  observed 
over  the  unshocked  sample.  However,  all  the  post  shocked 
samples  had  essentially  the  same  valve  for  all  the  elongations. 
The  microhardneos  values  "/ere  equivalent  for  both  matrix  and 
the  shear  bands. 

For  all  conditions  of  low  strain  irrespective  of  total 
elongation,  the  shear  bands  were  typically  characteristic  of 
shear  bands  observed  at  much  lower  strain  rates.  However,  above 
approximatly  5  %  strain  the  shear  bands  had  a  very  definite 
boundary  between  the  shear  band  and  the  matrix,  and  increased 
in  width  with  increasing  total  elongation.  This  increasing 
boundary  effect  can  be  directly  attributed  to  the  strain  heat 
which  is  not  homogeneous  and  is  highly  localized;  thus,  giving 
rise  to  this  boundary  effect,  it  is  similar  in  nature  to  the 
heat  effected  zone  in  welds. 
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Adiabatic  Shear  Bands  in  One  Dimension 
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ABSTRACT:  This  paper  summarizes  the  progress  made  at  the  Ballistic 
Research  Laboratory  to  date  in  analyzing  and  understanding  the  dynam¬ 
ical  processes  thfu  cause  the  formation  of  adiabatic  shear  bands  and  the 
influence  of  various  physical  parameters  on  their  formation. 


A  simple  one  dimensional  model,  which  may  be  thought  of  as  simulating  a  torsional 
Kolsky  bar  test  on  a  thin  walled  tube,  has  been  studied  extensively  by  both  numerical 
and  analytical  techniques.  For  a  body  in  simple  shear 

x=A >«(}',/),  y  =  Y,  z  =  (1) 

the  nondimensional  equations  of  thermoviscoplasticity  may  be  written  as 


Momentum:  vt 

F.nergy:  0, 

Elasticity:  s, 

Flow  Law:  a 


Work  hardening:  Kt 


-  SJP, 

-  +  .n,, , 
=  /<K  -  iv)  < 

-  M(k,0)x %  . 


(2) 


where  v  is  the  velocity,  0  is  tiie  temperature,  a  is  the  shear  stress,  k  is  the  work 
hardening  parameter,  7P  is  the  plastic  strain  rate,  and  sgn  indicates  the  algebraic  sign 
of  its  argument. 


Only  insulated  boundaries  with  constant  prescribed  velocity  arc  considered.  That  is, 
0V(±1,<)  =  0  and  v(±l,  t)  =  ±1.  The  initial  strain  rate  is  assumed  to  be  close  to  unity 
everywhere,  and  the  initial  stress  is  assumed  to  be  constant.  To  trigger  a  shear  band, 
only  perturbations  in  the  initial  temperature  will  be  considered.  0{y,  0)  =  0o(y),  where 
0o  is  small.  Mechanical  perturbations  can  also  initiate  a  shear  band,  but  for  brevity 
we  omit  them  here, 


From  a  variety  of  analytical  technique*  it  has  been  possible  to  understand  many  of  the 
primary  parametric  dependences  of  shear  band  formation,  at  least  in  the  quasi-static 
approximation.  Some  of  the  principal  conclusions  are  that  shear  band  formation  is  a 

fruj 


multistep  process,  that  instability  and  intense  localization  are  not  coincident  in  time, 
that  a  fully  formed  shear  band  often  behaves  like  a  boundary  layer  with  a  well  defined, 
calculable,  spatial  distribution,  and  that  for  perfectly  plastic  materials  a  theoretical 
“shear  band  susceptibility"  can  be  identified  and  calculated  from  macroscopic  labora¬ 
tory  measurements.  The  susceptibility,  which  appears  naturally  in  simplified  analyses, 
appears  to  he  the  key  quantity  in  assessing  stability,  early  growth  rate,  and  minimum 
possible  time  to  localization.  Thermal  conductivity  has  a  substantial  delaying  effect  on 
intense  localization  at  lower  applied  strain  rates.  Elasticity  has  little  effect  until  intense 
localization  occurs,  and  then  it  may  add  rather  elaborate  structure  to  the  late  stage 
morphology.  Work  hardening,  which  is  known  to  be  stabilizing  before  localization, 
has  a  partially  stabilizing  effect  after  localization  where  it  also  adds  structure  to  the 
late  stage  morphology.  Finite  element  analysis  indicates  that  inertia  may  also  have  a 
substantial  delaying  effect  on  severe  localization,  but  the  exact  parametric  dependence 
is  not  yet  well  understood. 
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Abstract 


An  investigation  is  presented  of  the  initiation  and  growth  of  adiabatic 
shear  bands  in  a  titanium  alloy  (T1-6A1-4V)  due  to  the  impact  of  hard 
wedge-ended  and  cone-ended  projectiles.  For  projectiles  of  these  shapes, 
geometrical  similarity  of  the  impact  deformation  is  maintained  over  the  range 
of  velocities  explored  (30-300  m  s  i).  The  mean  strain  e  associated  with  the 
impact  indentations  therefore  does  not  vary  [1J  with  velocity,  in  contrast  with 
the  case  of  sphere  impact  examined  in  previous  work  [2-4], 

Hardened  tool  steel  projectiles,  with  included  wedge  or  cone  angles  of  60, 
90,  120  and  150°  and  with  constant  mass  (0.88  g)  were  fired  at  normal  incidence 
at  thick  targets  of  the  titanium  alloy  at  velocities  up  to  300  ms  .  The 
impact  craters  were  sectioned  and  examined  metallographically  to  detect  the 
presence  of  adiabatic  shear  bends.  The  velocities  at  which  shear  bands  were 
observed  varied  with  wedge-angle  in  a  manner  shown  in  Fig.  1;  the  relationship 
tor  impact  by  cone-ended  projectiles  was  almost  the  same,  except  that 
deformation  failed  to  localize  at  all  for  impact  by  150°  cones.  The  dynamic 
hardness  of  the  target  metal  was  calculated  from  plots  of  the  crater  volume 
against  kinetic  energy  and  was  used  in  conjunction  ^:.th^a  rigid-plastic  model 
to  determine  the  mean  strain  rate  (1  x  10*  -  4  x  10  s  )  as  a  function  of  the 
projectile  geometry  and  tho  impact  velocity. 

The  shear  band  patterns  were  of  two  main  types,  consistent  with  the 
"cutting"  (Fig.  2a)  and  "radial  compression"  (Fig.  2)3)  modes  of  plastic 
deformation  observed  after  quasi-static  indentation  by  cones  (5,6]  and 
wedges  [6].  The  zone  containing  the  shear  band  deformation  was  restricted 
closer  to  the  crater  surface  after  impact  by  sharp  cones  compared  to  sharp 
wedges  [6],  particularly  at  low  velocities  (Fig.  3). 

Mean  strains,  e,  associated  with  quasi-static  indentation  of  metalG  by 
cones  have  been  reported  by  Atkins  and  Tabor  [5],  who  measured  e  *  0.30,  0.25, 
0.17  and  0.08  for  indentation  by  60°,  90°,  120°  and  150“  cones  respectively. 

In  the  present  work  the  shear  bands  were  found  to  initiate  at  0.<J8  <  e  <  0.^7 
in  the  cone  impact  experiments,  at  mean  strain  rates  from  1  x  10  to  4  x  10  s 
This  result  is  consistent  with  earlier  work  [2j;  shear  bands  were  found  to 
initiate  beneath  craters  formed  by  spheres  at  e  =  0.12  over  a  similar  range  of 
strain  rate.  The  technique  further  demonstrates  that  a  critijal_jtrain  is 
necessary  to  Initiate  shear  bands  at  strain  rates  of  about  10*  s  where 
plastic  deformation  is  effectively  adiabatic. 


Present  address:  Alcan  International  Ltd.,  Banbury  Laboratories,  Banbury, 
Oxon.  0X16  7SP. 
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Fig.  1  Results  for  inpact  of  wedge-anded  projectiles  or.  to  T1-6A1-4V. 


Fig.  2  Shear  bands  propagating  along  outwardly-directed  (a)  and 

predominantly  inwardly-directed  (b)  trajectories  formed  by  impact  of 
60°  and  120°  cones  respectively  at  284  ±  2  m  s  . 


Fig.  3  Comparison  of  shear  band  deformation  beneath  craters  formed  by  impact 
of  60°  wedge-ended  (a)  and  cone-ended  (b)  projectiles  at 
78  *  7  m  s  the  shear  bands  propagate  very  close  to  the  surface  in 
(b). 
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Pressurised  PE  pipes  can  be  prone  to  failure  due  to  "Rapid  Crack  Propagation"  (RCP). 
Hence,  these  pipes  should  be  designed  from  tough  PE  grades  whose  dynamic  crack 
resistance  GD,  ensures  immediate  crack  arrest.  Quasistatic  fracture  calibration  tests  fail  to 
generate  RCP  in  these  tough  PE  grades.  Only  impact  tests  successfully  initiated  the  required 
mode  of  fast  fracture  in  precracked  PE  samples.  However,  these  tests  have  other 
drawbacks :  exact  monitoring  of  their  load  and  crack-length  traces  is  most  complicated; 
their  analysis  hinges  on  sophisticated  models  taking  into  account  dynamic  effects;  inherent 
transient  initiation  effects  and  gross  vibrations  may  cause  the  generated  Gp  values  to  be 
geometry  dependent  and  thus  cannot  be  applied  to  other  geometries  such  as  pipes. 

These  problems  pointed  out  the  need  for  an  alternative  test  which  imitates  the  pipes'  fast 
fracture,  suppresses  gross  dynamic  effects  and  is  still  easy  to  monitor  and  analyse.  Such  a 
test  has  been  developed  using  a  unique  Glass-PE  (G-PE)  duplex  sample :  A  polished  glass 
slide  is  adhered  to  the  crack-mouth  of  a  razor  prenotched,  side-grooved  PE  sample  (see 
Figure).  Quasistatic  tensile  loading  of  the  G-PE  sample  results  in  the  glass  shattering, 
followed  by  a  fast  and  brittle  fracture  in  the  PE.  Evidently,  prior  to  gla:.s  breaking,  the  glass 
suppresses  the  singularity  at  the  crack  tip  (which  remains  sharp),  while  large  amount  of 
elastic  energy  is  stored  in  the  sample.  When  the  glass  reaches  its  ultimate  load  it  breaks, 
rapidly  loading  the  sharp  crack.  A  brittle  crack  is  initiated  and  propelled  across  the  entire 
sample  by  the  released  elastic  energy. 

The  new  technique  has  been  successfully  applied  to  two  PE  pipe-grades :  BP  Rigidex  002- 
40,  and  Philips  TUB-71.  The  samples  were  clamped  between  heavy  grips  which  were 
bolted  to  the  frame  of  a  tensile  loading  machine,  Half  of  the  samples  were  precooled  to  0°  C 
while  the  others  were  kept  at  ambient  conditions  prior  to  loading  at  room  temperature. 
Load-displacement  curves  prior  to  RCP  were  plotted  on  the  tensile  machine  recorder.  Crack 
history  was  monitored  simultaneously  by  an  Imacon  image  converter  camera  and  a  digital 
oscilloscope  connected  to  timing  lines  painted  along  the  crack  path.  Both  camera  and 
oscilloscope  were  triggered  by  the  breaking  of  a  conducting  line  painted  on  the  glass  slide 
(see  figure  enclosed).  The  pictures,  showing  the  crack  history  either  in  framing  or  "streak" 
modes,  were  in  good  agreement  with  the  timing  lines'  readings.  The  average  velocities 
ranged  between  50  —  80  m  s’!  at  room  temperature  and  230  —  270  m  s-1  in  precooled 
samples.  It  should  be  noticed  that  the  only  effect  of  an  increase  in  the  remote  load  prior  to 
glass-break  seem  to  be  a  shorter  delay  between  glass-break  and  crack  initiation.  After  the 
initiation  effects  subsided,  the  crack  velocity  profiles  seemed  to  maintain  a  constant  value. 

Careful  SEM  study  of  the  fractured  surfaces  revealed  the  almost  ideal  uniform  nature  of  the 
surface  morphology.  There  was  hardly  any  evidence  of  a  different  micromechanism  at 
initiation  or  along  the  free  surfaces,  or  of  any  changes  along  the  crack  .  The  relevance  of 


this  small  scale  test  to  the  real  pipe  problem  wus  supported  by  their  similar  microfeaturea  i.e. 
a  brittle  "flaky"  surface  with  traces  of  drawn-out  films  and  fibrils. 

A  dynamic  finite  element  program  was  employed  to  process  the  experimental  data  into 
valid  crack  resistance  terms.  The  fixed  boundary  conditions,  the  glass  breaking  almost 
instantaneously  with  practically  no  energy  Iosb,  and  the  planar  uniform  crack  are  aU  easy  to 
define  in  2-D  numerical  terms.  The  2-D  elastic  code  simulates  crack  propagation  by  a 
successive  gradual  "release"  of  nodes  along  the  crack-path  and  applying  a  decaying 
"holding-back"  force  on  the  last  node  to  be  released.  Two  node  release  algorithms  were 
used  :  a  "propagation"  and  a  "generation"  mode.  In  the  "propagation"  mode,  the  nodes 
were  released  as  the  nodal  forces  reached  a  critical  value  (the  crack  was  made  to  propagate 
according  to  a  known  value  of  Gq).  In  a  "generation"  mode  the  nodes  were  released 
according  to  crack  history  data  and  the  Gq  values  were  calculated  from  the  global  energy 
balance. 

The  propagation  mode  analysis  predicted  that  for  a  constant  Gq  material,  the  magnitude  of 
the  remote  load  at  glass-break  should  dictate  the  time  of  initiation  and  crack  velocity.  The 
lower  the  load,  the  longer  it  should  take  to  initiate  the  crack  which  should  propagate  at  a 
lower  velocity.  These  predictions  contradicted  the  expetimental  observation  of  constant 
crack  velocities  with  no  regard  to  the  magnitude  of  the  remote  load  or  initiation  time. 
Generation  mode  analysis  of  constant  crack  velocity  data  showed  0D  values  to  vary 
according  to  initiation  time  :  an  early  initiation  resulted  in  lower  initial  value  of  Gq  which 
increased  as  the  crack  propagated  until  it  reached  a  higher  level.  Late  initiation  resulted  in  a 
higher  initial  GD  value  which  decreased  until  it  reached  a  constant  low  level. 

These  results,  indicating  a  non-unique  crack  resistance  -  crack  velocity  relation,  are  currently 
being  checked  by  conducting  additional  experiments  and  by  using  an  improved  FH 
viscoelastic  analysis. 


Schematic  view  of  the  G-PE  duplex  sample 
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1.  INTRODUCTION 

Sinca  tha  pioneering  work  of  Mann  raap.  Manjoina  and  Nadai 
rotating  disk  impact  machlnaa  hava  essentially  contributad  to 
tha  exploration  of  tha  behaviour  of  materials  under  dynamic 
loading  conditiona,  eapacially  for  atrain  rates  between  10 
1/a  and  2000  1/a.  Today's  decreasing  ,  interest  in  the  rotary 
impact  teat  technique  is  rather  astonishing,  aa  rotating  disk 
impact  machines  cover  the  upper  range  of  strain  rates  which 
are  of  relevance  in  civil  applications  and  can  be  used  for 
tensile,  compressive,  bending  and  fracture  mechanics  tests. 
The  paper  focuses  on  tensile  and  fracture  mechanics  tests. 
Besides  a  standard  machine  (disk  diameter  720  mm,  maximum 
rotational  energy  600  kJ),  MPA  Stuttgart  disposes  of  a  high 
energy  machine  with  a  disk  diameter  of  2000  mm  and  a  maximum 
rotational  energy  of  33  MJ  at  a  circumferential  velocity  of 
150  m/s.  When  the  specified  testing  conditions  are  estab¬ 
lished,  the  specimen  is  computer-controlled  slewed  into  the 
path  of  the  claws  by  a  fast  acting  pneumatic  cylinder.  The 
lower  bound  of  possible  specimen  dimensions  is  defined  by  the 
bending  stress  resulting  from  the  slewing  motion.  The  upper 
bound  is  established  by  the  distance  of  the  striking  claws. 
At  the  moment  this  is  equivalent  with  a  test  section  diameter 
ranging  from  6  to  20  mm  for  round  bars  and  a  thickness  of  10 
or  15  mm  for  compact  tension  ( CT)-specimens . 

2.  TENSILE  TESTS 

Little  attention  is  generally  paid  to  the  fact  that  an  ideal 
strain  rate  function  e  «  const,  cannot  be  realized  by 
servohydraulic  rigs  or  rotary  impact  machines.  For  the  latter 
the  strain  rate  can  only  be  controlled  via  two  parameters: 

-  the  circumferential  velocity  of  the  disk 

-  the  load-deformation  characteristic  of  the  damper  which  is 
situated  between  claw  and  anvil  in  order  to  couple  both  via 
a  semi-plastic  impact. 

This  characteristic  depends  strongly  on  the  geometry  and  the 
material  of  the  damper  and  defines  at  which  time  during  the 
deformation  process  the  specified  strain  rate  is  reached  and 
to  what  extent  disruptive  oscillations  resulting  from  the 
acceleration  of  specimen  and  anvil  are  present  at  the  load 
measuring  points.  This  represents  an  optimization  problem,  as 
a  long  acceleration  period  is  coupled  with  a  reduction  of  the 
oscillations.  The  formation  of  an  upper  and  lower  yield  point 
provokes  a  sudden  increase  of  the  strain  rate  exceeding  the 
specified  value,  whereas  in  the  beginning  strain  hardening 


range  a  rather  constant  strain  rats  can  bs  found.  Whilst  ths 
slastic  strain  rat#  is  dominated  by  ths  damper-character¬ 
istic,  ths  plastic  strain  rats  is  govarnad  primarily  by  ths 
circumfarantlal  velocity  of  ths  disk.  Materials  with  a 
continuous  transition  into  ths  plastic  range,  a.g.  high 
strength  i  steels,  do  not  exhibit  these  problems.  On  the  other 
hand,  a  strain  rate  in  excess  of  1000  1/s  cannot  be  applied 
to  low  toughness  materials  as  fracture  already  initiates 
during  the  acceleration  stage. 

Investigations  of  several  fine-grained  ferritic  stasis  have 

frovan  thii  applicability  of  the  following  test  procedure.  The 
oad  appjlied  is  measured  at  the  specimen  itself,  as  load 
cells  or  external  dynamometers  cannot  provide  reliable  load- 
time-graphs.  Two  dynamometer  sections  are  instrumented  with 
two  strain  gages  on  opposite  faces  to  compensate  for  bending 
effects.  The  strain  is  measured  via  post-yield  strain  gages 
which  ca n  be  used  up  to  a  temperature  of  300* C  and  provide 
strains  in  excess  of  20 %  if  properly  installed,  in  the  strain 
hardening  range  the  strain  can  also  be  deduced  from  the 
displacement  of  the  anvil.  For  that  purpose  a  bar  code  is 
attached  to  one  side  of  the  anvil  and  a  fast  acting  scanner 
is  installed. 

3.  FRACTURE  MECHANICS  APPLICATIONS 

Besides  double-edge  notched  tensile  specimens  CT-specimens 
are  suited  to  be  tested  in  rotary  impact  machines  in  accor¬ 
dance  with  ASTM  E  399/813.  This  requires  information  about 
the  load  and  the  load  line  displacement. 

According  to  a  proposal  by  Krabiell  the  load  is  measured  via 
strain  gages  which  are  attached  in  parallel  to  the  crack  path 
on  the  top  and  bottom  side  of  the  specimen  and  which  are 
calibrated  quasi-statically  before  the  dynamic  test.  With 
growing  plastic  deformation  an  increasing  deviation  from  the 
linear  interdependence  between  the  strain  gage  signal  and  the 
load  applied  can  be  observed.  These  calibration  curves  can  be 
evaluated  experimentally  for  static  loading  conditions  only. 
However,  this  is  not  useful,  because  in  this  case  the 
'''influence  of  the  dynamic  stress-strain-curve  cannot  be 
modelled.  Therefore  the  dynamic  calibration  curves  are 
calculated  with  the  finite-element-code  ABAQUS  and  dynamic 
streBs-strain-curves .  The  analyses  made  obvious  that  the 
strain  maximum  is  located  perpendicular  to  the  crack  tip  and 
that  other  strain  gage  positions  exhibit  a  linear  relation 
between  load  and  strain  gage  output  up  to  higher  loads. 

For  the  measurement  of  the  crack  opening  displacement  in  the 
load  line  clip  gages  or  other  conventional  inductive  or  capa¬ 
citive  displacement  gages  cannot  be  used  because  of  their  low 
natural  frequency.  Therefore  the  load  line  displacement  is 
measured  via  an  opto-electronic  sensor  based  on  the  scanner 
principle.  It  offers  the  advantages  of  low  additional  mass, 
high  bandwith  and  low  costs.  According  to  a  proposal  by 
Giovanola,  strain  gages  are  attached  to  the  crack  tip  near 
field  for  the  detection  of  crack  initiation.  Test  results  on 
fine-grained  ferritic  steels  and  ductile  cast  iron  are 
presented. 
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ABSTRACT:  A  high  speed  digital  video  camera  Is  developed  by  using  a 
new  MOS  type  solid  state  Image  sensor,  The  digital  Image  taken  by  the 
high  speed  video  camera  Is  analyzed  by  an  Image  processing  program  on  a 
personal  computer.  In  order  to  analyze  strain  rate  distribution,  a  new 
grid  method  using  Fourier  transform  Is  developed.  By  using  this 
method,  the  displacement,  velocity,  strain  and  strain  rate  distri¬ 
butions  of  an  unloading  wave  propagating  In  a  rubber  tube  are  analyzed. 


1.  HIGH  SPEED  VIDEO  CAMERA 

A  system  of  a  high  speed  video  camera  Is  made  with  a  new  MOS  type  solid 
state  Image  sensor  (Hitachi  Ltd.  HE98246,  Image  size:  649x491  pixels).  It 
Is  possible  to  change  the  area  size  for  output  into  any  size  using  soft¬ 
ware.  As  the  selected  area  becomes  smaller,  the  frame  speed  becomes 
higher.  The  maximum  frame  speed  Is  about  300,000  (ideal  10,000.000) 
frames/sec.  In  this  case,  the  frame  size  Is  only  one  pixel.  The  frame 
speed  or  Image  size  of  this  video  camera  Is  controlled  by  a  personal 
computer.  The  signal  from  the  vertical  output  line  is  stored  In  IC 
memories  by  using  an  A/D  converter  or  stored  on  a  video  tape  recorder.  If 
the  Image  signal  is  sent  to  the  same  IC  memories  by  cyclically  rewriting, 
the  desired  Images  are  stored  by  sending  a  trigger  after  the  occurrence  of 
the  phenomenon. 

2.  PROGRAM  OF  IMAGE  PROCESSING  ON  MICROCOMPUTER  (PIMPOM) 

In  order  to  analyze  the  image  taken  by  a  video  camera,  a  wide  use  program 
on  a  personal  computer  (NEC  Corp.  PC-9801 VX,  CPU  80286)  has  been  devel¬ 
oped.  The  image  size  is  640x400  pixels  with  16  (4  bits)  of  4096  colors. 
It  is  not  necessary  to  use  any  special  hardware  for  Image  processing.  It 
has  a  lot  of  menus  such  as  region  set,  disk  access,  value  transform,  one 
and  two  dimensional  mask  filtering,  thinning,  geometric  transform,  printer 
output,  process  repeat,  graphics,  Image  Input  by  an  Image  scanner,  Fourier 
transform,  Hadamard  transform,  histogram  etc.  The  processing  speed  is 
very  high,  because  each  menu  is  written  by  a  machine  language  and  these 
menus  are  combined  by  a  BASIC  program,  For  example,  mask  and  geometric 
transform  Is  processed  In  about  6  seconds,  respectively.  Value  transform 
Is  processed  In  about  3  seconds.  This  processing  time  is  for  the  case 
that  the  Image  size  is  640x400  pixels.  The  processing  time  Is  In 
proportion  to  the  Inverse  of  the  area  size. 

3.  STRAIN  RATE  ANALYSIS  OF  UNLOADING  WAVE  PROPAGATING  IN  RUBBER  TUBE  BY 
USING  FOURIER  TRANSFORM 

This  high  speed  video  camera  is  applied  to  analysis  of  stress  wave 
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(#)  Original  Image  (b)  Fourier  spectra  of  Fig.  (a) 


(e)  Velocity  Distribution  (X)  (f)  Strain  Rate  Distribution  (1/sec) 

Fig.  1  Analysis  of  Unloading  Wave  Propagation  In  Rubber  Tube 


propagation  In  a  rubber  tube.  A  grating  is  painted  on  the  tube  surface. 
When  the  rubber  tube  is  stretched  and  cut,  an  unloading  stress  wave 
propagates  In  the  tube.  The  one  dimensional  behavior  of  the  unloading  wave 
Is  Imaged  by  the  high  speed  video  camera. 

In  order  to  analyze  t.he  displacement,  strain,  velocity  and  strain  rate 
distributions  of  the  tube,  new  moire  and  grid  methods  are  developed  by 
using  the  Fourier  transform  of  the  image  of  the  deformed  grating  on  a 
material.  The  basic  procedures  required  In  this  analysis  are  as  follows: 

(1)  Sampling  of  a  deformed  grating  by  a  TV  camera. (Fig,  1(a)) 

(2)  Discrete  Fourier  transformation  of  the  image. (Fig.  1(b)) 

(3)  Extraction  of  the  first  harmonic  of  the  spectra. 

(4)  Inverse  discrete  Fourier  transformation  of  the  first  harmonic. 
(Generation  of  the  complex  grating  pattern) 

(5)  Computation  of  arguments  of  the  complex  grating. 

(6)  Computation  of  displacement  from  the  argument. (Fig.  1(c)) 

(7)  Computation  of  strain  by  differentiating  the  displacement. (Fig.  1(d)) 

(8)  Computation  of  velocity  by  differentiating  the  displacement(Fig.  1(e)) 

(9)  Computation  of  strain  rate  distribution  by  differenciating  the 
velocity. (Fig.  1(f)) 


The  analysis  Is  completely  automated  by  digital  image  processing.  All  of 
the  laborious  and  subjective  procedures  required  in  the  conventional 
analysis  such  as  fringe  sign  determination,  fringe  ordering  and  fringe 
Interpolation  are  eliminated  permitting  objective,  fast,  and  accurate 
analysis. 
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Predictions  of  classical  and  recently  formulated  generalized  constitutive 
models  for  metallic  materials  have  been  up  to  now  scarcely  compared  with 
the  results  of  valid  multiaxlal  experiments.  We  refer  here  to  tension- 
tension  experiments  and  to  their  numerical  simulation,  realized  at  low  and 
medium  strain-rate  using  a  cruciform-specimen. 

The  experiments  have  oeen  performed  at  an  electromechanical  biaxial 
testing  device  with  strain-rates  range  [10“4,  10“' 3s ~ ^  as  well  as  at  a 
hydropneumatic  biaxial  testing  device  with  strain-rates  range  CIO"', 
1023s-1.  A  biaxial  cruciform  specimen  has  been  used  which  was  Instrumented 
by  a  photoetched  fine  grid  with  400  cross-prints  over  the  upper  side  of 
the  gauge  field  whose  displacements  have  been  recorded  by  perpendicular 
filming  during  tests.  From  the  bottom  side  at  the  centre  of  the  gauge 
field  a  thermocouple  as  well  as  a  three-fold  strain  gauge  rosette  have 
been  fastened,  A  code  used  either  total  strain  record  from  strain  gauges 
(up  to  2%)  or  from  grid  measurement  (from  2%  up  to  rupture;  as  well  as 
temperature  and  force  d.e.  stress)  measurement,  In  order  to  find  all 
components  of  plastic  strain  tensor  and  plastic  strain  rote  tensor  as 
functions  of  time,  as  measure  of  directionality  correctness  the 
experhuental  plastic  strain  rate  direction  In  X],X2  plane,  f.e.j 
•  • 

tga  -  Epn^p22 


has  been  compared  with  the  values: 

2TH'T22 


tg3t 


T1  rvT22 

T22=UTTf 


predicted  by  theories  of  Perjyna  (with  Wises  yield  function)  and  Cernocky- 
Krempl,  respectively,  where  Cp  plastic  strain  rate  tensor,  T  stress- 
tensor,  v  Poisson  coefficient. 

Owing  to  Inevitable  complexity  of  the  cruciform  specimen  used  In  the 
experiments  a  preliminary  finite  element  stress-strain/strain-rate 
analysis  has  been  done  for  a  number  of  proportional  strain-rate  controlled 
paths  using  the  computer  code  A13AQUS. 

The  finite  elements  calculations  predicted  a  homogeneous  strain 
distribution  over  the  constant  thickness  central  part  of  the  cruciform 
specimen.  These  predictions,  having  a  qualitative  character  due  to  the 
rough  approximation  of  the  material  models  Implemented  In  the  code,  have 
been  quantitatively  confirmed  by  the  experiments  at  low  and  medium  strain 
rate.  The  homogeneous  strain  distribution  over  the  constant  thickness 
gauge  section  during  the  whole  test  already  demonstrates  that  the  applied 
loading  is  uniformly  distributed  by  the  armfingers  of  the  cruciform 
specimen.  A  further  direct  verification  of  such  feature  has  been  obtained 
by  the  strain  measurements  on  the  armfingers,  which  remain  elastic,  during 
the  test.  The  armfingers  of  the  cruciform  specimen  distribute  uniformly 
the  applied  loading  to  the  constant  thickness  central  part  of  the 
specimen.  This  fact  and  the  homogeneous  strain  distribution  over  the 
surface  of  the  central  part  confirm  that  the  constant  thickness  central 


part  of  the  cruciform  specimen  deforms  In  a  homogeneous  tenslonal  plane 
stress  state.  Furthermore  the  experimental  set-up  permits  the  measurement 
of  the  components  of  stress  ant!  strain  tensors  along  the  principal 
directions.  Therefore  the  conditions  are  satisfied  for  the  correct  use  of 
an  equivalence  criterion  for  the  evaluation  of  the  material 
characteristics  from  biaxial  experiments.  The  Huber-Mlses  equivalence 
criterion  has  been  used  for  constructing  the  equivalent  stress-strain 
curves  which  have  been  determined  at  low  and  medium  strain-rate. 

The  equivalent  flow  curves  at  low  strain-rate  show  strain  hardening;  the 
uniaxial  and  the  equlolaxlal  flow  curves  are  practically  coincident  up  to 
large  strain  values  while  the  equivalent  flow  curves  corresponding  to  the 
other  straining  paths  lies  under  the  equlblaxlal  flow  curves  already  at 
small  strain  values. 

In  the  plane  of  the  stress  deviator,  the  Initial  yield  stresses, 
determined  both  from  temperature  minimum  and  from  0.2%  of  equivalent 
plastic  strain  definition,  practically  lie  on  Mlses  circles. 

Strain-rate  hardening  has  been  observed  In  case  of  equlblaxlal  straining 
at  medium  strain-rate. 

Divergence  has  been  observed  between  the  straining  directions  determined 
from  experiments  and  those  predicted  by  some  constitutive  models. 


Plastic  Flow  of  a  Ferritic  Mild  Steel  ana  a  High  Strength  Austenitic 
Steel  under  Dynamic  Biaxial  Loading 


K.  Stiebler,  H.-D.  Kunze,  E.  Staskewitsch 

Fraunhofer-Institut  fUr  angewandte  Materialforschung,  I FAM , 
Lesumer  Heerstr.  36,  D-2820  Bremen  77,  W.  Germany 


The  increasing  use  of  computers  in  the  area  of  technical  design  and  optimal 
utilisation  of  materials  require  a  better  understanding  of  the  material 
behaviour  at  its  specific  loading  types  and  In  addition  its  mathematical 
discription.  Especially  in  the  field  of  multiaxial  loading  under  high  rates 
of  deformation  were  only  limited  data  available.  To  expand  this  data  base 
two  testing  rigs  were  designed  and  built  to  investigate  materials  under 
these  loading  conditions. 

With  a  hydraulic  apparatus,  thin-walled  tubular  specimen  were  simultaneous¬ 
ly  loaded  by  two  normal  stresses  up  to  elastic  strain  rates  in  the  range  of 
10-1  s"l.  The  normal  stresses  wore  combinations  of  either  tension/compres¬ 
sion  and  internal/external  pressure.  Higher  strain  rates  of  102  s_1  were 
reached  with  a  modified  split.  Hopkinson  bar  which  allowed  simultaneous  com¬ 
bined  tension  and  torsion  stresses,  to  be  achieved. 

Both  experimental  set-ups  were  used  for  the  quasistatic  and  dynamic  inves¬ 
tigation  on  Ck  35  ferritic  mild  steel.  For  this  material,  yielding  is  de¬ 
fined  by  the  lower  yield  point.  Results  of  the  tests  with  two  normal  stres¬ 
ses  in  axial  and  circumferential  directions  show  for  both  elastic  strain 
rates  ip  =  10_5  s“l  and  te  -  10_1  s-1,  a  very  good  agreement  in  the  first 
quadrant  between  the  experimental  yield  points  and  the  ellipses  predicted 
by  the  von  Mlses  yield  criterion,  fig.  1.  The  results  in  the  second  auad- 
rant  differ  somewhat  from  this  criterion.  Under  combined  tension  and  tor¬ 
sion  loading  the  yield  points  fulfill  the  Tresca  yield  criterion  for  both 
the  quasistatic  (ee  -  2 * 10~5  s’l)  and  the  dynamic  Ue  =  3 * 10l  s'1)  biaxial 
tests,  fig.  2.  Consequently,  the  yield  criterion  changes  with  the  nature 
of  load.  The  yield  points  of  the  dynamic  tests  showed  higher  stresses  than 
the  quasistatic  one. 

The  yield  loci  of  the  quasistatic  and  dynamic  combined  tension-torsion- 
tests  of  an  austenitic  steel  are  described  by  ellipses  which  correspond  to 
different  proof  strain  definitions  of  yielding,  fib.  3.  This  figure  shows, 
that  the  size  imd  the  shape  of  the  yield  loci  are  "functions  of  strain  and 
strain  rate.  For  low  plastic  strains  and  plastic  strain  rates  of 
iP  •  10"3  s"1  the  yield  points  are  in  a  good  agreement  with  the  von  Mises 
yield  criterion.  Whereas  for  large  strain  and  strain  rates  UP  3  102  s”l) 
the  yield  loci  tend  towards  the  Tresca  yield  criterion. 


The  results  of  these  investigations  on  the  behaviour  of  materials  were 
used  in  the  calculations  to  allow  realistic  simulations.  For  the  descrip¬ 
tion  of  the  quasistatic  and  dynamic  tension-torsion  behaviour  of  the  auste¬ 
nitic  steel,  Perzyna's  constitutive  equation  was  choosen.  It  includes  a 
yield  criterion  as  a  function  of  strain  and  strain  rate  and  a  function 
which  incorporates  the  dependence  of  flow  stress  versus  strain  rate  based 
on  thermal  activation.  For  specific  measured  strains  and  strain  rates  the 
solution  of  the  constitutive  equation  results  in  the  associated  stresses. 
These  calculated  stresses  are  In  a  good  agreement  with  the  measured  stres¬ 
ses  for  all:  uniaxial  and  biaxial  tests. 


Fig.  1. 

Yield  loci  for  Ck  35  steel  loaded 
by  two  direct  stresses  at  two  strain 
rates. 
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Fig.  2. 

Yield  loci  for  Ck  35  steel  loaded 
by  tension-torsion  at  two  strain 
rates. 
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Fig.  3 

Yield  loci  of  the  quasistatic  and 
dynamic  combined  tension-torsion- 
tests  on  austenitic  steel  corres¬ 
ponding  to  different  proof  strain 
definitions  from  e P  =  0.2%  to  0.3%. 
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The  Ueo  of  Microscopically  Based  Constitutive  Models  in  Explosive  Metal 
Deformation  Modelling  Studios 


P  Church 
I  Cullis 


ABSTRACT 

The  numerical  modelling  of'  warheads  and  their  terminal  effects 
requires  an  accurate  description  of  material  behaviour  under  the  varying 
load  conditions  associated  with  shock  waves,  and  the  resulting 
hydrodynamic  and  plastic  deformation.  Choir  accuracy,  however,  is 
limited  by  the  constitutive  relation  and  thin  ability  of  the  equation 
of  state  to  predict  the  physical  state  variables  associated  with  these 
complex  stress  conditions. 

There  are  3  main  mechanisms  which  affect  the  way  a  material  deforms 
under  high  strain  rate  explosive  loading.  Those  are  work-hardening 
effects,  thermal  softening  effects  and  shock  hardening.  They  are  all 
influenced  by  grain  size,  strain  rate,  phase  changes,  and  impurity 
concentrations.  Ultimately,  therefore,  the  constitutive-  model  must 
be  able  to  describe  dynamic  path  dependant  plastic  deformation  over 
a  wide  range  of  temperatures  and  material  conditions. 

Recent  work  by  Pol  Jans  bee  on  copper  (19116)  end  Arms  l  rong-Xerllli  on 
iron  and  copper  (1986)  has  attempted  to  describe  material,  behaviour 
at  high  strain  rat.e  using  dislocation  mecnanios.  the  Armstrong-Zerilli 
model  derives  flow  stress  as  a  function  of  strain,  l.ormpefature ,  strain 
rate  and  grain  size  and  distinguishes  between  face-centred  cubic  (eg 
copper)  and  body  centred  cubic  materials  (eg  iron).  Account  is  also 
taken  of  the  influence  of  solute  material  and  the  original  dislocation 
density. 

The  Armstrong-Zer-il  li  ronstitutve  model  has  been  in  Explosively  Formed 
Projectile  (EFP)  calculations  and  compared  with  a  semi-empirical 
constitutive  model,  where  high  strain  rate  behaviour  is  described  by 
a  semi-empirical  work-hardening  curve  and  a  thermal  softening  curve. 

Two  EFP  designs,  using  an  Armco  iron  liner,  based  on  different  performance 
military  explosives  (Onto!  and  Composition  B),  have  been  modelled  on 
the  0YNA2D  hydrocode.  The  semi-empirical  model  and  lue  Armstrong- 
Zerilli  model  with  the  correct  grain  size  (2(3, urn)  give  the  same  order 
of  agreement  in  terms  of  describing  the  experimentally  observed  EFP 
projectile  profiles.  Changing  the  grain  size  to  l20^um  gives  better 
agreement  for  the  Armstrong-Zeri Hi  model  for  both  designs. 


For  the  , semi -empirical  model  the  detormi'iution  ol'  the  thermal  behaviour 
of  the  flow  stress  is  restricted  to  a  simple  comparison  with  an  experimental 
profile.  Sides  the  EFP  is  crucially  dependent  on  this  thermal  behaviour 
of  the  flow  stress,  this  represents  a  serious  drawback  if  a  predictive 
capability  is  required.  Therefore  the  semi-empirical  model  can  only 
be  oonsidered  a  first  oraer  material vftddei.;'" 

The  theoretical  modnls  do  try  and  account  directly  for1 the  complex  material 
mechanisms  by  dislocation  mechanics,  though  at  present  no  description 
is  provided  for  shock  hardening,  twinning  nr  rocrystallieation.  The 
results  so'  far  give  encouragomont  in  trying  to  describe  the1  complex 
internal  mechanisms  in  EPPs.  This  is  vital  if  one  ie  trying  to  understand 
the  sensitivity  of  material  properties  in  defining  design  envelopes. 

In'  conclusion,  therefore,  the  semi-empirical  models  are  op.ly  adequate 
for  first  order  effects  (experimental  profiled,  velocities),  whereas 
the  theoretical  type  models  show  great  promise  in  also' understanding 
second  order  effects  associated  with  tho  complex  internal  deformation 
mechanisms  in  EFPs.  More  work  is  required  in  refining  those  theoretical 
type  models.  '  , 
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Analysis  of  tha  Strain-Rats  Sonaitivtty  at  High  Strain  Rstsa  in  FCC  and 
BCC  Ms  tala 


Paul  S,  Follanabe* 

Los  Alamos  National  Laboratory,  Los  Alamos,  NM,  87545,  USA 


A  constitutive  formalism  for  deformation  over  a  wida  range  of  strain 
rates,  temperatures ;  and  strains  has  been  applied  to  several  pure  ftc 
metals  and  to  sevaral  fee  alloys.  The  formalism  is  leased  on  the  use  of 
internal  state  variables  and  a  phenomenological  understanding  of  the 
kinetics  of  dislocation  glide,  dislocation  storage,  and  dynamic  recovery- 
In  the  mechanical  threshold  stress  model,  the  applied  stress  (the  yield 
stress)  required  for  deformation  in  a  given  state  is  written  as  a 
combination  of  an  athermal  component,  o  ,  and  several  thermally  activated 
components  es 

'<*  ;  *fl)r  -  £[si<‘*T>  *i  ]t  •  (!) 

where  a.  gives  the  mechanical  threshold  stress  that  characterises  the 
interaction  of  dislocations  With  obstacle  i,  s^  specifies  the  kinetics  of 
those  interactions,  n  is  the  number  of  separate  obstacle  contributions, 
and  r  is  a  power  that  accounts  for  the  interactions  between  different 
types  of  obstacles  (Kr<2). 

Experimental  techniques  have  been  established  to  measure  the  internal 
state  variables  and  the  activation  energies  that:  characterize  disloca- 
t  iori/obstacle  interactions.  These,  techniques  involve  measurement  of  the 
temp',  rature  and  strain-rate  dependent:  yield  stress  in  samples  deformed 
according  tc  various  strain  rate,  strain,  and  temperature  histories. 
Because  we  are  mostly  interested  in  the  behavior  at:  high  strain  rates,  we 
restrict  our  yield  stress  measurements  to  poom  temperature  and  below. 

In  pure  fee  metals,  dislocation/dislocation  interactions  provide  the  only 
thermally  activated  contribution  to  the  fiow  stress,  in  addition  to  the 
small  grain  boundary  component,  which  is  considered  to  be  athermal.  This 
simplifies  the  analysis  in  that  n-1  and  thus  r-l  in  Eq .  (1).  The  study  of 
fee  alloys  introduces  the  complexity  of  n>l  in  Eq.  (1).  In  Ni-C  alloys, 
for  example,  carbon  goes  to  interstitial  sites  which  yields  a  more  ratc- 
dependent  strengthening  contribution  than  found  in  pure  fee  metals,  In 
austenitic  stainless  steels  the  addition  of  interstitial  atoms  (e.g, ,  C, 
N,  0)  and  substitutional  atoms  (e.g,,  Cr,  Ni,  Al)  leads  to  strongly  rate- 
dependent  deformation.  Dislocation  motion  in  a  bcc  lattice  is  opposed  by  a 
large  Peierls  barrier  which  leads  to  the  strongly  rate  dependent  yield 
stress  observed  in  metals  with  this  crystal  structure.  The  generation  and 
storage  of  dislocations  through  strain  hardening  then  gives  a  second 
contribution  to  the  flow  streas  in  pure  bcc  metals,  Thus,  even  pure  bcc 
metals  are  complicated  by  n-2  in  Eq.  (1).  The  combined  contribution  of 
several  strengthening  mechanisms  complicates  the  analysis.  However,  we 


can  begin  by  defining  two  cUmu:  those  contributions  ■  that  evolve  with 
utreivi  and  tjhoaa  that  are  preaoftt  in  the  annealed  condition  and  which  are 
assumed  to  remain  constant  with  strain.  Tha  former  contribution  ia 
considered  to  he  due  to  eho  stored  dislocation  donalty. 


Strain  hardening  leads  to  an  evolution  of  tha  state, 


an  increasing 
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(Uslocatlrtrv  dbnslty  Kvolvtlon  is  treated  .differentially  according  to  ! 
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where  o.  is  tha  mechanical  threshold  stress  that  characterizes  dialoca- 
tivm/dillocaeion  interactions,  3  is  the  temperature  and  strain-rate 
dependent:  aneuratiOn  value  of  3  ,  |nd  t  is  the  Stage  II  hardening  rate. 
The  Stage  II  hardening  rate  had  previously  boau  considered  to  be  a  con¬ 
stant  (of  order  m/20),  but  our  work  in  copper  (Follansbee  and  Kocka,  1988) 
loti  to  tho  conclusion  that  whan  the  strain  rate  is  raised  above  »10*  a  *.'* , 
tha  Stage  II  hardening  rate  becomes  strongly  strain-rate  dependent,  which 
results  in  rapid  hardening  and  high  flow  stresses .  While  the  origin  of  the 
strain-rate  dependence  of  f  remain*  open,  a  Contribution  of  .this  previous 
worh  is  the  observation  “that  the  behavior  at  high  strain  rates  is 
dominated  by  this  term,  rather  than  by  s(i,T)  in  Eq.  (1)  or  6  (i,T)  In 
Eq,  (2),  ** 


The  rapid  dislocation  generation  rates  observed  at  high  strain  rates  in 
copper  and  nickel  may  be  a  general  phenomenon  in' fee  metals,  although 
there  is  evidence  In  Nittonio  40  stainless  steel  that  deformation  twinning 
can  offset;  this  affect.  Work  in  bee  metals  is  just  beginning,  Similar 
trends  to  those  found  in  fee  metals  are  observed  in  iron  and  A3A0  steel. 
In  pure  iron  we  ohserve  that  the  rate  of  evolution  is  actually  lower  in  a 
sample  deformed  dynamically  than  in  a  sample  deformed  slowly. 
Metaliographic  examination  of  these  samples  showed  the  presence  of 
deformation  twins  in  the  former  but  not  in  the  latter.  The  description  of 
deformation  twinning,  particularly  the  associated  kinetics,  the  influence 
on  strain  hardening,  and  the  balance  between  deformation  by  twinning 
versus  that  by  glide represents  a  challenge  for  a  complete  deformation 
theory  for  the  high  strain  rate  regime. 
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CONSTITUTIVE  PROPERTIES  OF  COPPER  AND  TANTALUM 
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Tha  electromagnetic  launch  technique  [1]  has  been  used  to  study  the 
properties  of  oxygen  free  electronic  (QFE)  grade  copper  and  99. 9%  tantalum 
rings  expanded  at  peak  strain  rates  of  Urs“i.  Ouring  the  Tree 
expansion  phase,  measurements  of  tha  ring  deceleration  using  a  velocity 
interferometer  (VISAR)  yield  the  flow  stress,  strain,  and  strain-rate 
histories.  Specimen  temperature  is  inferred  from  launch  currants. 

The  results  for  OFE  copper,  processed  to  yield  uniform  grain  sizes 
between  10  ym  and  150-200  ym  with  minimal  texture,  show  a  sample-to- 
sample  variation  of  ±20  MPa.  Experiments  with  rolled  material,  for 
which  the  initial  amount  of  cold  working  is  uniform,  indicate  that 
approximately  half  of  this  uncertainty  can  be  ascribed  to  sampla-to-sample 
variations  in  the  mechanical  properties  produced  by  non-uniformities  in 
the  working  prior  to  final  heat-treatment.  The  remainder  apparently 
represents  a  limit  imposed  by  the  noise,  electronic  and  physical,  that  is 
characteristic  of  measurements  of  the  speed  of  the  external  surface  of 
the  specimen  ring  using  the  VISAR.  Despite  these  uncertainties,  however, 
the  constitutive  properties  of  the  10  ym  and  L50-200  ym  materials  are 
clearly  distinguished.  While  the  flow  stress  of  10  ym  OFE  copper 
increases  from  280  MPa  to  «C0  MPa  for  strains  between  13%  and  50%  and 
strain-rates  between  about  8000  s"l  and  3000  s_J-,  the  stresses  in 
150-200  ym  material  fall  75-100  MPa  lower  iri  the  same  ranges  of  strain 
and  strain-rate.  These  data  agree  well  with  Hopkinson  bar  data  at 
smaller  strains,  confirming  the  apparent  effect  of  grain  size.  Simple 
models  in  the  literature  [2,3]  do  not  provide  an  adequate  description  of 
the  ring  data,  and  comparisons  to  more  recent  dislocation  based  models 
[4,5]  are  presented. 


•Work  performed  under  the  auspices  of  the  U.S.  Department  of  Energy  by 
the  Lawrence  Livermore  National  Laboratory  under  contract  number 
W-7405-ENG-48. 
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Tantalum  specimens,  because  of  their  low  electrical  conductivity, 
must  be  launched  with  high  conductivity  "pusher"  rings  that  are  arrested 
after  ring  acceleration.  Ring  data  show  virtually  no  work  hardening  and 
limited  ductility  (30%  natural  strain  to  failure)  for  99,9%  tantalum 
having  a  grain  size  of  approximately  100  urn.  The  constitutive 
properties  obtained  from  ring  tests  are  compared  with  the  data  and  model 
of  Hoge  and  Mukharjae  [6]. 
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Strain  rate  dapandanoa  of  mechanical  bahaviour  and  ovaluation 
of  tha  microstructure  of  tungsten  monocrystala 

L.W.  Meyer,  IFAM-Bremen,  W. -Germany 
and 

C.y.  Chism,  ENSM-Nantea,  Frances 


Abstract 

Tungsten  monO'cryatais  wera  compressed  in  ( 110  ) -directions  at 
three  loading  rates  of  i  *  10"4,  10°  and  10“*  s"1  to  investigate 
tha  strength,  deformability,  fracture  mechanism  and  the  micro¬ 
structure.  Tha  vary  early  transition  from  elastic  to  plastic 
deformation  bahaviour  is  characteristic  for  the  behaviour  of 
tungsten  singla  crystals.  As  a  function  of  loading  velocity 
this  transition  is  reached  between  20  and  100  MPa. Initially 
high  deformation  hardening  is  achieved  up  to  3  %  deformation 
resulting  in  a  flow  stress  of  1000  MPa.  Increase  in  strength 
due  to  six  orders  of  magnitude  increase  in  strain  rate  results 

in  a  flow  stress  increase  of  about  300  MPa.  With  a  compression 
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strength  of  more  than  1500  MPa  at  a  strain  rate  of  10  s  the 
level  of  flow  stress  is  remarkable.  This  is  similar  to  that 
of  the  sintered  tungsten  heavy  metal  alloys. 

Transmission  electron  microscopy-discs  parallel  to  the  (Oil) -plane 
ware  taken  in  the  mid  of  the  specimen,  compressed  at  12  -  16  % 
reduction  of  height,  to  investigate  the  variation  of  the  dis- 
location  structure  as  a  function  of  the  strain  rate.  At  £  =  10  s 

the  dislocation  distribution  in  the  (Oil)  plane  is  very  uniform 
with  numerous  nearly  straight  screw  dislocations  parallel  to  the 
< 111 >  directions .  With  increased  rate  of  strain  to C  =  10°  s  *  in 
addition  to  the  screw  dislocations  edge  and  mixed  dislocations 
were  observed.  Also  at  this  strain  rate  many  of  the  dislocations 
were  tangled.  At  C  =  10  s  a  cell  structure  is  starting  to  form 
in  the  (Oil)  plane,  but  the  main  orientation  of  the  dislocations 
are  still  the  < 1 1 1 >  directions.  Dislocation  density  increases 
with  increasing  strain  rate.  Comparing  the  increase  of  the  flow 
stress  (due  to  the  strain  rate)  with  the  dislocation  density 
change  a  logarithmic  dependence  is  given. 


In  contrast  to  sintered  tungsten  alloys  under  dynamic  loading/ 
the  monocrystals  fail  at  12  to  14  %  reduction  of  height  by 
shear  failure.  The  shearing  occurs  in  a  planar  area  that  is 
diagonal  to  the  loading  direction.  It  ii  from  this  orientation 
that  the  TEM  specimen  were  taken.  In  the  following  loading/ 
high  stresses  in  the  radial  direction  arise  in  the  sheared/ 
fractured  segments.  The  monocrystal  is  not  able  to  sustain 
these  stresses  and  a  cleavage  crack  network  occurs  parallel 
to  the  main  crystallographic  planes. 


A  Constitutive  Description  of  the  Deformation  of  Alpha  Uranium  Baaed  on 
the  Use  of  the  Mechanical  Threshold  Strata  as  a  State  Variable 


P.  E,  Armstrong,  P.  S,  Follansbea,  T.  Zocco 

Los  Alamos  National  Laboratory,  Los  Alamos,  New  Mexico,  87543,  USA 


The  deformation  of  isotropic,  alpha  uranium  over  a  strain-rate  range  of 
10’4  s'1  to  104  s’*,  a  temperature  range  of  78K  to  900K,  and  a  strain 
range  to  0.6  is  modeled  using  a  constitutive  description  based  on  the  use 
of  evolving  internal  state  variables.  The  model  ia  able  to  fit  the 
experimental  results  and  data  from  previous  measurements,  although  it  does 
not  include  a  specific  contribution  from  deformation  twinning.  The  model 
assumes  thAt  for  a  given  atruoture,  the  Instantaneous  yield  stress  is  a 
function  of  tha  currant  temperature,  T,  strain  rate  t,  and  mechanical 
threshold  stress,  or 


3  n 

T  +  l 


(1) 


where  p  is  the  temperature  dependent  polycrystalline  shear  modulus,  o  is 
the  mechanical  threshold  stress  characterizing  athermal  dislocations 
interactions  and  o ^  is  the  mechanical  threshold  stress  characterizing 
thermally  activated  interactions.  The  functions  describe  the  tempera¬ 
ture  and  strain- rate  dependencies  of  each  of  these  thermally  activated 
interactions.  We  have  used  the  following  form  of  the  function  s  that  has 
been  applied  to  dlslocation/obstaclo  interactions; 


(2) 


where  k  is  the  Boltzmann  constant  (1 . 381x10" 29  MNm/K)  ,  b  is  the  Burgers 
vector,  g  is  the  total  normalized  activation  energy,  e  is  a  constant 
(typically  107  s‘*<  c  <1010  s'1),  and  p  and  q  are  constants  ( 0<p< 1 ; 
l<q<2).  To  evaluate  the  increase  in  mechanical  threshold  stress  with 
strain  (strain  hardening)  we  write 


(3) 


where  6  is  the  differential  increase  in  the  mechanical  threshold  stress, 
$  is  the  initial  (or  Stage  II)  strain-hardening  rate,  and  o_  is  the 
saturation  value  of 
ing  rate  is  usually 

obeyed  when  F-l  in  Eq .  (3),  Key  to  Eq .  (3)  is  the 
state,  represented  by  a  ,  which  is  the  temperature  and  strain-rate 
dependent  maximum  threshold  stress  reached  at  large  strains.  Because  the 
initial  dislocation  density  was  observed  to  be  very  low,  we  separate  the 
mechanical  threshold  stress  into  three  components, 


s 

the  mechanical  threshold  stress.  The  Stage  II  harden- 
a  constant,  roughly  equal  to  a/20.  The  Voce  law  is 

concept  of  a  saturation 


d  -  +  dt  ,  (4) 

whvr«  £-  is  the  Initial  value  (S.-600  MPa),  presumably  due  to  the  Pcierls 
barrier,  and  6  is  the  contribution  due  to  the  dislocation  density  which 
is  evolving.  Equation  (1)  therefore  becomes 

'J  o 

e(j ,T)  [S£  +  S£.|  ,  ,  (5) 

where  s  is  again  given  by  Eq,  ( 2)'.  In  Eq.  (4)  and  (5)  only  the  o  term 
needs  to  be  described  using  Eq.  (3).  An  integrated  form  of  Eq,  (3)  was 
used  with  the  further  assumption  of  Voce  law  behavior  < F— 1 > ,  giving 


The  temperature  and  strain-rate  dependencies  of  a  and  9  were  described 
with  an  equation  chat  had  been  used  previously  £n  foe  metals  to  describe 
the  temperature  and  strain- rata  dependency  of  d  , 

In  <d  )  -  In  (d  )  +  rr  In  (  4- )  ,  (7) 

CS  SO  A  £ 

C 

where  t  ,  A  and  o  are  constants.  This  equation  has  a  physical  inter¬ 
pretation  in  desc??bing  but  Its  application  to  as  in 

In  <#  >  -  In  (#0>)  4  k-|-ln  (4-)  (8) 

s 

is  purely  empirical. 

The  parameters  in  the  above  equations  were  determined  from  a  fit  to  200 
experimental  results  (yield  stress  as  a  function  of  temperature  and  strain 
rate  on  samples  prestrained  according  to  various  strain-rate,  strain,  and 
temperature  histories) .  Although  the  model  can  be  made  to  agree  with  the 
experimental  data,  its  physical  basis  in  uranium,  where  twinning  is  an 
important  deformation  mechanism,  is  questioned. 

Deformation  microstructures  are  characterized  using  optical  and 
transmission  electron  microscopy.  The  general  observation  is  that 
deformation  initially  occurs  primarily  by  deformation  twinning,  but  that 
with  increasing  strain  there  is  a  transition  to  deformation  primarily  by 
slip.  The  effect  of  this  transition  on  the  proposed  constitutive 
formulation  is  discussed  and  compared  with  the  results  ol  the  analysis. 
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Localization  Analysis  under  Dynamic  Loading 


Y.  Leroy*  an<J  M.  Ortiz 

Brown  University,  Providence,  RI  02912  USA. 


A  finite  element  method  proposed  by  Ortiz  et  al.  (1987)  is  used  to  study  shear 
band  formation  in  rate  independent  and  rate  dependent  pressure  sensitive  solids 
under  dynamic  loading. 

This  new  method  is  characterized  by  the  addition  to  the  strain  interpolation  of 
ar,  extra  mode  of  deformation  carrying  a  suitable  strain  discontinuity.  This 
added  mode  of  deformation  permits  to  alleviate  the  overly  stiff  response  of 
isoparametric  elements  observed  in  problems  involving  strain  localization. 

Results  concerning  localization  criteria  are  first  reviewed.  These  conditions  arc 
pertinent  to  the  set  up  of  the  added  modes  in  the  finite  element  method  used  here 
as  well  as  to  the  design  of  crossed  triangular  meshes. 

For  a  rate  independent  solid,  the  analysis  is  based  on  the  acoustic  tensor  defined 
for  Hill's  comparison  solid  (1958).  Localization  is  possible  as  soon  as  a  wave 
velocity  is  found  to  vanish  along  some  direction  which  defines  the  normal  to  the 
incipient  discontinuity.  The  wave  polarisation  vector  determines  the  mode  of 
localization.  For  rate  dependent  materials  wave  velocities  always  remain  real. 
Nevertheless,  localization  is  observed  for  such  material  in  the  form  of  unstable 
growth  of,  for  example,  shearing  mode  of  deformation,  Marchand  and  Duffy 
0987). 

Using  a  linearized  instability  analysis,  Molinari  (1987),  Leroy  and  Ortiz  (1989) 
have  determined  necessary  conditions  for  the  earliest  possible  localization 
instabilities.  The  main  outcome  of  such  analysis  is  that  the  first  possible 
instabilities,  i.e.  those  which  grow  infinitely  slowly,  can  be  recovered  by 
considering  the  inviscid  limit  of  the  constitutive  model. 

As  a  first  example  a  single  clement  problem  is  presented  in  order  to  motivate  the 
necessity  of  a  specialized  element  when  dealing  witli  localization.  Delay  in  the 
localization  process  when  rate  sensivity  of  the  solid  is  increased  is  also  observed. 
Next  we  turn  our  attention  to  the  dynamic  response  of  frictional  solids.  The 
constitutive  description  adopted  in  the  calculations  is  of  the  Drucker-Pragcr 
type  with  a  monotonic  hardening  law  on  the  friction  angle  up  to  a  saturation 
level.  The  dilatancy  angle  is  taken  to  be  null  throughout  the  analysis,  rendering 
the  plastic  response  non-associative.  Wc  consider  a  rectangular  sample  of 
material  constrained  to  undergo  plane  strain  compression.  Velocities  are 
prescribed  on  its  upper  surface  following  a  ramp  variation  in  time.  The  final 
impact  velocity  is  of  the  order  of  1/1000  of  the  elastic  longitudinal  wave  speed. 

*  Presently  at  Koninklijkc/Shell  Exploratie  en  Produktic  Laboratorium,  The 
Netherlands. 


Several  tests  are  conducted  for  a  rate  independent  model  and  a  rate  dependent 
model  based  on  a  linear  overstress  viscous  law.  These  results  are  compared  with 
the  static  and  quasi-static  analysis. 

The  most  striking  diiTcrcnce  between  the  quasi-static  and  the  dynamic  solutions 
concern  the  geometry  of  the  shear  band.  Whereas  under  quasi-static  loading  a 

single  shear  band  emerges  in  each  quarter  of  the  specimen,  a  network  of  bands 

is  seen  to  arise  under  dynamic  loading.  Shown  in  fig.  1  is  the  distribution  of 
effective  platic  strain  for  the  dynamic  rate  independent  case.  The  introduction 
of  rate  sensivity  has  several  effects  on  the  solution.  Substantial  delay  in  the 
localization  time  is  observed  as  well  as  a  broadening  of  the  shear  band  away 
from  its  nucleus. 

A  noteworthy  outcome  of  these  computations  is  the  fact  that  the  geometry  of 
shear  bands,  both  as  regards  their  thickness  and  spatial  distribution  is  quite 

sensitive  to  rate  dependence  and  inertia  effects.  These  results  raise  questions  on 
the  prospects  for  describing  the  thickness  of  shear  band  by  means  of  a  single 
parameter,  a  feature  common  to  many  non-local  and  generalized  continuum 

theories. 
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Fig.  1:  distribution  of  effective 
plastic  strain  for  the  rate 
independent  dynamic  test. 
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This  paper  describes  the  development  and  application  of  a  computer  program 
for  the  study  of  materials  response  to  short-duration  loading,  i.e., 
situations  where  wave  motion  in  the  material  due  to  impact,  explosive 
loading  or  energy  deposition  must  be  accurately  tracked. 

Wave  propagation  codes  were  originally  developed  to  solve  problems  charac¬ 
terized  by: 

a.  the  presence  of  shock  waves  (alternatively,  steep  stress  or  velocity 
gradients) 

b.  Localized  materials  response  (i.e.,  situations  where  the  overall  geo¬ 
metric  configuration  of  a  structure  is  of  secondary  importance  compared  to 
the  constitution  and  characteristics  of  the  material  in  the  vicinity  of 
the  applied  Load) 

c.  loading  and  response  times  in  the  submillisecond  regime. 

Characteristics  of  computer  codes  to  solve  these  problems  have  been  dis¬ 
cussed  by  Zukas  (1CI82,  1187)  among  others.  Excellent  results  have  been 
obtained  with  wave  codes  involving  many  different  materials  in  engineering 
designs  with  very  complex  shapes.  Good  accuracy  lias  been  achieved,  in 
comparison  to  exact  solutions  or  with  experimental  results  when  the 
materials  involved  ait:  very  well,  known  and  characterized.  When  this  is 
not  done,  when  the  codes  are  used  as  black  boxes  with  engineering  mat¬ 
erials  which  are  not  weJ  I.  understood  or  characterized,  it  is  not  uncommon 
to  obtain  qualitatively  incorrect  solutions,  i.e.,  garbage  in,  garbage  out.. 
The  greatest  limitations  on  accuracy  in  wave  propagation  calculations  are 
the  material  descriptions  embodied  in  t lie  constitutive  equations,  particu¬ 
larly  with  material  failure.  Equally  important  is  the  requirement  that 
the  material  data  (the  various  coefficients  of  the  constitutive  description) 
must  be  appropriate  to  the  characteristic  time  scale  of  the  problem  studied, 


The;  ZeuS  code  is  a  2-1)  explic  it  Lagrangiun  finite  dement,  code  originally 
designed  to  operate  in  the  PC)  or  MS-DOS  environment.  Because  of  its  ex¬ 
plicit  integration  scheme,  it  is  best,  suited  for  high  velocity  impact 
situations  although  nothing  in  the  physics  ol  the  formulation  inhibits  Its 
use  for  low  velocity  impacts.  As  a  general  rule,  though,  problems  which 
fall  into  the  structural  dynamics  category  -  situations  where  the  response 
is  measured  in  a  time  frame  of  milliseconds  to  seconds  and  where  both  the 
geometry  of  the  structure  as  well  as  the  material  constitution  must  be 
considered  -  are  best  addressed  with  codes  employing  implicit  integration 
schemes.  A  ZeuS  simulation  can  model  objects  involving  up  to  ten  isotropic 
inert  and  five  explosive  materials  simultaneously.  The  expandable  material 
library  maintained  by  the  ZeuS  Interactive  pre-processor  permits  material 
characteristics  to  be  stored  for  a  virtually  limitless  number  of  materials. 
Within  the  640  kbyte  environment,  of  MS-DOS,  ZeuS  permits  simuiat.ions  in¬ 
volving  2500  nodes  and  4500  elements.  The  node  and  element  state  variables 
are  contained  completely  in  RAM  memory,  thereby  eliminating  the  need  to 
engage  in  time  consuming,  inefficient,  disk-swapping  of  data.  Pre-  and 
post-processing  modules  are  a  standard  part  of  the  ZeuS  package,  and  make 
the  tasks  of  defining  a  simulation  and  graphically  displaying  results 
straightforward.  The  post  processor  supports;  graphical  display  on  the 
standard  IBM  (and  compatible)  graphics  adapters:  the  OGA ,  EGA  and  VGA. 

The  normalized  speed  ol:  computation  varies  with  several  factors,  including 
the  type  of  PC  machine,  the  amount  of  contact  processing  required  by  the 
simulation,  and  the  size  of  Lhe  data  bus  anticipated  by  the  compiled  ZeuS 
code  (e.g.,  8,  16,  or  32  hit:  code).  At  the  low  end,  a  <t,77  MHz  IBM-PC/XT 
(8  bit  code)  runs  at  a  normalized  speed  of  processing  inherent  in  the 
computation.  On  the  other  hand,  ZeuS,  when  compiled  for  32  bit  80386 
machines,  and  running  at  20  MHz  processor  speeds,  can  achieve  normalized 
speeds  under  .003  se(./(cyc!e  note).  Such  computat  ional  speeds  for  codes 
of  this  variety  are  typically  found  only  on  mainframe  computers.  Addition¬ 
ally,  when  employing  trie  PC's  16  or  32  bit  protected  modes  of  operation, 

RAM  storage  of  state  variables  may  exceed  640  kbytes,  thereby  allowing 
simulations  to  be  run  which  use  much  more  than  2500  node’s;  and  4500  elements 

Results  of  calculations  are  presented  for  hy pet  velocity  impact,  perforation 
of  finite  plates  by  long  rod  projectiles  at  ordnance  velocities  and  Taylor 
anvil  impact  situations.  Good  agreement  is  shown  with  e> perimental  data, 
where  available. 


Stability  Conditions  for  Simple  Shearing 

A,  S,  Douglas 
H.  Tz.  Chen  and 
R,  Malek-Madani 

Since  few  physical  events  lead  to  simple  shearing  with  homogeneous 
fields,  it  is  important  to  investigate  the  stability  of  shearing  motions  which 
are  inhomogeneous.  Experiments  which  examine  shear  localization  using 
torsion  of  thin  cylindrical  specimens  usually  have  thick  sections  adjacent 
to  the  thin  test  section  in  order  to  grip  the  specimen.  These  thick  sections 
act  as  large  thermal  masses,  causing  significant  heat  transfer  from  the  test 
section  into  the  thicker  section.  This  cooling  of  the  edges  prevents  local¬ 
ization  at  the  interface  but  it  also  ensures  a  nonuniform  temperature  field 
and  (because  of  the  nature  of  the  material  constitutive  law)  velocity  field. 

The  problem  considered  here  is  that  of  the  shearing  of  a  thermo-visco¬ 
plastic  material  in  which  the  effects  of  heat  conduction  and  material  inertia 
are  included  explicitly.  Just  as  in  experiments,  the  inhomogeneous  fields 
arise  because  heat  is  conducted  from  the  boundaries  giving  rise  to  a  non- 
uniform  thermal  field,  which  dictates  a  non-uniform  velocity  distribution. 

A  perturbation  method  is  used  t,o  examine  the  stability  ot  the  steady 
solutions  to  the  governing  equations.  Since  the  resulting  set  of  partial  dif¬ 
ferential  equations  have  spatially  dependent  coefficients,  an  integral  method 
is  used  to  establish  criteria  for  the  existence  of  steady  solutions.  Three  dif¬ 
ferent  materials  laws  are  examined  and  conditions  for  stable  steady  motions 
are  obtained. 


MICROSTRUCTURE  TRANS  FORMAT  IONS  IN  HIRE  FOLYCRISTALUNK  a-TITANIUM  UNDER 
STATIC  AND  DYNAMIC  LOADING. 
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To  understand  t  hi'  mechanism  of  deformation  in  metals  of  hexagonal 
crystal  .structure,  compress  ion  tests  have  been  performed  at.  strain  rates 
in  ttie  range  from  10-4  to  10' J  s’  1  in  commercially  pure  titanium 
annealed  one  hour  at  700"C.  Two  deformation  modes  have  been  observed  to 
bo  active:  slip  und  twinning.  The  object  of  this  study  is  to  analyse  the 
effect,  of  strain  rate  and  specimen  orientation  on  slip  and  twinning 
occurence.  Transmission  electron  arid  optical  microscopy  observations 
have  I, een  performed  in  order  t.o  establish  a  relationship  between 
macroscopic  and  microscopic  phenomena. 

MECHANICAL  PROPERTIES 

The  st  ress-st  rain  curves  are  similar  under  static  and  dynamic 
conditions.  Nevertheless,  the  flow  stress  required  to  obtain  a  given 
deformation  is  lower  under  static  conditions  than  under  dynamic  ones, 

The  specimen  orientation  appears  to  have  an  influence  on  t.he 
mechanical  properties  under  dynamic  as  well  as  static  conditions. 
Compressions  along  ST  and  LT  directions  lead  to  the  same  flow  stress; 
flow  stress  is  lower  when  the  loading  axis  is  parallel  to  the  L 
d  i  rec t.  ion . 


M1CKOSTKUCTUKAL  CHARACTER 1 Z AT 1 ON 

The  microstruotural  study  reveals  that  twinning  is  an  important 
deformation  mode  even  at  room  temperature.  Fraction  of  twinned  grains 
increases  with  deformation  and  strain  rale,  The  number  of  twins  iri  a 
given  grain  increases  also  with  strain  rate. 

Transmission  electron  microscopy  allows  to  analyse  the  dislocation 
network.  For  both  static  and  dynamic  specimens,  dislocations  are 
arranged  in  bands.  The  mean  distance  between  two  walls  is  about  4  pm. 


C?<-J 


CONCLUSION 


The  flow  stress  increase  with  strain  rate  can  be  explained  by  a 
change  in  the  mechanism  which  controls  dislocations  motion;  twinning  can 
also  have  an  influence  but  its  participation  to  the  deformation  is  not 
yet  well  understood. 

Active  twinning  or  slip  modes  is  not  the  same  for  all  the  specimen 
orientations.  Easy  glide  are  predicted  to  be  predominantly  prismatic  In 
L  specimens,  basal  and  pyramidal  in  the  LT  and  ST  specimens.  This  can 
explain  the  different  mechanical  behaviour  of  the  L  specimens. 
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SEQUENCE  OF  THE  DEVELOPMENT  OF  DEFORMATION  STRUCTURE  OF  A1  AND  Cu 
SINGLE  CRYSTALS  UNDER  0.5  MBar  and  1  MBar  SHOCK  LOADING 

M. A. Mogilevsky,  L.S.Buahnev 
Lavrentyev  Institute  of  Hydrodynamics 
Siberian  Division  of  the  USSR  Academy  of  Sciences 
Novosibirsk  630090  USSR 

1.  The  plastic  deformation  under  shock  loading  may  take  place 
at  the  compression  front  on  exposure  in  the  state  of  peak  compressi¬ 
on  and  unloading.  Some  peculiarities  of  shock  deformation  were  reve¬ 
aled  earlier  |1|  by  examining  the  twinning  structure  (compression  and 
unloading: systems,  sequential  activation  of  different  systems,  a  full 
time  of  compression  at  the  front).  The  dislocation  structure  develop¬ 
ment  has  not  been  investigated  systematical ;y. 

2.  The  shock  loading  of  single  crystals  of  Cu  under  1  MBar  and 
A1  under  0.1  to  0.5  MBar  at  77*  with  a  1  to  2  %  residual  deformation 
has  been  realized. 

3.  As  far  as  the  dislocation  density  is  concerned,  the  well- 
known  results  on  oscillographic  shock  front  width  in  A1  under  40  GPa 
and  in  Cu  under  100  GPa  |2|,  being  about  X  to  3  ns,  allowed  estimating 
the  density  of  moving  dislocations  at  the  front.  It  is  in  a  signifi¬ 
cant  excess  of  a  full  density  of  dislocations  in  the  recovered  speci¬ 
mens,  that  is  indicative  of  the  pronounced  rearrangement  of  the  struc¬ 
ture. 

4.  The  polygonal  structure  formation  in  A.1  at  the  final  compres¬ 
sion  stage  and  during  exposure  is  attended  with  dislocation  reactions 
and  decrease  in  dislocation  density.  The  role  of  interaction  between 
dislocations  of  different  systems  is  also  confirmed  by  the  analysis 
of  the  dislocation  structure  of  Cu  under  5  GPa  and  by  a  computer  si¬ 
mulation  of  the  behaviour  of  the  crystalline  lattice  under  a  high 
one-dimensional  deformation  |3|. 

5.  In  the  shock-loaded  A1  and  Cu,  there  were  revealed  specific 
defects,  such  as  long  (  up  to  hundreds  of  microns)  strips  of  defor¬ 
mation  formed  during  rarefaction,  where  an  intense  deformation  was 
attended  with  disorientation).  A  lower  microhardness  of  Cu  in  the  strip 
and  lower  dislocation  density  make  it  possible  to  interpret  them  as 
adiabatic  shears,  which  were  not  previously  observed  in  copper. 

Ill  Mogilevsky  M. A.  ,  Physics  Reports , 1983 ,  v.97,  No.  6,  357, 

| 2 |  Chhabildas  L.C.,  Asay  J.R.,  J .Appl .Phyo. ,  1979,  v.50,,  p.2749. 

1 3 j  Mogilevsky  M.A.,In:  Proc.  Int.  Conf.  on  Impact  Loading  and  Dynamic 
Behaviour  of  Materials,  Bremen,  1987. 
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MECHANISMS  OF  DISLOCATION  MOTION  IN  T075-T73  ALUMINIUM  ALLOY 
AT  STRAIN  RATES  AROUND  10  V1. 


C  0  BURSTOW,  M  C  LOVELL  AND  A  L  RODGERS 


Tha  Royal  Military  Collage  of  Science  (Cranfield),  Shrivenham,  SWINDON, 
Wiltshire,  England. 


The  testing  of  aluminium  alloys  at  strain  rates  below  10*5* 1  has  led  to  the 
consensus  that  the  dominant  deformation  mechanism  is  thermal  activation  of 
dislocation  movement.  While  at  higher  strain  rates  there  is  some  indication  of 
viscous  drag  effects.  Our  results  on  7075-T73  do  not  fit  this  picture. 

Results  for  this  alloy  at  strain  rates  around  10 V1  were  reported  by  us  at 
Bremen  in  May  1987,  These  were  obtained  using  a  Direct  Impact  Hopkinson  Bar 
(DIHB).  The  DIHB  is  however  limited  at  high  strain  rates  by  its  time  resolution 
of  about  lps.  To  extend  our  results  to  higher  strain  rates  we  have  constructed  a 
Plate  Impact  Facility  (PIF).  Using  a  compressed  gas  gun  and  quartz  piezoelectric 
stress  transducers  a  time  resolution  of  12ns  has  been  achieved. 

To  obtain  as  much  information  as  possible  from  each  shot  the  following 
arrangement  is  used.  A  loading  plate  strikes  a  specimen  of.  the  aluminium  with 
a  quartz  gauge  mounted  on  its  rear  surface.  The  stress  wave,  induced  by  the 
impact,  is  measured  by  the  gauge  after  a  single  transit  across  the  specimen.  The 
separation  of  the  elastic  and  plastic  wavefronts  during  this  transit  allows  the 
Hugoniot  Elastic  Limited  (HEL)  to  be  observed. 

Measurements  of  the  HEL  for  7075-T73  have  been  made  over  the  strain  rate 
range  6.6  x  10V  1  to  4.9  x  10‘s*1  at  ambient  temperature.  The  results  correlate 
closely  with  our  earlier  DIHB  data  and  continue  to  show  no  significant  strain 
rate  dependence.  It  is  concluded  that  for  strain  rates  up  to  5  x  10‘s'1  and 
temperatures  of  20°C  or  above  the  limiting  factor  controlling  dislocation  motion 
is  not  thermal  activation  but  the  athermal  back  stress. 


EFFECTS  OF  GRAIN  SIZE  AND  STRAIN  RATE 
ON  THE  MECHANICAL  RESPONSE  OF  COPPER 


David  H.  Lassila 
University  of  California 
Lawrence  Livermore  National  Laboratory 
Livermore,  CA  94550 


ABSTRACT 

The  flow  stress  of  copper  is  known  to  increase  with  decrease  in  grain  size  and  this  behavior 
can  be  described  reasonably  well  by  the  Hall-Petch  expression; 

o(e)=»o6(e)4-k(e)d-lrt  m 


where  0(e)  is  the  tensile  flow  stress,  fl0(e)  is  a  constant  (sometimes  referred  to  as  frictional 
stress),  k(«)  is  the  Hall-Petch  constant,  and  d  is  the  average  grain  diameter.  The  Hall-Petch 
relationship  is  empirical,  however,  some  authors  nave  suggested  various  mechanisms  to 
explain  the  dependence  of  the  constants  on  material  chemistry  and  microstxucture  as  well  as 
loading  conditions.  Thus  far  all  work  in  this  area  has  indicated  that  the  constant  k(e) 
should  be,  as  a  first  order  approximation,  independent  of  deformation  rate. 

In  this  study  the  effects  of  grain  size  and  deformation  rate  on  the  flow  stress  of  OFE  copper 
were  examined.  Testing  was  performed  in  compression  at  ambient  conditions  at  strain 
rates  from  104  $•'  to  5  x  103  s->.  Testing  at  different  strain  rates  allowed  the  dependence  of 
the  Hall-Petch  parameters  on  strain  rate  to  be  examined.  The  effect  of  various 
crystolographic  textures  on  mechanical  behavior  was  also  examined.  The  findings  of  the 
study  are  as  follows: 

•  Decrease  in  the  grain  size  caused  an  increase  in  flow  stress.  The  observed 
dependence  of  flow  stress  on  grain  size  was  found  to  be  consistent  with  the  Hall-Petch 
relationship.  The  Hall-Petch  Parameter  k(<)  was  found  to  be  independent  of  strain  up  to  a 
strain  level  of  20%. 

•  Results  of  the  study  suggest  that  strain  rate  does  not  have  an  influence  on  the 
proportional  relationship  between  How  stress  and  grain  size,  i.e,,  the  Hall-Petch  parameter 
k(e)  is  not  strain  rate  dependent. 

•  ,  Crystallographic  texture  had  no  discemable  effect  on  ouer.s-stmin  response. 


Effects  of  Grain  Size  and  Temperature  on  the  Mechanical  Properties  of 
OFHC  Copper 


D  J  Parry  and  A  G  Walker 

Department  of  Physics,  Loughborough  University  of  Technology, 
Loughborough,  Leicestershire,  LE11  3TU.  UK 

Abstract)  A  comprehensive  investigation  has  been  undertaken  of  the 
dynamic  mechanical  properties  of  OFHC  copper  over  a  strain  rate  range 
of  about  2  x  10" ’a*'  to  6  x  lO’s*'  with  testing  temperatures  of  20CC 
to  600°C  and  grain  sites  from  20pm  to  240pm.  The  low  strain  rate 
region  was  studied  using  an  Instron  screw  machine,  while  at  the 
higher  rates  compressive  (to  6  x  1 0 3 s " ' )  and  tensile  (to  1.5  x 
lO’s*1)  split  Hopkinson  pressure  bar  systems  were  employed. 

The  results  can  be  divided  into  two  distinct  regions  according  to  the 
strain  rate.  For  rates  below  about  10 3 s " 1  the  flow  stress  increases 
slowly  with  strain  rate  in  a  manner  consistent,  with  a  thermally 
activated  mechanism.  The  strain-rate  sensitivity  X  ■  9o/91og£  (at 
constant  strain)  is  found  to  be  independent  of  strain  rate  and  grain 
size  but  increases  with  strain  and  temperature.  There  is  a 
progressive  reduction  in  the  flow  stress  with  increasing  grain  size 
and  temperature,  until  at  600°C  the  differences  in  flow  stress  for 
all  the  grain  sizes  becomes  insignificant. 

At  strain  rates  in  excess  of  approximately  U)3s''  the  flow  stress 
increases  much  more  rapidly  with  strain  rate  for  all  grain  sizes  and 
temperatures.  The  room  temperature  results  are  In  agreement  with  a 
linear  viscous-drag  type  of  relationship.  At.  elevated  temperatures 
there  are  insufficient  results  to  allow  any  difinito  conclusions  to 
be  made  regarding  the  form  of  the  strain  rate  dependence. 

The  validity  of  the  Hall-Petoh  relationship  has  also  been 
investigated  by  plotting  flow  stress  as  a  function  of  the  reciprocal 
of  the  square  root  of  the  grain  diameter  for  five  grains  sizes  over 
the  wide  range  of  strain  rate  and  temperature  indicated  above.  The 
linear  dependence  predicted  by  the  relationship  was  found  to  be 
valid  for  most  combinations  of  strain-rate  and  temperature.  However, 
at  20°C  and  the  low  strain  rate  of  2  x  10"3s*',  the  results  show 
evidence  of  a  bilinear  relationship  for  strain  above  1%,  with  a 
transition  grain  size  of  32pm. 
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J.  Duffy  and  S.  Suresh 

Division  of  Engineering,  Brown  University,  Providence,  Rl  02912,  U.S.A. 


The  potential  advantages  of  ceramic  matrix  composites  (CMC)  arc 
well-known;  CMC's,  for  instance,  are  the  only  material  of  reasonable  cost  able  to 
withstand  very  high  temperatures.  Their  great  disadvantage  at  present  is  their 
brittleness.  Tests  must  be  devised  to  evaluate  their  fracture  behavior  when  subjected 
to  impact  loading.  Our  results  indicate  that  this  is  feasible  if  the  specimen  is  first 
properly  fatigue-precracked  and  carefully  loaded. 

Experiments  are  described  in  which  stress-wave  loading  is  employed  to 

fracture  a  notched  CMC  specimen  with  a  prefatigued  crack  at  the  root  of  the  notch. 

Tests  are  performed  both  in  Mode  I  and  Mode  III,  with  specimens  of  a  fine 
polycrystalline  AJj03  and  with  an  AljOj  containing  25v/o  SiC  whiskers.  Quasi-static 
tests  are  also  performed  and  we  arc  thus  able  to  measure  the  fracture  toughness 
values  both  statically  and  under  truly  dynamic  conditions  and  for  two  fracture 
modes.  The  critical  advance  comes  in  specimen  preparation:  our  development  of 
techniques  making  possible  the  growth  of  fatigue  cracks  in  CMC’s  by  cyclic 
compressive  loading,  in  the  dynamic  fracture  test,  the  stress  wave  is  initiated 
explosively.  The  geometry  and  loading  insure  plane  strain  conditions  at  the  crack  tip. 
A  full  finite  element  analysis  is  available  confirming  the  accuracy  of  results.  Values 
of  can  be  calculated,  in  addition,  the  experimental  method  has  been  adapted  to 
dynamic  Mode  III  fracture.  The  present  experimental  results  for  Ai.O,  reveal  that 
the  fracture  initiation  toughness  values  are  about  50%  higher  under  dynamic  loading 
in  Mode  I  than  quasi-statically,  The  corresponding  ratio  is  30%  for  Mode  III.  In 
Af2Oj  •  SiCw,  these  ratios  decrease  to  about  30%  in  Mode  t  and  only  about  5%  in 

Mode  III.  Scanning  electron  microscopy  reveals  that  Mode  I  fracture  initiation  occurs 

predominantly  by  an  intergranular  failure  mode  irrespective  of  whether  the  loading  is 
quasi-static  or  dynamic.  Dynamically,  however,  there  u  evidence  of  a  relatively 
greater  percentage  of  transgranular  failure.  In  distinct  contrast  with  the 
macroscopically  flat  fracture  surfaces  observed  under  tensile  loading,  a  liigh!< 
tortuous  and  rough  fracture  surface  was  observed  for  Mode  III  fracture  in  both  the 
quasi-static  and  dynamic  tests. 

Microcracks  running  dynamically  in  a  brittle  solid  influence  both  the  clastic 
moduli  and  strength,  (f  such  a  material  contains  a  population  of  microcracks,  the 
velocity  of  the  microcracks  and  the  conditions  for  microcrack  coalescence  in  front  of 
the  macro-crack  tip  will  determine  the  measured  values  of  tensile  fracture  toughness. 
A  model  is  proposed  which  determines  an  estimate  of  the  stiffness  loss  caused  by  a 
dilute  population  of  microcracks  aligned  in  a  direction  normal  to  the  far-field  tensile 
stress  axis.  This  solution  is  then  used  to  calculate  the  stiffness  loss  in  brittle 
polycrystalline  solids.  The  dynamic  to  static  fracture  toughness  ratio  is  predicted  as  a 
function  of  microcrack  velocity  and  elastic  properties.  Experimentally  measured  values 
of  dynamic  to  static  fracture  roughness  ratio  in  A 1 203  and.  Alj03  with  SiCw  fall 
within  the  bounds  predicted  by  the  analysis. 
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Structure  *  Property  Characterization  of  a  Shock-Loaded  Boron  Carbide- 
Aluminum  Cermet 


W.R.  Blumenthal  and  G.T.  Gray  lit 

Materials  Science  And  Technology  Division,  Los  Alamos  National  Laboratory, 
Los  Alamos,  New  Mexico  87S45  USA 


ABSTRACT 

The  propagation  of  shock  waves  through  metals  and  alloys  is  known  to 
induce  metallurgical  changes  such  as  the  formation  of  dislocations, 
deformation  twins,  point  defects,  and  phase  transformation  products  in  the 
microstructure.  Consequently,  the  mechanical  properties  of  metals  and 
alloys  (i.e.  hardness  and  yield  strength)  are  highly  dependent  on  their 
shock  histories,  in  particular,  the  shock  pressure  and  shock  duration. 
Under  quasi-static  compression  loading  most  ceramics  and  cermets  initially 
respond  elastically  and  subsequently  fail  at  a  critical  stress  with  little 
or  no  macroscopic  yielding.  The  mechanical  response  is  controlled  by  the 
elastic  properties  of  the  material  followed  by  the  intervention  of  one  or 
more  microstructural  flaws.  Microstructura'I  studies  on  shock- recovered 
minerals  (quartz,  anorthite,  and  periclase)  have  revealed  planar  features, 
"shock  lamallae",  implying  inhomogeneous  local  plastic  flow  above  the 
Hugoniot  elastic  limit  ( HEL )  as  suggested  by  Grady  (1977).  Also  grain 
boundary  flaws  and  microcracks  can  lead  to  fragmentation  of  poly¬ 
crystalline  ceramics  at  compressive  shock  pressures  below  the  HEL.  Hence, 
the  development  of  high  fracture  toughness  cermets  is  of  interest  since 
fragmentation  is  a  serious  limitation  to  impact  performance. 

Shock-recovery  experiments  on  a  65  vol%  boron  carbide-35  vo1%  aluminum 
cermet  supplied  by  the  University  of  Washington  were  performed  as  a 
function  of  shock  pressure.  The  purpose  was  to:  a)  develop  shock-recovery 
techniques  for  cermets  and  ceramics  which  reduce  the  destruction  of  the 
microstructure  by  accompanying  tensile  stress  waves  enabling  intact 
recovery  of  the  shock-loaded  samples  ard  b)  investigate  changes  in  the 
microstructure  and  mechanical  properties  of  a  high  fracture  toughness 
aluminum-boron  carbide  cermet  as  a  function  of  peak  shock  pressure.  The 
investigation  also  compares  quasi-static  and  dynamic  (Hopkinson  Bar) 
uniaxial  compression  results.  The  objective  in  this  case  was  to  determine 
the  effect  of  moderate  strain  rate  loading  on  the  mechanical  properties 
(modulus  and  compressive  strength)  and  failure  behavior  of  the  cermet. 

The  material  was  made  by  infiltrating  pure  liquid  aluminum  into  a  pre¬ 
sintered  porous  boron  carbide  skeleton  resulting  in  a  finely  inter¬ 
connected  microstructure  with  an  average  B4C  phase  size  of  7  microns  and 
a  porosity  level  below  1.5%.  TEM  characterization  of  the  as-received 
material  revealed  the  substructure  of  the  cermet  to  consist  of  pre¬ 
dominantly  defect-free  B4C  grains  interspersed  with  the  aluminum  regions. 


Shock- loading  experiments  were  performed  on  the  cermet  as  a  function  of 
shock  pressure.  Samples  were  recovered  largely  Intact  due  to  the  use  of 
confinement  and  "soft-recovery"  techniques  (Gray  1988).  Post-shock 
examination  by  optical  and  transmission  electron  microscopy  (TEM)  showed 
the  B*C  structure  to  be  unchanged  up  to  5  GPa,  while  the  A1  phase 
exhibited  a  high  density  of  randomly  distributed  dislocations.  After  a 
shock  of  10.6  GPa  some  of  the  B^C  grains  displayed  dislocation  debris,  but 
the  majority  showed  no  evidence  of  plastic  deformation.  Uniaxial 
compression  testing  of  the  unshocked  cermet,  was  performed  at  two  nominal 
strain  rates,  10‘3  s'1  and  103  s'1  iF’9ure  !)■  Dynamically,  the  material 
showed  limited  yielding  after  about;  1.5%  elastic  strain  and  subsequently 
accumulated  about  1%  Inelastic  strain  (2.5%  total  strain)  prior  to 
failure.  Quasi -statically,  the  cermet  began  yielding  at  about  1%  strain 
and  only  accumulated  between  0.5%-0.9%  inelastic  strain  (1.5%-1.9%  total 
strain)  prior  to  failure.  The  unshocked  cermet  showed  a  9% - 24%  increase  in 
strength  when  the  strain  rate  was  Increased  from  10'3  s'1  (1.94  GPa)  to 
103  s'1  (from  2.1  GPa  to  2,4  GPa)  . 

!n  summary:  1)  momentum  trapping,  confinement,  and  "soft"  arrest 
techniques  are  critical  to  successful  shock  recovery  and  subsequent 
mlcrostructural  analysis.  2)  The  predominant  response  of  the  cermet  to 
both  the  5  and  10.6  GPa  shocks  Is  elastic,  although  the  aluminum  phase 
behaved  plastically.  3)  No  localized  damage  was  observed  at  shock 
pressures  of  either  5  GPa  or  10.6  GPa.  4)  Inelastic  deformation  and 
strength  improve  with  strain  rate  under  uniaxial  compressive  loading. 
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Figure  1:  Stress-strain  plot  for  both  quasi-static  and  Hopkinson  bar 
compression  loading  of  a  65%  B4C  -  A1  cermet. 


Impact  Strength  of  Ceramics  at  High  Temperatures 
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Introduction 

This  report  aims  to  present  an  outlook  of  mechanical  properties  of  some  engi¬ 
neering  ceramics  under  slow  impact  load  and  their  testing  methods.  Three 
kinds  of  popular  ceramics  (alumina,  SiC  and  TTZ)  are  chosen  for  characteri¬ 
zation^  impact  properties  up  to  1200"C.  Conventional  low  rate  tests  were 
also  performed  to  assess  the  rate  effect.  The  present  work  consists  of  (1) 
the  strengtn-data  are  characterized  at  RT  through  three  point  bend,  Brazilian 
and  tensile  tests.  The  bending  strength  %  i.s  assessed  and  compared  with 
tensile  data  by  the  use  of  Weibuli’s  failure  probability  theory.  Brazilian 
test  data  are  also  compared  with  these  data  and  (2)  high  temperature  strength 
of  these  ceramics  is  characterized  by  three  point  bend' tests.  Character¬ 
istics  of  fracture  pattern  arc  also  reported  in  conjunction  with  deformation 
rate  and  temperature. 

Specimen 

Tested  ceramics  are  99. S",  purity  alumina,  95T>purity  SiC  and  TTZ  (transforma¬ 
tion  toughened  ZrQ.).  Bend  specimens  are  5x5mm  square  cross-sectional  bar 
(span  length  is  25 . 4-'-25 . 8 ) .  Brazilian  specimens  are  0^x5* .  Tensile 

specimens  of  alumina  are  3.2 'l'x3Cmm. 

Experimental  Method  and  Data 

Three  kinds  of  Hopkinson  bar  type  impact  tests  were  carried  out;  three  point 
bend  test  (Figure  1),  Brazilian  test  (Figure  2)  and  tensile  test  (Figure  3). 
The  strength-data  through  these  tests  are  shown  in  Table  1  and  Figure  4. 

The  ratio  between  three  point  bend  and  tensile  strength  at  RT  ,  is  1.6 

for  Impact  test  and  2.0  in  low  rate  test.  On  the  other  hand,  the  strength 
ratio  throuqh  Wei  bull  failure  probability  theory  is  1.6IK1.75  for  the  modu¬ 
lus  ttp'KIO  (the  »</-value  is  from  bending  data  in  low  rate  tests  at  RT.) 

( o j ;  the  tensile  strength  in  tensile  test,  (oj/.;  the  tensile  strength  by 
Brazilian  test. ) 


Figure  1  Experimental  set-up  of  bend  test 
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Figure  2 
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ABSTRACT 


The  dynamic  properties  of  monolithic  ceramics  are  usually  altered 
by  the  presence  of  a  second  phase.  Impedance  mismatch  is  the  primary 
factor  responsible  for  microstructural  damage  in  the  region  of  the 
phase  boundaries.  Planar  shock  recovery  experiments  have  shown  that 
microcracks  in  the  alutaina/glass  interface  are  generaly  initiated  at  a 
stress  level  below  the  obsereved  Hugoniot  Elastic  Limit  (HEL). 

It  turns  that  the  spall  strength  decreases  with  increasing  of  the 
impact  pressure.  These  studies  also  showed  that  the  onset  of  the  loss 
of  spall  strength  is  almost  the  same  as  the  microcracks  initiation 
stress  level.  This  indicates  that  spall  strength  is  associated  with 
microstructural  damage  in  two  phase  ceramics. 

In  this  paper  we  also  calculate  che  threshold  shear  scress  at 
which  microcracks  are  initiated  using  simple  arguments  relating  to  the 
different  shock  properties  of  che  two  phases.  Thus  we  assume  that  the 
mismatch  in  particle  velocity  between  che  two  phases,  is  the  cause  to 
Che  Interfacial  shear  stresses  at  the  phase  boundaries.  From  our 
calculation  we  conclude  that  Che  shear  threshold  stress  for 
microcracks  formation  depends  on  grain  size,  the  Hugoniot  mismatch  and 
the  impact  pressure.  These  calculations  are  also  in  good  agreement 
with  the  average  shear  stresses  in  the  shocked  specimen  as  determined 
with  longitudinal  and  Transvers  Manganin  gauges. 
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In  recent  years,  there  has  been  a  growing  interest  in  the  response  c>f  ceramics  to  ballistic 
impact.  This  paper  describes  a  study  of  the  impact  behaviour  of  a  range  of  ceramics  when 
impacted  by  steel  spheres.  High  speed  photography  with  a  Hadland  790  MACON  camera 
was  used  to  observe  the  modes  of  failure  of  the  ceramics,  damage  to  the  projectile  and  to 
obtain  quantitative  data  on  the  kinetics  of  impact.  The  projectiles  were  accelerated  by  a 
single  stage  gas  gun  of  the  double-diaphragm  type.  The  velocity  range  covered  was  30  to 
1000m  s*L 

A  large  number  of  ceramic  types  and  projectile  materials  have  been  investigated.  It  is  only 
possible  to  present  a  few  examples  in  the  paper  and  Table  1  gives  information  on  their 
relevant  physical  properties.  The  projectiles  were  5mm  diameter  hardened  steel  spheres  in 
all  cases  except  where  stated.  The  ceramics  were  in  the  form  of  50mm  x  50mm  tile 
specimens  with  thicknesses  in  the  range  3  to  20mm. 

Table  1 .  MateriaLpropertics  of  target  and  projectile  materials 


Material 

Density 

Hardness 

Hv/GPA 

Fracture 

Shear 

Bulk 

Young 

Poissons 

p/103kg  in'3 

Toughness 
Klc/MN  m-3/2 

»'aUS 

Ws 

Ratio 

\) 

Soda-lime 

2.50 

4.5  ±  0.4 

0.75  ±  0.03 

27.6 

42.5 

79.2 

0.19 

glass 

±  0.8 

±  2.7 

±  3.0 

±0.02 

Glass 

2.47 

9.2  ±0.5 

1.77  ±0.12 

34.8 

45.5 

83.3 

0.20 

ceramic  A 

±  0.9 

±  3.4 

±  4.0 

±  0.03 

Glass 

2.59 

8.2  ±  1.1 

0.42  ±  0.05 

34.3 

59.0 

92.0 

0.24 

ceramic  B 

±0.9 

±2.6 

±4.0 

±0.01 

Alumina 

3.69 

11. 9±  1.6 

2.95  ±0.17 

108 

147 

260 

0.21 

±4 

±  10 

±  10 

±0.03 

Hardened 

7.80 

10.0±  0.5 

77.4 

163.7 

201 

0,30 

steel 

(a)  Impacts m, alumina 

In  this  case,  the  hardness  of  the  target  is  greater  than  that  of  tire  projectile  (see  Table  1).  For 
impact  velocities  up  to  ~  1 80m  the  sphere  rebounds  from  the  target.  The  high  pressures 
generated,  due  to  the  high  hardness  of  tl.e  alumina,  result  in  extensive  plastic  deformation 
of  the  projectile  and  the  formation  of  a  Hat.  If  the  impact  velocity  is  increased,  fracture  of 
die  projectile  occurs. 

The  impact  of  a  hardened  steel  sphere  on  an  alumina  tile  specimen  produces  a  Hertiuan-type 
cone  failure  over  the  velocity  range  investigated  (50  to  900m  s~r).  For  impact  velocities 
greater  than  ~220m  s-1,  the  cone  crack  propagates  completely  through  an  8.6mm  thick 
specimen  resulting  in  a  conical  fragment.  The  semi-apex  angle  of  the  cone.  0,  is  less  than 
that  obtained  by  static  loading.  As  the  velocity  of  impact  increases,  0  increases  reaching  a 
value  of  -57®  at  the  highest  velocities.  The  cone  has  a  smooth  upper  region  and  a  lower 
roughened  region,  The  roughened  reginn  is  caused  by  the  interaction  of  the  propagating 


cone  crack  with  the  reflected  tensile  wave  from  the  rear  face  of  the  specimen.  The  change 
of  angle  for  the  lower  part  of  the  cone  can  also  be  explained  by  considering  the  stress  wave 
interaction. 

(b)  Impacts  on  glass.£sramiatypc,A 

In  this  case,  the  hardness  of  the  target  is  less  than  that  of  the  projectile  (sec  Table  1),  and 
very  little  damage  of  the  sphere  results  even  at  the  highest  velocities,  Jmpuets  of  -40m  s“* 
generated  Hertzian-type  ring  cracks  in  this  material.  Above  55m  s_1,  the  cone  cracks  reach 
the  rear  of  10mm  thick  specimens  producing  a  conical  fragment.  As  the  velocity  of  impact 
increases  the  semi-apex  angle,  0,  decreases.  Lateral  cracks  produced  from  the  deformed 
zone  at  the  impact  site  form  at  impact  velocities  greater  than  80m  s~s  Spall  failure  was 
evident  in  this  glass  ceramic,  for  10mm  thick  specimens,  at  velocities  above  ~120m  s-1. 
At  impact  velocities  >300m  s“*.  a  second  cone  crack  was  formed  in  the  material  with  a 
much  larger  semi-apex  angle  than  the  initial  cone  crack.  The  radial  cracks  which  are 
formed  by  bending  ate  evident  at  velocities  above  -60m  s~  * .  The  number  of  radial  cracks 
increases  with  impact  velocity.  Penetration  of  10mm  thick  specimens  takes  place  at 
velocities  above  -450m  s“*. 

(c)  Inipact5„on.5Qdarlimc  glass 

The  damage  modes  with  soda-lime  glass  are  similar  to  those  in  the  glass  ceramic  though 
they  take  place  at  lower  velocities.  For  example,  the  Hertzian  cone  crack  ejects  at 
velocities  above  -30m  s-1.  The  hardness  of  soda-lime  glass  is  much  less  than  that  of 
hardened  steel  (see  Table  1)  and  negligible  damage  is  suffered  by  the  projectile.  At  high 
velocities,  failure  modes  such  as  the  Hertzian  cone,  lateral  and  median  cracking  and  (with 
thin  plates)  radial  fracture  which  can  all  be  obtained  by  quasi-static  loading  are  increasingly 
dominated  by  stress  wave  controlled  fracture. 

(d)  Discussion 

The  relative  hardnesses  of  projectile  and  specimen  are  of  great  importance.  In  the  one 
category,  are  soda-lime  glass  and  the  glass  ceramics  with  lower  hardnesses  than  that  of  the 
piojectile  and  which  do  not  deform  the  projectile.  In  the  second  category,  are  aluminas, 
boron  carbide  etc...  with  higher  hardnesses.  They  deform,  break  up  and  eventually  shatter 
the  hardened  steel  projectile.  This  different  behaviour  has  a  great  affect  on  the  energetics  of 
impact.  The  toughnesses  of  the  materials  have  a  relatively  small  effect  on  the  behaviour. 

Impacts  onto  ceramics  which  arc  harder  than  the  projectile  show  that  at  300m  s'1  about 
80%  of  the  impact  energy  is  consumed  in  plastic  deformation  of  the  projectile  and  about 
15%  is  taken  by  the  rebounding  fragments.  Only  a  small  amount  of  energy  is  needed  to 
produce  the  fractures  in  the  specimen  and  the  conical  fragment  which  ejects  at  only  about 
4m  s'1'  Impacts  with  brittle  specimens  which  arc  softer  than  the  projectile  do  not  deform  it. 
Again  little  energy  is  consumed  in  fracturing  the  specimen.  A  large  amount  of  energy  is  left 
in  the  projectile  and  the  fragments  which  eject  from  the  rear  of  the  ceramic. 

High  speed  photography  has  proved  of  great  value  in  identifying  the  different  modes  of 
failure  of  both  the  specimen  and  the  projectile  during  the  ballistic  impact  of  ceramics.  A 
camera  such  as  the  IMACON  with  its  ease  of  synchronisation,  ability  to  record  at  modest 
light  levels  and  production  onto  polaroid  film  has  allowed  each  experiment  to  be 
photographed  routinely.  Combined  with  accurate  impact  velocity  measurement  and  the 
results  from  the  momentum  balance  we  now  have  a  great  deal  of  data  on  a  wide  range  of 
ceramics. 
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ABSTRACT!  ■■  T1"* ,, 

.Two  new  experimental  methods  were  used  to  Investigate  model 
micronofnposites  for  quasi  a  tat.  ip  and  dynamic  fracture  behaviour. 

fracture:’ modes  at,  low  and  very  nigh  strain  rule  Impact  experiments 
were  studied  using  model  unidirectional  composites.  The  model 
miorpromposites  consist  of  a  set  of  single  fibers  embedded  in  an  epoxy  film 
such  thd't  the  inter  fiber  distance,  number  of  fibers,  interface  and  resin 
properties  can  b«>  control  led.  Micronomposi te  monolayers  can  demonstrate 
basic  failure  modes  of  reinforced  composites. 

The  preparation  of  controlled  geometry  model  composites  was  possible 
due  to  a  new  method  developed  by  Wagner  et  a,i.  [1,2,3]. 

The  quaslstatic  d'e format  ion  and  fracture  modes  were  tested  using  a 
specially  designed  custom-made  miniature  tensile  testing  machine  [2'J.  The 
sequences  of  the  fracture  events  were  recorded  using  a  colour  video  camera 
fitted  to  a.  polar l sea  light  stereo  zoom  microscope.  The  results  of  tensile 
testing  are  Summarised  in  table  no.  1 

The  dynamic,  tools  at  ultra  nigh  strain  rat'}  were  performed  using  short 
pulsed  laser  induced  shock  waves.  This  new  method  developed  by  Gilath  ct. 
al.  [II-8J  for  metal  folia,  was  appl  led  to  study  the  fracture  mode  of  thin 
tapes  of  model  mlcrooompoaltcs.  To  avoid  laser  ablation  of  the  thin 
samples,  the  shock  wave  was  delivered  through  a  thin  aluminum  foil.  Due  to 
the  small  diameter  of  the  laser  beam,  high  shock  wave  pressure  and  very  fast 
energy  deposition,  short  pulsed  laser  experiments  are  the  most  suitable  to 
study  dynamic  fracture  for  small  samples  such  as  microcomposit.es. 

The  second  advantage  of  laser  induced  shock  wive  experiments  is  the 
possibility  to  obtain  very  close  energy  density  steps  at  the  amplifier 
stages  of  the  laser  system  thereby  enabling  the  intermittent  monitoring  of 
the  various  failure  modes  in  the  sample. 

The  following  values  wore  measured,  set.  table  no. .2: 

a)  The  energy  density  necessary  to  lift  the  epoxy  film  from  the  aluminum 
support. 

b)  Tne  energy  density  necessary  lo  induce  damage  in  epoxy  film  when 
reinforcement  is  undamaged. 

o)  The  energy  density  for  complete  rupture  of  epoxy  and  fibers  ir. 
microoompositfi. 


&<=~l 


Material 

Young  Modulus 
E(GPa) 

Failure  strain 

% 

Strength 
1.  MPa] 

Kevlar  29 

58.9 

'i.O 

i’b'io 

CYa&3/HY956 

2.0 

3.0 

31 

8  fiber  monolayer 

2.8 

3.2 

65 

Table  1.  Mechanical  properties  of  single  fiber,  matrix  and  microcomposlte. 


Experiment 

Energy  density 
delivered,  J/orrr 

Pressure  in 
aluminum  at 
interphase, 

Khar  (simulation) 

Pressure  transmitted 
to  epoxy,  kbar 

A 

75110 

9 

1 

J 

Cr.  i 

B  ; 

206+10 

15 

5.7 

C 

300+10 

20 

7.6 

Table  no,  2,  Energy  densities  and  pressures  obtained  in  pulsed  laser 
experiments. 
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As  part  of  a  wider  programme  to  study  the  impact  response  of  hybrid  woven 
reinforced  carbon/glass  laminates  it  has  been  necessary  to  make  an  estim¬ 
ate  of  the  effect  of  strain  rate  on  the  failure  strengths  of  the  individual 
reinforcing  plies  in  tension  and  compression  and  of  the  interlaminar  shear 
strength  on  the  planes  between  neighbouring  plies.  This  requires  tests  to 
be  performed  at  a  quasi-static  and  an  impact  rate  on  non-hybrid  specimens 
reinforced  with  woven  layers  of  either  carbon  or  glass.  The  present 
paper  describes  the  experimental  techniques  and  specimen  designs  used  at 
both  rates  of  loading  in  tests  on  the  carbon  reinforced  specimens. 

For  loading  in  both  tension  and  compression  the  same  design  of  specimen,  a 
thin  strip  walsted  in  the  thickness  direction,  see  Fig.l,  was  used  at  both 
rates  of  strain.  The  quasi-static  tests  were  per- 
formed  on  a  standard  screw-driven  Instron  testing 
/  machine  while  the  impact  tests  used  a  modified  split 

!  Hopkinson's  pressure  bar  apparatus.  Typical  stress¬ 
's  strain  curves  obtained  are  shown  in  Fig. 2.  A 

|  marked  increase  in  the  failure  strength  and  some 

_  increase  in  both  elastic  modulus  and  failure  strain 
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b)  Compression 


Fig. 2  Effect  of  Strain  Rate  on  Stress-Strain  Curves 
in  tension  and  compression 


Fig.l  Tension/Compression  Specimen 
(all  dimensions  in  mm) 


THE  IMPACT  ON  A  COMPOSITE  PLATE  : 

A  BOTH  THEORETICAL  AND  EXPERIMENTAL  APPROACH 
OF  THE  IN-PLANE  AND  TRANSVERSE  EFFECTS. 
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INTRODUCTION 

The  impact  ol  a  sphere  on  a  carbon  liber  reinforced  plastic  plate  (CFRP)  may  be  considered 
as  an  Impactor's  kinetic  energy  transfer  to  the  plate.  This  energy  separates  into  various 
types  of  energy,  the  relative  importance  ol  which  depends  on  various  parameters  such  as 
the  Initial  velocity  ol  the  indenter  and  the  energy  release  rate  ot  the  material. 

The  evaluation  ol  the  relevance  ot  each  type  ol  energy  has  been  previously  performed  by 
several  in-situ  and  post-mortem  observations.  The  fact  that  post-mortem  analyses  only 
shows  final  Mates  o(  the  Impacied  plates  and  not  the  chronological  evolution  duo  to  the 
various  mecnanisms  has  already  been  stressed  [Beaumont  N.  1988a]. 

EXPERIMENTS 

The  impact  experiments  where  performed  using  a  gas  gun.  The  spherical  indenter  (d  :  6.36 
mm)  is  made  ol  steol.  A  high  speed  camera  and  strain  gages  whore  used  to  perform  in  situ 
observations.  After  the  impact,  several  post-mortem  observations  were  pertormod  using 
US-SCAN  and  optic  microscopy.  All  thoso  methods  have  been  deveioppod  in  another  paper 


MOOELISATION 

A  complete  model  has  been  established  in  the  frame  ot  non-equilibrium  thermodynamics 
The  sphere-plate  system  is  assumed  to  be  fully  described  by  a  number  ot  external  and 
internal  variables  and  by  two  thermodynamical  potentials.  (Rq.  l) 

Experiments  results  have  led  to  impose  the  tollowing  model-nypotheses  : 

The  indenter  penetrates  into  the  plates.  1  fiat  is  to  say  that  the  Hertzian  contact 
energy  can  be  neglected  [Shivakumar  t9U5| 

-  The  plate  has  an  elastic-brittle  behaviour.  That  Is  to  say  that  the  plasticity  ol  Ihe 
matrix  is  neglected  and  the  energy  release  rate  remains  constant  during  the  whole 
dalamination. 

-  The  delamination  ot  the  plate  is  tully  governed  by  the  kinematic  of  the  impact  and 
spreads  at  Ihe  indentation  spreading  velocity  noted  U 

-  The  preferential  delamination  direction  is  correlated  with  the  flexion  moduli  ot  the 

plate. 


-  The  inertia  etleets  are  neglected 

CONCLUSION 

Several  in-situ  observation  systems  have  been  deveioppod  in  order  to  have  more 
information  on  the  behaviour  ot  a  CFRP  plate  impacted  by  a  spherical  steel  indenter,  The 
physical  understanding  ol  the  phenomenon  lead  to  the  construction  of  a  model  established  in 
the  frame  of  non-equilibrium  thermodynamics 

Thanks  to  this  model,  it  is  possible  to  predict  the  residual  speed  ol  a  spherical  indenter 
impacting  a  CFRP  plate  with  a  quasi -isotropic  lay-up  in  the  case  of  an  elastic  brittle- 
behaviour  ol  the  plate.  In  the  case  ot  strongly  anisotropic  layers,  another  model  has  to  be 
considered  taking  into  account  the  lay  up  effects. 

For  any  lay-up  and  any  material,  as  long  as  delamination  takes  place  during  impact  and  has 
not  reached  saturation,  ii  is  possible  to  prodict  tiro  shape  and  surface  ot  the  delaminated  zone 
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AR5IBACI 

We  wish  to  discuss  the  phenomena  of  low  velocity  elastic  impact  In  AS4/3501-6  graphlte/epoxy 
composite  plates.  To  accomplish  this  objective  an  Investigation  combining  computational, 
experimental  and  precise  manufacturing  methods  was  developed.  Computationally,  a  finite 
element  -  mixture  theory  formulation  is  employed  to  capture  mlcrostructural  detail  at  the 
composite  Interface.  Experimentally,  dynamic  moire  Interferometry  is  developed  and  employed 
to  study  transient  mlcrostructural  displacement  fields  on  the  composite  free  edge  under  Impact. 

Four  types  of  composite  plates  were  manufactured  using  the  AS4/3501  -6  fibers.  The  specimens 
consisted  of  alternating  layers  of  0  and  90  degree  fiber  orientations  with  constant  plate 
dimensions  equaling  50  mm  wide,  7  mm  thick  and  250  mm  In  length.  The  layup  sequence  was 
varied  In  tho  plates.  Specifically,  the  first  series  of  plates  consisted  of  pure  90  degree  fibers; 
the  second  sequence  consisted  of  [Oie/SOie/Oieli the  thlrd  sequence  was  [02/902)3  for  24  plies 
and  the  final  sequence  was  [0/90]8  for  48  plies. 

The  Impact  process  occurred  by  a  striker  bar  impacting  the  composite  plate  parallel  to  the 
layer  direction  at  11  m/s.  Moire  interferometry  displacement  fields  of  the  pulse  propagation 
on  the  composite  free  edge  wore  obtained  using  a  20  n3.  pulsed  ruby  laser.  These 
displacement  fields  are  obtained  both  near  and  far  field  from  impact.  Figure  1  details  the 
longitudinal  displacement  field  far  field  from  Impact.  The  bulk  fringes  are  the  result  of  the 
bulk  lamb  wave  propagation  The  reticulation  or  microstructure  within  the  bulk  fringes 
correspond  to  the  layering  of  the  composite.  Figure  2,  shows  that  this  microstructure  Is  the 
result  of  the  differing  wave  velocities  In  the  0  and  90  degree  fiber  orientations,  8  and  3  km/s 
respectively  Under  stronger  elastic  Impacts  the  interface  shear  of  figure  1  would  grow 
developing  into  a  delamination,  Simitar  displacement  fields  exhibiting  microstructure  are  found 
on  the  plate  free  edge  for  the  24  and  3  ply  specimens.  Uniform  displacement  fields  are  found 
on  the  pure  90  degree  specimens.  Finally  displacement  fields  for  a  natural  delamlnatlon 
occurring  at  the  Interface  In  the  3  ply  specimen  will  be  presented.  This  lecture  will  compare 
the  non-delaminated  displacement  fields  to  the  finite  element  •  mixture  theory  formulation  to 
asses  a  broader  understanding  of  the  complete  parameters  causing  Interface  delamlnatlon  under 
low  velocity  Impact. 


Dynamic  Moire  Interferometry 

G/E  Composite  Stress  Wave  Behavior 


Quasi-Longitudinal  Phase  Velocity  Data  Versus 
Angle  of  Propagation,  <t>  in  the  Xi  -  X2  Plane  for 
a  Unidirectional  Composite  Laminate* 
(AS4/3501-6)  (100  x  304  x  8.3  mm) 


*L.  Pearson,  B.  Murri, 
D.  Gardiner 
Hercules  Aerospace 
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Figure  2,  the  phase  velocity  verses  angle  of  propagation  relationship  for  AS4/3501-6 
graphite  epoxy  prepreg  (courtesy  8.  Murri  and  D.  Gardiner  Hercules  Aerospace). 


Dynamic  indentation  of  two  dimensional  ring  systems 


T  Y  Roddy,  A  3  Kaddour  and  S  R  Reid 
Department  of  Mechanical  Engineering,  UMIST, 
PO  Box  88,  Manchester,  England 


ABSTRACT 

Dynamic  crushing  of  single  rings  (Reid  and  Reddy  (1879)),  ring  chains  (Reid  and 
Reddy  (1983)  and  Reid  and  Ball  (1984))  as  well  aa  that  of  two  dimensional  ring 
ayatams  (Stronge  and  8hia  (1987,1988))  have  been  studied.  This  paper  ia 
concerned  with  dynamic  indentation  of  hexagonally  close  packed  systems  made  of 
rings  (19.05  mm  diameter  and  12.7  am  length)  under  impact  conditions,  to  extend 
the  understanding  of  the  energy  absorbing  characteristics  of  such  systems  as  well 
as  those  of  cellular  structures  and  solids.  Typical  systems  were  288  ms  (15 
rings)  wide  and  267  ms  (16  rows)  deep  and  were  enclosed  in  a  rigid  frame  with  a 
central  opening  of  96.5  ms  wide  on  one  side. 

A  rigid  mass  of  0.16  kg  (95.25  am  wide  and  12.7  an  height)  was  propelled  along  a 
guided  track  by  a  cartridge  gun  to  impact  the  ring  systems  centrally  at; 
velocities  in  the  range  of  36  as-1  to  72  ms-1.  Deformation  process  was  recorded 
by  using  a  high  speed  camera  with  a  film  speed  of  about  6000  frames  per  second. 

TWo  different  deformation  patterns  depending  on  the  ratio  of  the  energy  absorbing 
capacity  of  the  struck  rings  to  the  kinetic  energy  of  the  striker  defined  as  the 
Q  factor  by  Stronge  and  Shim  (1987)  were  observed.  At  low  values  of  Q  the 
deformation  started  with  an  initial  crushing  mode,  in  which  the  rings  were 

directly  crushed  (with  little  lateral  expansion)  as  if  they  were  constrained  at 

the  sides.  This  was  followed  by  a  shear  type  of  deformation  in  which  the  rings 

were  crushed  to  deform  at  an  oblique  angle  to  the  direction  of  impact.  There  was 
considerable  lateral  expansion  of  the  rings  crushed  in  the  shear  mode.  The  sheer 
bands  started  at  the  corners  of  the  striker  and  formed  a  wedge  shape.  Continuing 
deformation  resulted  in  the  growth  (widening)  of  these  shear  bands.  Such 

behaviour  was  also  observed  in  quasi-static  indentation  tests  on  similar  systems. 
At  higher  impact  velocities,  i.e.  for  high  Q  values,  the  crushing  mode  dominated 
in  causing  deformation.  The  deformation  front  was  cylindrical.  Rings  adjacent 
to  the  side  walls  and  distal  end  support  were  found  to  deform  due  to  the 
reflection  at  these  supports  stress  waves  produced  during  impact. 

In  all  the  systems  tested  the  initial  decelerations  were  found  to  be  higher  than 
the  average  decelerations,  the  ratio  of  these  two  values  being  higher  for  higher 
values  of  Q.  Considering  the  strain  rate  effects  and  using  the  quasistatic 
failure  mode  of  rings,  an  elementary  energy  balance  procedure  is  adapted  to 
predict  the  variation  of  the  striker  velocity.  Considering  the  siaplicity  of  the 
analysis  the  predictions  compnre  reasonably  well  with  experimental  observations. 
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A  DESIGN  GUIDELINE  FOR  PREDICTING  MAXIMUM  CENTRAL  DEFLECTIONS 
OF  IMPULSIVELY  LOADED  PLATES 

G  N  Nurick  -  Applied  Mechanics  Research  Unit, 
University  of  Cape  Town. 

An  extension  of  Johnson ' a  ( 1 )  damage  number  is  refined  for  use 
as  a  design  guideline  for  predicting  the  maximum  central  def¬ 
lection  of  plates  subjected  to  impulsive  loads.  This  predic¬ 
tion  compares  favourably  with  experimental  evidence.  In 
addition,  using  the  experimental  results  an  upper  bound  on  the 
loading  conditions  may  be  determined,  i.e.  the  maximum  deflec¬ 
tion  before  tearing  of  the  plate  material  takes  place.  This 
is  illustrated  for  circular  and  rectangular  plates  of  varying 
dimension  and  loading  conditions. 


1.  Johnson  W  1972  Impact  Strength  of  Materials  p.303.  Edward 
Arnold. 


On  the  Effects  of  Materials  and  Radius  to  Thickness  Ratio  in  Dynamic 
Plastic  Buckling  of  Circular  Tubes  under  Axial  Impact 


K  Kawata,  S  Matsumoto  and  E  Yoshimizu 

Science  University  of  Tokyo,  Node,  Chiba-ken  278  Japan 

ABSTRACT:  In  dynamic  plastic  buckling  of  circular  tubes  under  axial 
impact,  the  effect  of  materials  in  the  mean  plastic  buckling  load 
(MPBL)  vs.  loading  velocity  relation  for  bcc  and  fee  metals  and  the 
effect  of  radius  to  thickness  ratio  in  MPBL  are  clarified,  using  a 
new  precise  dynamic  properties  characterization  system  based  upon 
the  KHKK  one  bar  method.  Experimental  values  are  fairly  explained 
by  a  theoretical  expression  derived  by  the  authors,  composed  of  three 
mechanisms  of  bending,  rolling  and  toroidal  deformation. 

INTRODUCTION 

In  designing  the  crashworthiness  of  structural  elements  such  as  circular 
tube  under  axial  impact,  the  effect  of  materials  in  the  MPBL  vs.  loading 
velocity  relation  and  the  effect  of  radius  to  thickness  ratio  in  MPBL 
should  be  known,  basing  upon  the  precise  experimental  result.  The  precise 
curve  of  dynamic  plastic  buckling  load  vs.  displacement  seems  not  be  hither¬ 
to  satisfactorily  determined.  Such  data  acquisition  could  not  be  carried 
out  until  precise  high  velocity  mechanical  properties  characterization 
system  with  no  reflected  wave  disturbance,  such  as  the  KHKK  one  bar  method 
(Kawata  et  al  1979)  appeared.  These  experimental  data  are  compared  with 
the  theoretical  expression  newly  derived  by  Kawata  et  al  (1988). 

CONCLUSIONS 

1)  A  new  theory  of  dynamic  plastic  buckling  of  circular  tube  with  wide 
range  of  R/t  under  axial  impact,  is  derived  basing  upon  a  model  with 
bending,  rolling  and  toroidal  deformation  mechanisms.  The  theoretical 
values  well  coincides  with  experimental  ones  for  specimens  of  R/t  ranging 
from  about  20  to  120  and  of  various  metals,  tested  here. 

2)  Dynamic  plastic  buckling  load  vs.  displacement  relation  can  be  satis¬ 
factorily  determined  by  tha  new  precise  dynamic  property  characterzation 
system  baaing  upon  the  KHKK  cne  bar  method  without  reflected  stress  wave 
disturbance. 

3)  Dynamic  MPBL  can  be  calculated  for  various  loading  velocity  by  con¬ 
sidering  the  dynamic  mechanical  properties  of  the  materials,  giving  fair 
coincidence  with  experiment. 

A)  The  effect  of  materials  in  the  MPBL  vs.  loading  velocity  relation  is 
clarified  for  both  of  fee  and  bcc  metallic  materials. 

5)  The  effect  of  R/t  ratio  in  MPBL  is  clarified  experimentally  and  is 
fairly  explained  by  the  above  mentioned  theory. 


Qooj 


Dynamic  Coupling  of  Plastic  Deformations  in  Multi-Component  Structures 


W  J  Stronge  and  D  Shu 
Department  of  Engineering 
University  c?  Cambridge 


ABSTRACT 


The  distribution  oZ  energy  dissipation  in  multi-component  structures 
that  are  deformed  by  impact  depends  or.  the  structural  arrangement,  the  relative 
compliances  and  the  distribution  of  mass  in  the  components.  A  collision 
against  a  rigid-perfectly  plastic,  multi-component  structure  results  in 
at  laaat  two  distinct  stages  of  deformation;  a  transient  followed  by  a  modal 
stage  of  deformation.  Deformation  during  the  transient  stage  ia  characterised 
by  travelling  hinges  or  hinge  lines  while  the  modal  stage  hnB  a  time-independent 
velocity  distribution.  The  total  impact  kinetic  energy  is  disbursed  throughout 
the  structure  during  these  stages  of  deformation.  In  general  the  fraction 
of  the  total  energy  dissipated  in  each  component  depends  primarily  on  the 
structural  arrangement  and  the  relative  compliance;  the  mass  of  tho  component 
is  of  secondary  importance.  The  structural  arrangement  can  be  classified 
as  either  a  parallel  or  series  system  in  analogy  with  electric  circuits. 

The  distinction  between  these  systems  is  based  on  displacement  constraints 
that  define  the  connectivity  between  elements. 

In  the  parallel  arrangement,  all  components  have  the  same  displacement 
at  the  impact  point;  the  dynamic  load  is  imparted  to  all  components  simultaneously. 
In  a  multi-component  structure  with  this  arrangement,  the  fraction  of  the 
total  energy  absorbed  in  a  component  during  the  transient  phase  increases 
with  the  relative  compliance  and  decreases  with  the  relative  mass  ratio. 

During  the  modal  stage  of  deformation,  the  fraction  of  the  total  energy 
dissipated  increases  in  proportion  with  the  relative  compliance.  The  deformation 
of  a  component  is  independent  of  the  relative  mass  ratio  provided  the  mode 
of  deformation  ia  not  changed  when  the  mass  ratio  changes.  To  conclude, 
the  heavier  and  stronger  components  of  a  structure  absorb  a  larger  part 
of  the  initial  kinetic  energy  in  a  parallel  structural  arrangement. 

In  the  eeries  r.rrangement,  only  the  dynamic  forces  that  have  passed 
through  an  element  are  transmitted  to  successive  elements.  With  this  arrangement 
the  fracture  of  the  energy  absorbed  in  a  component  during  the  transient 
stage  increases  with  relative  mass  and  decreases  with  relative  compliance. 

Again  the  stronger  components  absorb  a  larger  fraction  of  the  initial  energy 
imparted  to  the  structure. 

Applications  of  these  results  to  the  analysis  of  ship  structures  composed 
of  frames  covered  by  plates  will  be  discussed. 


Tha  Use  of  Simulation  Techniques  to  Study  Perforation  Meohanisms  In 
laminated  Netallio  Composites 


X.G.CrouchC*)  and  R.L. Woodward 

Materials  Researoh  Laboratory,  Defence  Soienoe  and  Technology  Organisation, 
Asoot  Vale,  Victoria  3032,  Australia 

EXTENDED  ABSTRACT 

As  a  class  of  oomposites,  metallic  laminates  offer  reoognised  benefits  in 
terms  of  structural  properties,  especially  enhanced  toughness  (Kaufman, 
1967?  Taylor  et  al,  1976?  Throop  et  al,  1979)  and  improved  fatigue 
performance  (Taylor  et  al,  1976;  Vogelsang,  1966).  This  is  particularly 
true  with  the  weakly-bonded  systems,  such  as  the  adhesively-bonded, 
aluminium  laminates.  The  addition  of  low-denalty,  polymerio  fibres  to 
relnforoe  the  M'tk  Interlayer,  as  in  the  case  of  ARAM,  a  hybrid  laminate 
ourrently  being  marketed  commercially  by  Alooa,  increases  the  specific 
strength  values  whilst  enhancing  fatigue  orr.ok  toleranoe  even  farther. 
These  attrsotlve  speoific  properties  have  led  to  a  number  of  full-soale 
trials  within  aircraft  structures  (Vogelsang,  1986),  for  vhloh  weight 
aavinge  of  about  30*  are  olelmed  possible.  Unfortunately,  little  is 
reported  of  the  through-thiokneos  properties  of  these  laminates,  especially 
under  impact  conditions,  for  it  is  these  properties  which  need 
investigating  before  auoh  structural  laminates  could  be  oonsidered  as 
suitable  materials  for  containment  applioations  or  in  situations  where 
perforation,  at  high  rates  of  strain,  is  a  possibility. 

The  energy-absorbing  c&paolty  of  metallic  laminates  has  reoently  been 
highlighted  by  Crouch  (1988)  who  compared  the  perforation  characteristics 
of  adhesively-bonded,  aluminium  laminates  with  medium-strength,  monollthio, 
aluminium  alloy  plates.  It  was  shown  that  suoh  laminates  suppress  failure 
by  the  through-thickness  shearing  mechanism  known  as  plugging  and 
consequently,  under  certain  conditions,  oan  absorb  more  kinetic  energy. 
The  perforation  of  metallic  laminates  by  blunt  objects  is,  however,  a 
complex  terminal  ballistio  event  and  attempts  to  mathematioally  model  the 
prooesa  hava  proved  extremely  difficult  because  little  is  known  about  the 
plastlo  flow  and  yield  criteria  associated  with  metallic  laminates.  A 
solution  for  simple,  unbended,  two-layer  laminates,  whioh  treats  both  the 
plug  acceleration  and  dishing  failure  aspects  of  the  process,  has  been 
developed  (Woodward  et  al,  1988)  by  examining  the  deformation  of  the  two 
component  layers  separately,  for  a  multi-layered  configuration  a  continuum 
approaoh  is  required,  necesaitating  either  a  better  understanding  of 
oompoaite  behaviour  or,  alternatively,  approaches  whioh  sample  properties 
in  a  form  similar  to  the  actual  stress  state  in  a  ballistic  experiment. 

The  cbjeot  of  this  work  was  to  establish  such  a  set  of  quasi-statio, 
meohanioal  teat  procedures  whioh  would  accurately  simulate  part  of  the 
oomplex  perforation  prooeas.  Data  and  derived  parameters  from  these  teats 
oan  then  be  used  as  realistlo  input  data  in  mathematical  models  of  laminate 
perforation. 

(*)  Current  address:  Royal  Armaments  Researoh  and  Development  Establishment, 
Chobham  Lane,  Cherteey,  Surrey,  ET16  OEE,  U.K. 


The  perforation  mechanisms  of  weakly-bonded  matallio  laminates  are  compared 
with  the  failure  mode  of  monollthlo  aluminium  plates  when  impaoted  with 
blunt-nosed  projeotlles.  Figures  la  and  1b  illustrate  the  comparative 
damage  features  associated  with  both  the  monolithic  and  laminated 
materials:  (A)  compression,  (B)  shearing,  (C)  tensile  stretching  and  (D) 
delamination.  A  family  of  quasi-static  simulation  tests  have  been 
developed  for  laminates  which  generate  relevant  and  meaningful  meohanioal 
property  data  in  a  form  suitable  as  input  data  for  mathematical  modelling 
of  the  perforation  prooess.  Figures  2a  and  2b  show  how  constrained 
compression  tests  have  been  ueed  to  simulate  the  compression  stage  in  the 
perforation  of  both  olaaaes  of  materials  by  blunt-nosed  projeotlles. 


The  properties  described  in  this  work  are  ourrently  being  utilised  in  a 
simple  sequential  model  whioh  treats  plug  aooeleration  and  target 
delamination  as  two  basic  stages  in  the  perforation  of  laminated  metallio 
composites.  Using  data  derived  from  these  simulation  tests,  good  estimates 
of  oritloal  ballistic  Impaot  velooities  are  Deing  obtained  and  the  effeots 
of  parameter  changes  on  performance  are  being  oorreotly  described. 
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Instabilities  at  High  Strain  Kate  of  AISI  316  Damaged  by  Creep,  Fatigue 
and  Irradiation.  Discussion  of  Adiabatic  Temperature  Oscillations. 

C.  Albertlnl,  A.M.  Elelche*  and  M.  Montagnanl 

CEC-JRC,  Ispra  Establishment,  21020  Ispra  C Va )  -  Italy 

“Visiting  scientist  at  JKC-Ispra,  on  leave  from  the  Cairo  University 

Nuclear  reactor  components  are  submitted,  In  the  case  of  a  hypothetical 
severe  accident,  to  dynamic  loading  generating  strain  rates  of  up  to  1000 
s~1.  The  structural  behaviour  during  severe  accidents  can  be  assessed 
provided  that  the  mecnanlcal  properties  of  the  material  at  the  time  of  the 
accident  can  be  predicted  by  constitutive  equations.  Therefore,  the 
tensile  mechanical  properties  of  the  nuclear  grade  austenitic  stainless 
steel  AISI  316H  were  measured  at  strain  rates  ranging  from  10-3  to  103 
s-1 ,  which  steel  was  previously  submitted  to  Increasing  amounts  of  pre¬ 
damage  by  creep,  low  cycle  fatigue  and  Irradiation  at  a  temperature  of 
550*C. 

The  following  programme  of  pre-damage  treatments  was  performed  on  the  as- 
received  AISI  31 6H  stainless  steel. 

-  Irradiation  to  2  dpa  (E  >  0.1  MeV)  In  the  High  Flux  Reactor  In  Petten, 

In  sodium  at  550°C; 

-  low  cycle  fatigue  at  550*C  (performed  by  P.N.C.  Tokyo  at  the  University 
of  Tckyo),  cycle  strain  range  J.6%  and  1%,  cycle  ratio  n/NF  of  0.2,  0.4 
and  0.6; 

-  Low  cycle  fatigue  as  above  with  superposition  of  Irradiation  t.o  2  dpa; 

-  creep  at  550°C,  300  MPa  (performed  at  JKC-Ispra),  at  the  fractions  of 
0.2,  0.4  and  0.6  of  the  2000  h  life-time. 

From  the  typical  flow  curves  regarding  creep,  low  cycle  fatigue,  low  cycle 
fatigue  and  Irradiation  damaged  materials,  respectively,  we  observe  that 
the  single  damage  process  (creep  alone,  low  cycle  fatigue  alone, 
Irradiation  alone)  Is  a  hardening  process  which  provokes  an  Increase  In 
the  flow  stress  at  a  given  strain  and  a  reduction  of  ductility  with 
respect  to  the  as-received  material.  We  observe  also  that  up  to  large 
strain  values  the  strain  rate  Increase  causes  a  flow  stress  Increase  at  a 
given  strain  In  each  of  the  material  conditions  considered  alone.  Dynamic 
strain  ageing  serrations  are  observed  along  the  low  strain  rate  curves  of 
as-received  and  damaged  materials,  while  at  high  strain  rate,  larger 
instabilities  of  the  stress-strain  curve  of  the  damaged  materials  are 
present.  Furthermore,  we  observe  that  the  superposition  of  Irradiation  and 
low  cycle  fatigue  enhances  the  effect  of  flow  stress  Increase  at  a  given 
strain  and  that  of  ductility  reduction  with  respect  to  the  as-received 
material,  and  Introduces  a  trend  to  strain  rate  softening. 

In  an  attempt  to  find  a  correlation  between  the  mechanical  response  of  the 
damaged  materials  and  the  amount  of  damage,  two  parameters  were  defined. 
One  parameter  (0.2%  yield  stress  ratio  of  damaged  to  as-received  material) 
measures  the  Increase  of  flow  stress  with  respect  to  the  as-received 
material.  The  other  parameter  (fracture  strain  ratio  of  damaged  to  as- 
received  material)  measures  the  reduction  of  ductility  with  respect  to  the 
as-received  material.  The  numerical  values  of  these  two  parameters  show  a 
correlation  which  Is  worth  further  verification  because  Its  consistency 
would  permit  a  simplified  Implementation  of  the  effects  of  different 
damage  processes  In  the  constitutive  equations.  Nevertheless,  a  more 
complete  theoretical  modelling  of  the  damaged  material  properties  for 
structural  calculations  must  take  Into  account  also  the  instabilities 


which  modulate  the  averaged  mechanical  response  of  the  damaged  materials 
with  oscillations  of  the  flow  curves.  The  amplitude,  frequency  and  nature 
of  these  oscillations  depend  on  strain  rate,  temperature  and  stress  state. 
The  low  strain  rate  flow  curve  of  the  creep-subjected  material  shows  the 
serrations  of  dynamic  strain  ageing.  Some  Informations  about  the  nature  of 
dynamic  strain  ageing  can  be  deduced  by  observing  the  records  of  the 
stress-tline  and  of  the  temperature  Increase  of  the  specimen-time  which 
show  osscll latlons  at  corresponding  times.  This  fact  means  that  dynamic 
strain  ageing  Is  a  thermally  activated  phenomenon. 

The  medium  strain  rate  flow  curve  of  the  creep-subjected  material  Is 
smooth  without  oscillations;  also  the  stress-time  and  the  temperature 
Increase  of  the  specimen-time  curves  are  smooth  without  oscillations. 

The  high  strain  rate  flow  curve  of  the  creep-subjected  material  shows  an 
Initial  upper  and  lower  yield  point  followed  by  large  oscillations  due  to 
Instabilities.  Both  the  stress-time  and  the  temperature  Increase-time 
curves  at  high  strain  rate  show  also  oscillations  at  the  same  time.  The 
oscillations  of  the  high  strain  rate  flow  curve  are  assisted  by  thermal 
activation.  The  constraint  to  slip  due  to  the  short  time  and  the 
consequent  dislocation  locking  are  responsible  for  the  remarkably  high 
strain  rate  hardening  with  superimposed  large  Instabilities. 

These  facts  could  mean  that  the  straining  modes  at  grain  level,  developed 
during  the  creep  and  tensile  tests  at  low  strain  rate,  are  of  a  similar 
nature,  mainly  by  micro-flow  along  grain  boundaries,  while  they  develop 
new  mechanisms  at  high  strain  rate,  characterized  by  slip  Cor  twinning) 
Inside  the  grains.  The  same  phenomenon  Is  also  confirmed  by  the 
experimental  evidence  that  the  monotonic  tensile  curves  at  low  and  high 
strain  rate  coincide  with  the  correspondlntg  flow  curves  of  the 
Interrupted  test  performed  at  the  same  strain  rate.  Also  the  strain  rate 
softening  and  the  higher  ductility  reduction  at  low  strain  rate,  shown  by 
the  pre-fatlgued  and  Irradiated  material,  can  be  explained  as  being  due  to 
different  mechanisms  of  deformation  at  low  and  high  strain  rate. 

A  first  attempt  to  establish  an  equation,  based  on  atomistic  models,  able 
to  predict  the  onset  of  serrated  yielding,  taking  into  account  the  strain 
and  strain  rate  histories,  was  made  by  Elelche  et  al.  C19a7). 

Nevertheless,  this  approach  does  not  take  Into  account  the  directionality 
dependence  Introduced  by  the  stress  state.  In  fact,  recent  experiments  on 
the  same  AISI  316H  tested  in  tension  at  550"C,  showed  the  disappearance  of 
the  dynamic  strain  ageing  Instability  along  the  flow  curves  determined 
under  torsional  stress  state.  In  the  same  experiments  campaign  It  was 
stated  that  under  torsional  stress  state  by  Increasing  the  damage  through 
sequential  reverse  prestraining,  the  yield  stress  decreases  and  the 
residual  ductility  Increases..  These  results  In  torsion  are  compared  with 
those  obtained  In  tension  and  here  reported,  from  which  they  differ 
substantially. 
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Extended  Abstract 


The  present  paper  attempts  to  investigate  direct  shear  failure  in  the 
roof  slabs  of  shallow-buried  RC  boy-type  structures  (  Figure  1  )  under 
impulsive  pressure  loading.  Since  direct  shear  failure  is  governed  by  the 
buildup  of  the  transverse  shear  force  at  the  roof  support,  an  elastic 
Timoshenko  bean  model  is  adopted  for  the  RC  roof  slab.  The  finite  element 
equations  of  motion  for  the  Timoshenko  beam  are  derived  from  Hamilton's 
principle  and  solved  using  a  step-by-step  Integration  method. 

The  finite  element  procedure  was  applied  tc  obtain  the  dynamic  response 
of  the  RC  slabs  or  beams  under  distributed  impulsive  pressures.  The 
experimental  data  used  for  the  analysis  come  from  a  series  of  the  dynamic 
tests  on  the  shallow-buried  RC  box-type  structures  conducted  by  Kiger  and 
Slawson  (1984),  Each  test  structure  was  covered  with  a  shallow  layer  of 
soil  and  was  loaded  with  a  high-explosive  induced  blast  pressure  p(t), 
which  was  approximated  by  a  triangular  pulse  with  a  zero  rise  time.  The 
main  characteristics  of  the  example  RC  beam  are  summarized  in  Table  1. 
Because  of  symmetry,  only  half  of  the  RC  beam  was  considered. 

Figures  2(b)  and  2(c)  show  plots  of  normalized  support  shear  (V/Vu)  and 
normalized  support  moment  (M/Mu)  versus  time  for  two  different  peak 
pressures  for  the  example  RC  beam  with  fixed  ends.  The  present  results 
agree  reasonably  well  with  the  normal  mode  s'’utions.  It  should  be 
remarked  that  the  curves  in  Figure  2  are  theoretical  because,  as  soon  as 
either  ratio  V/Vu  or  M/Mu  exceeds  a  value  of  one,  the  RC  beam  is  presumed 
to  have  reached  its  strength  capacity  and  is  no  longer  elastic.  The  time 
at  which  V/Vu  ■  1  is  denoted  as  t8,  and  the  time  at  which  M/Mu  »  1  is 
denoted  as  t^.  The  occurrence  of  failure  in  either  direct  shear  or 
flexure  depends  on  whether  the  failure  threshold  (V/Vu  and  M/Mu  *  1) 
Intersects  these  curves  at  an  early  or  a  later  time,  respectively. 

The  transition  from  a  predicted  direct  shear  failure  to  no  shear  failure 
occurs  when  t8  “  tm  .  If  this  transition  for  different  combinations  of 
peak  pressure  Pg  and  rise  time  tf  is  plotted  in  the  Pg-t  domain,  a 
series  of  points  produces  a  "failure  curve.”  Typical  failure  curves 
determined  for  two  beam-end  conditions  are  given  in  Figure  3.  Points  that 
lie  in  the  region  above  a  failure  curve  define  a  loading  condition  for 
which  analysis  Indicates  a  direct  shear  failure,  while  points  that  lie  in 
the  region  below  the  failure  curve  describe  loading  parameter 
combinations  which  will  cause  either  a  flexural  failjre  or  no  failure  in 
the  RC  beam.  Parametric  studies  demonstrate  that  the  failure  curves  are 
sensitive  to  loading  parameters  (  Pq  and  t_  )  and  structural  parameters  ( 
R,  ft  and  L/d  ),  but  are  not  significantly  affected  by  load  duration  t,j. 
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Figure  If  Failure  modes  of  RC  box-type  structures  subjected  to  transverse, 
uniformly  distributed,  impulsive  loads:  (a)  flexural  damage;  (b)  combined 
flexural  ana  shear  damage;  (c)  collapse  :nder  direct  shear 


Table  1:  Main  characteristics  of 
example  RC  beam 


Parameter 

value 

Beam  length  (L) 

1.14  m 

Beam  width  ( b) 

0.0251.  m 

Beam  dt  ih  (h) 

0.184  m 

Effective  beam  depth  (d) 

0.162  m 

Young' s  modulus  (E) 

33.1  GPa 

Poisson's  ratio  (v) 

0.2 

Shear  coef f i c ient  (k ' ) 

0.822 

Hass  density  (q) 

2,600  kg/m5 

Ultimate  shaar  capacity  (Vu) 

6,160  kg 

Ultimate  moment  capacity  (Mg) 

370  kg-m 

Strength  enhancement  factor  (ft) 

1.6 

Figure  2:  Normalized  support  shear  and  Figure  3:  Effect  of  beam-end 
moment  versus  time  for  example  RC  beam  restraint  on  failure  curve 
:  (a)  fixed  ends;  (b)  peak  pressure  Pq 
■  35  MPa;  (c)  peak  pressure  Pq  *  14  MPa 


Dynamic  properties  of  construction  materials  using 
a  large  diameter  Kolsky  bar 
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ABSTRACT 

High  strain  rate  testing  at  the  University  of  Sheffield’s  Dynamics  Testing 
Laboratory,  near  Buxton,  Derbyshire,  has  been  developed  using  38mm  diameter 
maraging  steel  pressure  bars  and  a  small  explosive  charge  to  produce  the 
the  stress  pulse.  The  large  diameter  pressure  bars  are  of  particular 
interest  to  engineers  because  they  admit  the  use  of  thicker  discs  of 
specimen  in  the  test,  on  the  basis  of  Davies  and  Hunter’s  (1963)  geometric 
criterion. 

The  diameter  of  the  pressure  bar  was  chorren  as  38mm  for  the  purpose  of 
testing  construction  materials  with  particle  sizes  up  to  2mm.  The  pressure 
bars  were  aligned  vertically  and  supported  in  a  steel  frame.  The  dynamic 
strains  produced  in  the  input  and  output  pressure  bars  were  measured  using 
four  1mm  electrical  resistance  strain  gauges  (ERSG)  connected  to  a 
Wheatstone  bridge,  and  the  amplified  output  was  captured  on  a  Gould  OS  4050 
digital  storage  oscilloscope  (see  fig  1)  and  subsequently  transferred  to  an 
Olivetti  M24  personal  computer. 
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fig  1  i  Details  of  the  38mm  diameter  DTD5212  Kolsky  bar  apparatus 


The  recorded  stress  pulses  were  corrected  for  dispersion  using  an  FFT  with 
Bancrofts  dats.  The  correction  gave  very  good  agreement  with  experimental 
data. 


The  rod  velocity  of  specimens  was  found  using  dynamic  photoelasticity  to 
observe  the  progress  of  a  stress  pulse  through  the  specimen  (See  fig  2). 
The  camera  used  was  the  Barr  and  Stroud  CPS  rotating  minor  camera,  which 
can  operate  at  2,000,000  frames  per  second. 

This  method  gave  a  useful  visual  check  on  the  progress  of  a  stress  pulse 
through  a  spicimen  but  the  experiment  required  a  camera  setting  which  gave 
an  interframe  time  of  1.9  microseconds.  The  accuracy  of  the  rod  velocity 
found  using  this  method  was  limited  by  the  camera  speed. 


fig  2  i  Photoelastic  method  for  finding  the  tod  velocity  of  specimens 

Stress  /  strain  in  the  specimens  was  found  using  Lindhom  and  Yeakleys 
method.  Initial  tests  were  carried  out  on  Perspex  discs,  (38mm  diamter)  of 
different  height.  The  yield  stress  and  yield  strain  for  different  heights 
of  Perspex  were  plotted  together,  and  it  can  be  seen  (fig  3)  that  a  marked 
change  of  behaviour  of  the  specimen  occurs  between  the  12.5mm  ^ud  15mm 
specimen  heights.  The  critical  height  of  the  specimen  marks  the  transition 
between  two  modes  of  specimen  behaviour  as  a  result  of  the  influence  of 
boundary  conditions.  The  smaller  specimens  failed  at  greater  stress  and 
strain  values,  because  the  friction  coefficient  of  the  Perspex  cylinder 
ends  was  more  significant  than  radial  inertia.  Using  the  correct  height  of 
specimen  (14.5mm)  the  error  due  to  friction  was  of  the  order  of  +5Z  yield 
stress.  As  the  height  of  the  specimen  was  increased  beyond  14.5mm,  the 
inertia  term  due  to  radial  strain  acceleration  for  the  specimen  increased, 
reducing  yield  stress 


allllilMla 

°  5  *  to  *  It, a  a  is  x  to  t  ao 


fig  3  i  Yield  stress  and  yield  strain  calculated  from  Kolsky 
bar  data  for  specimens  of  different  height. 

A  number  of  construction  materials  were  tested  in  the  Kolsky  Bar,  and  two 
sets  of  results  are  presented  in  the  paper  i  a  high  quality  brick  pavior, 
and  sand /cement  mortar. 


D/oJ 


ULTIMATE  STRENGTH  OF  BOX  COLUMNS 
SUBJECTED  TO  IMPACT  LOADING 

to  be  pretested  at  the  Fourth  Oxford  Conference 
on  Mechanical  Properties  of  Materials  at  High  Rates  of  Strain, 
University  of  Oxford,  March  19-22,  1909 

by 

Tomass  Wierzbicki 
Department  of  Ocean  Engineering 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139 
nnd 

Leonhard  Recke 

BMW,  Research  and  Development  Center 
8000  Munich  40 


Atotract 

Buckling  of  plates  and  shells  under  dynamic  loading  have  been  studied  in  the  literature 
for  several  decades.  A  comprehensive  literature  review  on  this  subject  was  published, 
among  others  by  Singer  (1981),  Lindberg  and  Florence  (1986)  and  Jones  (1989).  It  was 
found  that  the  dynamic  buckling  load  increases  with  the  rate  of  loading  and  is  always 
greater  than  the  corresponding  static  buckling  load.  In  the  case  of  plates  subjected  to  a 
mass  impact  the  rise  in  the  plate  stiffness  was  found  to  be  attributed  to  the  transverse 
inertia  of  the  plate. 

A  related  problem  of  the  ultimate  strength  of  plates  under  longitudinal  impact  does  not 
appear  to  have  been  studied  in  thu  literature.  It  is  expected  though  that  the  trend  will 
be  similar  to  the  case  of  a  buckling  load.  The  problem  of  the  determination  of  the 
maximum  strength  of  compression  members  under  impact  loading  is  of  great 
importance  in  many  industrial  applications,  most  notably  in  the  impact  energy 
absorption  and  collision  protection  of  vehicles.  For  example,  thin-walled  box  columns 
loaded  statically  exhibits  a  sharp  peak  load  followed  by  subsequent  lower  peaks 
corresponding  to  folding  and  crushing  of  the  column’s  walls.  The  magnitude  of  the 
peak  load  (often  referred  in  the  literature  as  ultimate  or  crippling  load)  depends  on  the 
cross-sectional  dimensions  of  the  column,  elastic  constants  and  the  yield  stress  of  the 
material.  The  peak  load  of  the  same  column  loaded  dynamically  is  larger. 

The  objective  of  this  communication  is  to  present  a  simplified  method  of  predicting  the 
ultimate  strength  of  prismatic  multi-cornered  columns  subjected  to  impact  loading.  The 
impact  velocities  considered  in  this  paper  are  relatively  low  and  correspond  to  typical 
collision  velocities  of  automobiles  and  airplanes.  The  impacting  mass  is  assumed  to  be 
larger  than  the  total  mass  of  the  column. 

The  analysis  consists  of  two  parts.  First,  a  new  two-degree-of-freedom  deformation 
model  of  a  single  corner  element  of  the  multi-cornered  column  is  developed  and  the 


post -buckling  stiffness  and  ultimate  load  is  evaluated.  This  is  derived  in  a  closed  form 
from  the  first  yield  criterion.  The  calculated  strength  of  n  perfect  plate  agrees  within 
IS  %  with  the  exact  solution  (Rhodes  I98S).  The  accuracy  is  further  increased  by 
taking  into  account  the  effect  of  initial  imperfections. 

In  the  second  part  of  the  analysis,  the  lateral  inertia  of  the  column  wails  is  incorporated 
into  the  governing  equation.  The  axial  inertia  is  neglected  for  the  considered  class  of 
problems  (large  mass,  low  velocity).  The  formulation  is  based  on  the  concept  of  the 
growth  of  initial  imperfection  with  time.  The  resulting  initial  value  problem  for  a 
system  of  ordinary  differential  equations  is  solved  numerically.  It  is  shown  that  the 
dynamic  ultimate  strength  depends  on  the  magnitude  of  the  striking  mass  as  well  as  on 
the  impact  velocity.  A  correlation  of  the  predicted  strain-time  diagrams  are  made  with 
the  test  results  reported  by  Singer  (1981). 

The  present  results  will  help  interpret  the  drop  tower  test  on  box  columns  performed  in 
many  laboratories  around  the  world.  In  this  type  of  experiment,  the  force-time  history 
at  the  mass/column  interface  is  determined  from  the  reading  of  the  accelerometers.  The 
readings  are  usually  affected  by  the  inaccuracies  in  differenciating  the  displacements  to 
get  accelerations  and  by  the  superposed  elastic  vibrations  of  the  system.  Therefore,  the 
experimentally  determined  peak  load  is  not  uniquelly  defined  and/or  is  affected  by  the 
filtering  frequency.  Our  analysis  provide  a  much  needed  insight  into  the  physics  of  the 
impact  problem. 

At  the  present  stage  we  have  studied  the  inertia  effect  only.  The  peak  load  is  also 
influenced  by  the  strain  rate  effect  and  the  delayed  yield  phenomenon  (Campbell 
1957).  Consideration  of  the  above  factors  is  left  to  a  future  continuation  of  this 
rssearch. 
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ABSTRACT 

A  raw  experimental  facility  is  described.  This  comprises  a  two  stage 
light  g as  gun  and  associated  diagnostics  which  will  be  used  to  carry  out 
impact  experiments  and  measure  high  strain  rate  properties  of  Materials. 
The  gss  gun  launch  tube  has  sn  8mm  bore,  firing  into  an  evacuated 
experlawmtal  chamber  and,  to  date,  the  saxisua  recorded  velocity  is  3*2 
km /s  using  a  one  gram  projectile.  Diagnostic  equipment  Includes  high 
speed  photography,  impact  pressure  measurement  with  piezoelectric  gauges 
connected  to  a  fa at  digitiser  and  projectile  velocity  measurement. 

Areas  of  research  Include  debris  distribution  and  fragmentation  of 
targets  and  projectiles,  accuracy  and  applicability  of  existing  high 
strain  rate  data  end  eicrostrucfcural  effects  in  dynamic  fracture.  A 
schematic  of  the  gun  is  shown  below. 
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ABSTRACT 


This  paper  addresses  the  problem  of  constructing  a  numerical  model 
of  the  entire  hydraulic  shocK  and  damage  sequence  of  events.  This 
embraces  the  penetration  of  a  supersonic  missile  fragment  into  a  thin- 
walled  structure,  followed  by  the  deceleration  of  the  projectile  in  the 
fluid.  The  cavitation  is  modelled  and  the  consequent  pressure  rise  due 
to  fluid  compression.  The  constitutive  law  is  discussed.  Finally  the 
fluid  pressure  is  applied  to  the  structure  (or  a  pressure  alleviation 
device)  and  the  consequent  response  modelled. 

The  numerical  and  financial  consequences  of  modelling  an  entire 
three-dimensional  wing-box  are  stated,  together  with  a  cheaper  alternative. 

Finally  the  sparsity  of  data  is  outlined  for  representing  the 
composite  matrix  toughness  and  its  dependence  on  strain  rate.  The 
evidence  indicates  a  rapid  deterioration  in  resistance  to  interlaminar 
shear  failure  with  strain  rates  in  excess  of  100  per  second. 
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ABSTRACT:  This  paper  presents  a  simple  analytical  model  for  predicting  the  impact 
behaviour  of  Composite  Material  armours.  It  can  be  used  to  calculate  trie  armour  ballistic 
curve  and  other  data  up  to  the  time  of  penetration,  the  total  force  retarding  the  projectile, 
the  tension  in  each  layer,  the  displacements  and  velocities  of  the  plate  and  projectile,  the 
strain  rate  supported  by  each  layer,  and  the  stresses  and  strains  in  the  plate . 


1.  TECHNIQUE 

1 . 1  Behaviour  of  single  yams  under  transverse  impact 

The  problem  considered  is  as  follows:  a  rigid  projectile  of  mass  m  and  velocity  v  impacts  a 
plate  of  lengh  L,  width  W  and  thickness  H  and  made  of  n  layers  of  fabric. 


The  stress  analysis  is  based  on  the  general  assumption  that  the  fibre  behaviour  is  linear  elastic 
up  to  the  fracture  (Harding,  1987). 

When  die  projectile  reaches  the  n+1  layer  the  equation  system  is  (Parga, 73.,  1988): 
vn  =  vn-l  '  an(tn  '  1 ) 


an  * 


(n-1  +  1)  h 

c  (tn  •  tj. i) 


Ti  =  2|£(tn-‘n.r2ti-i)E 


ri  t  (n-i+l)  h  ,2  j  (  (n_if  Q  h  Y»-t 

LJV  c  (!„-£,) )  c(tn-tn.7H  J 


h  *  vn-l  On  "  *n-l) '  On  '  ln-l)^ 


which  is  a  system  of  n+3  equations  and  n+3  unknowns:  vn,  an,tr  and  Tj,  i  =  l,..,n 
The  fracture  criterion  adopted  has  been  that  of  maximum  strain. 

In  the  case  of  Composite  Material  armours  it  is  assumed  that  as  soon  as  the  projectile  strikes 
the  impregnated  yam,  fibre  debonding  is  produced.  So  two  deformation  waves  will  propagate 
one  through  the  yam  and  the  other  through  the  resin.  When  yarns  are  impregnated  with  a 
matrix  content  below  the  critical  content,  the  resultant  wave  speed  proposed  is: 


0 >n 

.  .  .  _  / 


C  *  Cjh  +  ^^(Vf+K-l) 

where:  Vf  1-K,  cm  is  the  matrix  wave  velocity,  Cf  is  the  fibre  wave  velocity,  K  is  the  critical 
resin  content  per  unit,  Vf  is  the  fibre  volume  fraction  per  unit, 

1.2  Behaviour  of  fabrics  under  transverse  impact 

If  the  only  mechanism  operating  between  the  primary  and  the  secondary  yams  is  friction,  the 
relation  between  the  primary  ycm  tension  supporting  th  ;  projectile  Ti  and  the  tension  T 

supported  by  the  same  yam  after  n  secondary  yams  are  passed  is1  T  *»  Tj  e‘PPn  whcre:p  is 

the  friction  coefficient  between  primary  and  secondary  yams,  P  is  the  coutact  angle  between 
primary  and  secondary  yams,  n  is  the  number  of  secondary  yams  considered.  This  fact 
modificaKs  the  total  area  affected  by  the  impact. 

2.  RESULTS 

To  take  experimental  scatter  into  account  in  determining  the  dynamic  modulus,  simulations 
were  performed  using  two  different  values.  Ballistic  limit  curves  were  obtained,  looking  for 
the  lowest  velocity  that  perforates  an  armour  plate  for  a  given  surface  density.  Then  »he  surface 
density  of  the  armour  was  increased  by  half  a  unit  and  the  process  repeated. 


Fig.  Ballistic  curves  vs-  surface  density  for  a  Kevlar  29  soft  armour  plate  when  impacted  by 
a  9mm  Parabellum  projectile.  Triangles  are  the  experimental  values  of  V50  (Dupont  1983b). 


3.  CONCLUSIONS 

1.  The  impact  behaviour  of  soft  armour  plates  is  simulated  by  a  unidimensional  model 

2.  The  model  permits  the  simulation  of  the  penetration  process  into  Composite  Materials  with 
much  shorter  computer  times  than  with  conventional  impact  codes. 

3.  Predicted  ballistic  curves  match  experimental  values  fairly  well. 
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Structures  consisting  of  two  or  more  layers  joined  by  rivets  or 
adhesively  bonded  are  common  in  aeroengines  and  airframes.  Under  normal 
service  conditions  the  shear  stresses  between  adjoining  layers  are  usually 
considerably  lover  than  tha  tensile  or  compressive  stresses  carried  by  the 
layers  themselves  and  the  only  mode  of  failure  that  needs  to  be  considered 
by  designers  is  the  one  associated  vith  direct  stresses.  Impact  loading  on 
the  other  hand  is  capable  of  inducing  relatively  high  shear  stresses  that 
may  cause  tha  layers  to  separate  vith  a  consequent  reduction  of  residual 
strength  after  the  impact.  In  this  report,  the  problem  is  illustrated  by 
means  of  dynamic  photoelasticity,  using  a  tvo-layer  contilevar  beam  that  is 
analysed  by  solving  the  stress  vave  equations  directly  and  by  the  ABAQUS 
computer  code.  It  is  shown  that  the  ratio  of  shear  stress  to  tensile  stress 
can  be  about  four  times  higher  under  dynamic  conditions  than  under  static 
conditions.  Frictional  effects  are  also  shown  to  be  important  in 
dissipating  some  of  the  kinetic  energy  of  the  impactor.  The  effects  of 
these  findings  on  the  design  of  dual  layered  shields  are  finally  discussed. 
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Abstract  -  An  experimental  method  involving  coupled  pressure  bars  (CPBs)  [1], 
preloaded  in  tension,  has  been  used  to  measure  dynamic  fracture  initiation, 
propagation  and  arrest  properties  of  A533B  pressure  vessel  steel  at  temperatures  up  to 
50°C  (77°C  above  the  Nil  Ductility  Temperature).  Two  compact  specimens  are 
simultaneously  loaded  by  releasing  the  energy  stored  in  the  pressure  bars  by  the 
controlled  fracture  of  a  brittle  notched  specimen.  Loading  characteristics  and  crack 
opening  displacements  are  deduced  by  analyzing  strain  histories  measured  on  one 
pressure  bar,  as  well  as  by  using  eddy  current  displacement  transducers.  The  dynamic 
crack  initiation  event  is  characterized  by  a  stress  intensity  rate  (k.)  of  3  x  106 
MPa  dn  s  .  The  crack  then  rapidly  propagates  under  a  constant  crack  opening 
displacement  rate  inducing  crack  velocities  ranging  from  200  to  600  m/s.  A  precise 
analysis  of  the  experiment  has  been  performed  using  a  dynamic  viscoplastic  finite 
element  computer  program.  The  resolution  of  the  procedure  has  been  investigated  by 
obtaining  replicate  measurements  of  three  quantities  versus  time-speeiflcally,  crack 
opening  displacement,  at  two  locations,  and  crack  growth.  The  well-defined 
displacement  boundary  conditions  in  the  experiment,  along  with  the  analysis 
procedure,  enable  accurate  evaluation  of  dynamic  crack  initiation  and  propagation 
toughness  of  high  toughness  material  using  small  specimens. 


[1]  H.  COUQUE,  S.J.  HUDAK,  JR.  and  U.S.  LINDHOLM,  J.  Pbys.,  Colloque  C3,  Sup. 
n°9,  T49,  C3-347,  1988. 
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Abstract.  The  paper  alms  at  the  description  of  the  Influence  of 
thermal  effects  on  shear  band  localization  failure. 

In  technological  dynamical  processes  failure  can  arise  as  a 
result  of  an  adiabatlo  shear  band  localization  generally  attributed 
to  a  plastio  instability  generated  by  thermal  softening  during 
plastic  deformations.  Recent  experimental  observations  /of. 
H.H.Orebe,  H.-R.Palc  and  M.A. Meyers,  Me  tall.  Trans.  ,  1 6  A,  76 1-775 . 
1985;  K. A. Hartley,  J. Duffy  and  R.H. Hawley,  J. Mech.Phys. Solids , 

21>  283-301,1987/  have  shown  that  the  shear  band  prooreates  in  a 
region  of  a  body  deformed  where  the  resistance  to  plastio  de¬ 
formation  is  lower  and  the  predisposition  for  bond  formation  is 
higher.  Shear  bands  nuoleate  due  to  the  presence  of  a  looal 
inhocogeneity  or  defects  oausing  enhanced  looal  deformation  and 
looal  heating.  The  shear  band  once  prooreated  behaves  differently 
than  matrix  global  body. 

.  Basing  on  available  experimental  observation  results  an  analysis 
of  Important  phenomena  for  proper  description  of  the  final  failure 
during  dynamic  deformation  processes  have  been  investigated.  It 
has  been  found  that  such  phenomena  as  the  inelastio  flow  process, 
the  Instability  of  flow  prooess  along  localized  shear  bands, 
the  micro-damage  process  whioh  oonsists  of  the  nuoleation, growth 
and  ooaleaoenoe  of  miorooraoks  and  the  final  mechanism  of  failure 
are  very  much  influence  by  thermal  effeots. 

This  conclusion  suggests  to  consider  thermomeohanioal  coupling 
and  to  treat  the  dynamio  deformation  process  as  adiabatic. 

It  has  been  postulated  that  the  response  of  a  material  within 
shear  band  region  is  different  than  in  adjacent  zones.  Within  shear 
band  region  the  deformation  prooess  is  characterized  by  very 
large  strain  rates  and  very  serious  temperature  changes. 

During  the  first  stage  of  the  flow  prooess  a  thermo-elastio- 
plastio  model  of  a  material  with  internal  mioro-damage  effeots 
is  applied.  A  criterion  for  localization  of  the  plastio  defer- 
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nation  within  shear  banda  la  examined.  Particular  attention  la 
focuaaed  on  the  Investigation  of  thermal  effects  and  Induced 
anisotropy  on  localisation  phenomenon. 

Postorltioal  behaviour  within  shear  banda  is  modelled  by  a 
thermo- elastic- viaooplaotlo  response  of  a  material  with  advanoed 
micro-damage  process. 

Both  models  are  developed  within  the  framework  of  thermodynamic 
theory  of  the  rate  type  with  internal  state  variables.  The  main 
role  in  these  models  is  played  by  a  set  of  the  internal  state 
variable#.  This  set  oosists  of  the  lsotropio  hardening-sof toning 
parameter  ot  ,  the  residual  stress  tensor  oc  and  the  porosity  ^  , 

A  set  of  the  evolution  equations  proposed  for  the  shear  band 
region  is  temperature  and  rate  dependent. 

Both  models  introduced  are  consistent  with  the  requirement 
that  the  thermo-alas  t,lo-plast±o  response  of  a  porous  solid  can  be 
obtained  as  a  limit  aase  for  quasi-statio  processes  of  the  thermo- 
-elastio-viaooplaatio  behaviour. 

The  dynamic  failure  criterion  within  shear  bands  is  proposed. 

The  anisotropic  softening  of  the  material  implied  by  thermal 
effects  and  the  micro-damage  process  is  described.  A  simple  mioro- 
meohanioal  model  of  final  failure  with  influence  of  thermal  and 
anisotroplo  effeots  is  proposed.  Discussion  of  the  mechanical 
results  is  given. 
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The  plastic  constitutive  and  damage  laws,  as  well  as  the  fracture  toughness  of  a  35 
NCD16 steel  (0,378 C-  0,39 SI  -0.! 6 (In-  0.003 S  -  0.008  P  -  4~1I  Ml  -  1.82  Cr 
-  0.44  Mo  -  In  Wfc )  heat  treated  to  3  different  yield  strength  to  ductility  ratios,  hove 
been  studied  under  tensile  ctynamlc  loading  .  Tests  have  been  performed  on  a  dynamic 
tensile  apparatus,  the  principle  of  which  is  based  on  that  of  split  Hopkinson  bars.  Smooth, 
notched  and  fatigue  precracked  axlsymmetricai  specimen  heve  been  used  .  8oth  the  load  at 
the  specimen  ends  and  the  notch  opening  displacement  or  the  specimen  elongation  have  been 
experimentally  recorded . 

Dynamic  elasto-p lastic  finite  element  analyses  of  the  specimens  have  been 
performed  It  is  shown  how  the  plastic  constitutive  laws  have  been  determined  by 
comparing  finite  element  simulations  end  experimental  results,  in  spite  of  stress  and 
strain  gradients  that  occur  in  smooth  bars  under  dynamic  loading  .  The  results  are 
compared  to  those  obtained  under  dynamic  compression  and  under  static  solicitation  both 
in  tension  and  compression  . 

The  stress  and  strain  distribution  in  the  round  notched  bars  under  dynamic  loading 
have  been  obtained  by  finite  element  calculations  as  a  function  of  time  for  3  notch 
geometries  leading  to  different  stress  triaxi’ity  ratio  distributions  .  Prom  these  numerical 
results  and  experimental  results,  multlaxtal  damage  laws  similar  to  Rice  and  Tracey's 
ductile  growth  of  voids  model ,  have  been  shown  to  applicable  Rupture  criteria,  based  on 
critical  void  growths,  have  also  been  determined  .  Differences  between  experimental  and 
predicted  loads  and  meen  diametrical  strains  at  rupture  have  been  interpreted  in  terms  of 
the  existence  of  a  critical  threshold  plastic  strain  for  void  initiation  wicli  has  not  been 
taken  into  account  in  the  model . 

The  experimental  procedures  followed  to  determine  the  fracture  toughness  from 
tests  performed  on  axisymmetrtcally  cracked  specimens  have  been  justified  on  the  basis  of 
the  finite  element  simulations  ,  Experimental  and  calculated  variations  with  time  of  the 
load  In  the  reduced  section  and  the  crack  mouth  opening  displacement  have  been  shown  to 
be  in  good  agreement  .  Fracture  toughness  have  been  found  to  be  well  predicted  by 
integrating  the  damage  low  identified  on  notched  bars,  and  by  assuming  that  cracks  Initiate 
as  soon  as  the  rupture  criterion  is  reached  over  a  critical  microstructural  distance  from 
the  crack  tip .  The  model  can  explain  the  increase  of  fracture  toughness  with  loading  rate 
in  35  NCD  1 6  steel  for  which  an  increase  of  the  flow  stress  with  strain  rate  is  observed, 
and  for  which  dynamic  rupture  appears  to  proceed  by  ductile  void  growth  mechanisms  and 
not  by  cleavage . 
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In  recent  years  considerable  attention  hae  been  focusBed  on  to  the 
behaviour  of  cellular  materials.  Materials  with  u  regular  honeycomb 
structure  have  been  studied  under  conditions  of  transverse  loading  (e.g. 
Klintworth  and  Stronge  (1988)  and  Gibson  et  al  (1982))  and  various  features 
of  their  modes  of  deformation  have  been  discussed  in  detail.  In  particular 
the  tendency  of  auch  materials  to  fail  locally  at  well  defined  stress  levels 
has  been  identified  as  a  major  characteriatic.  Similar  behaviour  is  to  be 
found  in  many  other  man-made  (e.g.  polymeric  foams)  and  naturally 
occurring  (e.g.  wood)  cellular  materials.  When  cylindrical  specimens  of 
cellular  materials  are  crushed,  the  deformation  results  from  the  growth  of 
these  localized  failure  zones,  the  deformation  taking  place  at  a  constant  or 
gradually  increasing  stress  (the  crushing  stress).  For  a  constrained  or 
low  density  material,  thit  continues  until  the  length  of  the  specimen  has 
been  engulfed  by  the  deformation  fronts  at  which  point  rapid  hardening 
occurs  as  the  crushed  ct  11b  are  themselves  compressed  in  a  consolidation 
phase.  For  an  unconstrained  higher  density  specimen  the  deformation 
ceaseB  when  the  specimen  fractures.  Fig.  1  shows  typical  load-deformation 
behaviour  for  samples  of  Yellow  Pine  and  Oak  tested  parallel  to  the  grain 
and  perpendicular  to  it. 

The  characteristic  shape  of  the 
load/displacement  (stress/mean  strain)  300 
curves  makes  cellular  materials 
attractive  as  potential  materials 
for  use  in  impact  energy  absorption.  130 
This  has  been  discussed  by  Gibson 
and  Ashby  (1988)  in  their  interesting 
new  text  on  the  topic  of  cellular  100 

materials.  Wood  has  traditionally 
been  used  for  energy  absorption 
and  containment  and  applications  33 

can  be  found  ranging  from  old  wooden 
warships  (Johnson  1986a, b)  to  present 
day  designs  for  nuclear  transport  0 

flasks.  The  use  of  wood  in  such  0  1  0  3  *  3  6  7  *  *  is  tt  u  is  i« 

dynamic  loading  applications  and  Displacement  mm 

more  especially  the  need  to  model 
impact  events  involving  wood 

presupposes  a  quantitative  understanding  Fig.1  Guasi-static  load-deflection 
of  the  response  of  wood  under  high  rate  curves  for  cylindrical  wood  specimens 
loading  conditions.  To  the  author's  45  mm  dla.  x  18  mm  thick, 

knowledge  no  systematic  examination  of 
the  dynamic  behaviour  of  wood  has  appeared 
in  the  literature. 

A  Bample  of  the  data  for  two  woods,  Yellow  Pine  and  American  Oak, 
produced  in  an  extensive  set  of  impact  (Peng  1987)  and  penetration  tests 
will  be  presented.  The  first  set  of  tests  involves  an  apparatus  in  which 
cylindrical  specimens  of  wood  were  attached  to  an  aluminium  backing  disc 
and  push  rod  and  fired  from  a  cartridge  gun  against  the  end  of  a  load  cell 


at  velocities  of  around  100  ms-1.  The  behaviours  of  specimens  loaded 
parallel  to  the  grain  and  across  the  grain  were  investigated.  The 
deformation  of  some  of  the  specimens  was  filmed  and  the  initial  localised 
nature  of  the  deformation  observed.  The  impact  loads  were  measured  using 
a  long,  strain-gauged  load  cell  which  measures  force  puleos  in  the  same 
way  as  in  a  Hopkinson  bar  apparatus.  Enhancement  of  the  initial  crush 
stress  by  factors  of  up  to  2.5  for  loading  parallel  to  the  grain  and  5.0  for 
loading  acrosB  the  grain  were  observed. 


The  penetration  of  flat  ended 
cylindrical  projectiles  fired  normally 
into  a  block  of  wood  produces  distinctive 
deformation  patterns  which  depend 
primarily  on  the  orientation  of  the  grain. 
The  extremes  are  sketched  in  Fig.  2. 

This  shows  a  localized,  heavily  crushed 
plug  beneath  the  indenter  when  the  grain 
Is  parallel  to  the  direction  of 
penetration  and  a  more  extensive 
frustum-shaped  cavity  behind  a  region  in 
which  the  fibres  have  been  crushed 
transversely  and  bent  when  the  loading 
is  across  the  grain.  As  wall  as  describing 
the  change  in  the  nature  of  deformed 
region  as  the  impact  velocity  is  increased, 
penetration  force  traces  are  presented 
which  were  acquired  by  installing  an 
accelerometer  in  the  cylindrical  - 

projectile.  At  impact  velocities  of 
approximately  50  m/s  the  initial 
penetration  stress  increases  by  a 
factor  of  four  or  more  compared  with 
the  quasi-static  value  for  a  short 
period  following  contact  before  falling 
to  a  lower  plateau  which  is  about  twice 
the  quasi-static  penetration  force. 

In  both  types  of  test  the  localised 
nature  of  the  deformation  process  is 
evident  and  will  be  discussed. 
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Fig. 2  Deformation  modes  for  Yell 
Pine  following  cylindrical  penet 
(a)  along  the  grains  (b)  perpend 
to  the  grains. 
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Abstract  -  A  major  advance  in  the  toughening  of  ceramics  has  been  achieved  by  means 
of  ceramic  fiber  reinforcement.  The  principal  mechanism  by  which  benefit  is  realized 
is  through  pullout  of  fibers,  whose  frictional  resistance  creates  sufficient  crack 
closure  to  significantly  raise  the  fracture  toughness.  However,  all  studies  to  date 
have  involved  relatively  slow  rates  of  loading;  there  is  some  question  as  to  whether 
this  mechanism  w'll  obtain  at  rapid,  or  impact,  loading  rates.  This  paper  reports  on 
the  results  of  tests  aimed  at  evaluating  the  dynamic  versus  static  fracture  behavior  of 
SIC-fiber  reinforced  glass-ceramic  composite. 

The  dynamic  fracture  tests  have  been  performed  using  an  adaptation  to  brittle 
materials  of  a  recently  developed  dynamic  fracture  experiment  [1].  Hopkinson 
pressure  bars  preloaded  in  tension  are  utilized  to  load  simultaneously  two  compact 
type  specimens  with  notch  of  different  radius.  The  experiment  is  characterized  by 
well-defined  specimen  boundary  loading  conditions,  and  by  stress  intensity  rates  {tp 
on  the  order  of  10°  MPe/m  s  .  Crack  velocities  are  measured  by  using  a  crack  gage 
technique. 

Static  fracture  stresses  are  compared  with  the  stresses  required  for  dynamic  crack 
initiation,  and  results  are  interpreted  in  terms  of  fracture  mode.  It  will  be  shown  that 
even  under  dynamic  loading  conditions,  the  influence  of  the  weak  axially-oriented 
fiber-matrix  Interfaces  manifest  themselves  in  the  fracture  path.  Furthermore,  the 
observed  multiple  fracture  locations  appear  to  be  related  to  the  loading  mode  involved 
In  the  compact  type  specimen  geometry. 


[1]  H.  COUQUE,  S.J.  HUDAK,  JR.  and  U.S.  LINDHOLM,  J.  Phys.,  Colloque  C3,  Sup. 
n°9,  T49,  C3-347,  1988. 


